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PREFACE 


New England Division (NED) participated in Landsat Follow-on 
Investigation No. 22510 to demonstrate the construction and oper- 
ation of a regional data downlink and to explore operational uses 
of Landsat imagery for hydrology. It was shown that a relatively 
simple, low cost, automatic ground receive station could be con- 
structed and operated for data collection; and a study was made 
of matters related to the operation of the station, such as data 
reduction and management, personnel requirements, daily operating 
procedures, reliability, maintenance, and costs. A separate part 
of the investigation was carried out by the U.S. Cold Regions 
Research and Engineering Laboratory (CRREL) to explore computer 
analysis of digital images for extraction of hydrologic informa- 
tion and to develop techniques that will lead to the quantifica- 
tion of the water content of snow. 

Participation in the Landsat Follow-on was stimulated by several 
factors. The original investigation, documented in NED’s Landsat-1 
Final Report (Cooper et al, 1975), suggested several avenues of 
further research for making imagery analysis and data acquisition 
operational for agencies such as the Corps of Engineers. Oppor- 
tunities and resources were made available in the present investi- 
gation to gain experience in the technology of data collection by 
satellite and to make the comparisons and judgments necessary for 
developing fully operational systems. Alternative data collection 
systems have been examined, and based upon critical comparisons of 
cost and reliability, satellite methods are rated favorably. It 
has been possible to define precisely the range of downlink per- 
formance that users can expect with respect to maintenance, reli- 
ability, and personnel requirements. In the 1972-1974 Landsat 
investigation, the reliability of the overall data collection sys- 
tem (DCS) was proven, and the follow-on has demonstrated that an 
automatic ground receive station can also perform with high 
reliability. 

The overall program was administered and directed by Saul Cooper, 
Principal Investigator. The portion of this report dealing with 
data collection and the ground receive station was prepared by 
Timothy D. Buckelew, Hydrologist, New England Division. Much 
credit and appreciation are due to many others who have assisted 
in various phases of the data collection activities, but especially 
to Mr. Paul Hetu (NED) who has been responsible for the management 
and field installation of equipment; to student assistants 
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Richard Colburn and John Burnham; to several offices of the U.S. 
Geological Survey for their cooperation and sharing of equipment 
in stream gaging stations; and to NASA for supporting such a pro- 
gressive investigation. Special thanks go to Jacqueline Izzi for 
typing this manuscript. 

The portion of the study dealing with imagery was under the 
direction of Dr. Harlan McKim (CRREL) and was prepared by 
Carolyn J. Merry, Research Geologist, Earth Sciences Branch, 
Research Division. 

The authors also express their appreciation to the following per- 
sonnel who assisted with CRREL 1 s investigation: Blanchard Pratt 

and Richard S. Guyer for sensor interface development for the 
Landsat data collection platforms; to Gregor E. Fellers for devel- 
opment of the computer programs used to decode the data from the 
DCP f s; to Eleanor Huke for her assistance in layout and arrange- 
ment of the computer displays and figures; to Dr. Stephen G. Ungar 
(NASA Goddard Institute for Space Studies) for development of the 
computer algorithms used in the analysis of the Landsat digital 
data; to Thomas L. Marlar and David A. Gaskin for field support; 
to Lawrence W. Gatto and Dr. Samuel Golbeck for technical review 
of sections of the snow cover analysis; to Ronald T. Atkins for 
technical review of the DCP sensor interface description; to 
Michael A. Bilello and Harold O f Brien for technical review of the 
snow cover analysis section. 

The contents of this report are not to be used for advertising or 
promotional purposes. Citation of trade names does not constitute 
an official endorsement or approval of the use of such commercial 
products. 
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1 . 0 INTRODUCTION 


This report documents progress made by New England Division, 
Corps of Engineers, in utilizing Landsat data and imagery products 
in watershed management, and it includes a description of sensors 
and interfaces developed for measuring water related parameters. 

The project was the result of the continuing need for reliable, 
timely, hydrologic data, and it evolved as a sequel to an earlier 
investigation into the feasibility of using satellite data for 
reservoir management (Cooper et al, 1975). 

1 . 1 Organization of this Report 

Sections 1 through 5 and Appendix A cover data collection 
activities, while sections 6 through 10 and Appendix B concern com- 
puter analysis of Landsat imagery and development of sensors. 

What follows in this introduction is a precis of the earlier report 
giving the major results and recommendations which led to the 
present follow-on. 

1 . 2 Earlier Work with Landsat 


Based on the 1972-1974 experience with a network of 26 
data collection platforms, NED found real time data collection by 
orbiting satellite to be both reliable and feasible. For large 
scale systems of national scope, it is possible to design orbiting 
satellite systems that are more flexible, easily maintained and 
less expensive than conventional ground-based means. Even though 
the frequency of Landsat appearances (four to six times daily) is 
not suited to watershed management needs in all areas, it should 
be understood that the Landsat system is experimental; and one of 
the principal objectives is to test the feasibility of data col- 
lection by orbiting satellite. 

In 1974 NED endorsed the institution of a satellite data 
collection system on a Corps-wide basis or a nation-wide system 
with sharing among other Federal and State agencies. A further 
recommendation was that any operational satellite configuration 
should include ground receive stations at major user locales for 
direct receipt of critical data, rather than the relay of data by 
land lines from a single station. 

Experience with imagery in the former study indicated 
that .Landsat photos may be enlarged about five times, or to a 
scale of 1:200,000 which is sufficient for defining gross feature 
patterns, such as the depiction of floods from the larger rivers 
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in New England. This was judged to be only marginally useful for 
reservoir regulation purposes. Ice is readily detectable on the 
imagery as is the boundary between ice and open water. This was 
found useful for monitoring developing ice cover, especially over 
remote areas. Winter snow cover patterns are readily obtainable 
with excellent accuracy from Landsat imagery; however, photographic 
imagery provides only snow location, not water equivalent which is 
the important parameter from an operational viewpoint. Landsat 
imagery appears able to distinguish areas previously but no longer 
flooded, for periods of several months after flood recession. 

NED concluded in the 1972-1974 report that the coordinated 
use of all data available to an operational reservoir control center 
should include the interaction between real time imagery and point 
data sources, such as ground truth obtained by the Landsat DCS. 
Before this interaction could become a reality it was considered 
necessary to provide some means of real time relay of Landsat 
imagery to an operational control center. In addition, the report 
recommended the continued development of an interactive system to 
better utilize computer compatible tapes to depict changes in 
hydrologic conditions. 

2.° HISTORICAL BACKGROUND 

2.1 Corps* Mission 

NED’s involvement in space technology for imagery and data 
collection is tied to the mission of the Corps and more subtly to 
the manner in which New England has developed in the last century. 
Since the Industrial Revolution in the 1800’s, the rivers of New 
England have been developed to supply water for power and transpor- 
tation. As new means of transportation became more economical, 
both railroad and highway systems were built along the banks of the 
rivers to serve the expanding needs of the industrial, commercial 
and urban centers. Structures, such as buildings, roads, bridges 
and dams have restricted floodways to such an extent that consider- 
able property and environmental damages have occurred even during 
only moderate floods. Notable floods of November 1927, March 1936, 
September 1938 and August 1955 have demonstrated the need for flood 
control to prevent these natural catastrophes. 

At the direction of Congress, the U.S. Army Corps of 
Engineers developed a comprehensive plan of flood protection for 
each river basin. The protective works recommended generally con- 
sist of a combination of channel improvements, dikes or floodwalls 
at major damage centers augmented by upstream flood control reser- 
voirs. Many of these reservoirs contain additional storage 
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reserved for other uses such as water supply, conservation and 
recreation. The Corps has built 35 flood control reservoirs, 44 
local protection projects and four hurricane barriers in New 
England at a total investment of $350 million. 

2.2 Data Collection at NED 


To gain maximum benefits from this comprehensive protec- 
tion system, the New England Division requires hydrologic data 
such as river, reservoir and tidal levels; wind velocity; baro- 
metric pressure; and precipitation. 

In the past this data was collected by field observations 
and relayed via telephone or voice radio. For a scan of New 
England it took several hours to compile and assess the data in 
this manner. As the flood control system expanded and the need 
for timely and reliable information increased, the Corps pioneered 
new methods of data collection. 

In 1970 the Automatic Hydrologic Radio Reporting Network 
(AHRRN) was placed in operation. This ground-based radio relay 
system consists of 41 remote reporting stations, and a central con- 
trol at Division Headquarters in Waltham, Massachusetts. This 
network, with the assistance of a computer, automatically collects 
and reduces the real time data that is essential for reservoir 
regulation. The remote reporting stations of this system are 
strategically located in five major river basins and at key coastal 
positions, with each contributing to a comprehensive hydrologic 
picture. 

2.3 Entry into Satellite Data Collection 

In June 1972 the Corps contracted with NASA for an exper- 
iment to study the feasibility of using the Earth Resources Tech- 
nology Satellite for collection of environmental data from data 
collection platforms (DCP ? s). Of the 26 B€P f s installed in loca- 
tions throughout New England, many are situated at existing U.S. 
Geological Survey gaging stations. Between July 1972 and 
September 1975, uandsat relayed river stage, precipitation, and 
water quality data from DCP f s to the Goddard Space Flight Center, 
whence it was sent to the New England Division within one or two 
hours via teletype. 

As early as 1973 NED personnel envisioned a ground receive 
station, because it was apparent that real time data was needed 
without danger of disruption from regional power or telephone 
failures. Table 1 lists chronologically the major steps related 
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to development of the ground receive station from 1973. Included 
in the table are significant events associated with the construc- 
tion of the ground receive station, such as presentations, meetings 
and a special briefing of the Chief of Engineers in January 1976 
which led to a Corps-wide examination of the use of satellite data 
collection. 

The present investigation began in 1975; and by late 1975 
NED had constructed an inexpensive, semiautomatic, and easily 
maintained ground receive station as a follow-up to its o7:iginal 
study. The Division is now able to receive hydrometeorological 
data from data collection platforms in the field directly at its 
headquarters in Waltham, Massachusetts with no time delays. The 
satellite tracking system operates unattended automatically at all 
times, with a computer controlling all processes. 

2.4 Imagery Experience 

From 1972 to the present, NED, in cooperation with the 
University of Connecticut and the Cold Regions Research and Engin- 
eering Laboratory, studied computer compatible tapes of Landsat 
scenes for potential hydrologic applications. In the Landsat-1 
investigation it was concluded that gross outlines of waterbodies 
and snow or ice could be distinguished, but little advantage would 
be gained in watershed management without refined interactive 
computer techniques. Since that time Landsat photographs and 
taped scenes have been collected for selected areas of New England 
for hydrologic analysis. Work is continuing on development of 
computer techniques and regression analyses to relate multi- 
spectral "signatures" with systems involving snow water content , 
ground cover types, and varying slopes and elevations. 


TABLE 1 


CHRONOLOGY OF THE LANDS AT FOLLOW-ON 
INVESTIGATION AT NEW ENGLAND DIVISION 


1973 

September - 
December 

May 

1974 

February 

August 

September 

October 

November 

1975 
January 

February 

March 

April 

April - May 


Preliminary discussions among personnel of NED, GSFC, 
NASA Wallops. 

Satellite Workshop at Wallops Station, Virginia. 

NASA/NED exploratory meeting at Waltham. 

Funds transferred from OCE to NASA. 

Final system specifications drafted by NASA for ground 
receive station. 

Ground receive station site preparation 
Inspection of DCS demodulator/decoder at Wallops. 
Conference at GSFC on satellite tracking software. 


Technical proposal on ground receive station submitted 
by Scientific Atlanta. 

Landsat Follow-on contract signed. 

Coordination meeting among personnel of NED, CRREL 
and University of Connecticut. 

Meeting at Waltham with NED, GSFC, Wallops and 
Scientific Atlanta personnel present. 

Started acquiring 9-day imagery coverage of New England. 

Data sharing between NED and Saint John (New Brunswick) 
River Flood Forecast Center. 
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1975 (cont.) 

Hay Coordination meeting among personnel of NED, CRREL and 

University of Connecticut. 

June Type II Report No. 1. 

Satellite tracking programs debugged. 

July Acceptance trial of antenna and pedestal at Scientific- 

Atlanta in Atlanta, Georgia. 

September Type II Report No. 2. 

Arrival of antenna in Waltham. 

October Antenna installation. 

Presentations on data collection and imagery at Tenth 
International Symposium on Remote Sensing of the 
Environment, Ann Arbor, Michigan. 

Presentation before Atlantic Fisheries Biologists 
meeting, Newagen , Maine. 

Presentation before International Telemetering Con^ 
ference (Silver Springs, Maryland) on ground receive 
station and imagery analysis. . »-■ 

November Initial tracking of Landsat by ground receive station. 

Presentation before American Water Resources Conference 
at Baton Rouge, Louisiana. 

Tulsa, Oklahoma workshop on convertible data collec- 
tion platforms. 

December Type II Report No, 3. 

1976 

January Briefing of General Gribble, Chief of Engineers, on 

satellite data collection. 

Meeting at CRREL on cooperative remote sensing pro- 
grams, among personnel of NASA GISS, and CRREL. 


1976 (cent.] 
February 

March 

April 

May 

June 

July 

September 

October 

November 
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Landsat DCS demonstration held at Boston USGS Regional 
office. 

DCS data sharing with Saint John (New Brunswick) River 
Flood Forecast Center during spring runoff. 

Type II Report No. 4. 

Coordination meeting at CRREL on CCT algorithm. 

Arrival of LaBarge convertible DCP's at NED. 

Consultation with NED/CRREL on use of Landsat imagery 
for detection of red tide by Mr. Jerry McCall of the 
Massachusetts Department of Environmental Quality 
Engineering. 

Installation of time code generator in tracking system. 
Type II Report No. 5. 

Meeting at Sugarloaf Mountain to plan emplacement of 
thermocouple cable. 

Video taping of DCS equipment for public information. 
Type II Report No. 6. 

Demonstration DCP assembled at ground receive station. 

Snow pillows installed at NED and northern Maine. 

Thermocouple chain interfaced with DCP and installed 
on Sugarloaf Mountain, Maine. 

NASA briefing before discipline specialists. 

CRREL /GISS meeting to plan ongoing Landsat digital 
analysis. 

Informal Earth Resources Program Review at GISS, New 
York, New York. 

ORIGINAL PAGE* 

OF POOR QUALITY 
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February Presentation on analysis of water equivalent of snow 
using Landsat imagery, before Eastern Snow Conference 
Belleville, Ontario. 

April Presentation at Meteorological Satellite Workshop on 

data collection systems, White Sands, New Mexico. 

May OCE briefing and survey conducted at NED jointly with 

NASA, NESS and OKI. 

Presentation on permafrost in New England at a seminar, j 

University of Maine, Institute of Quaternary Studies, 

Department of Geological Sciences. 



3.0 


SCOPE AND OBJECTIVES OF THIS INVESTIGATION 


3. 1 Evaluation of the DCS 


The first principal objective of this investigation was to 
evaluate the effectiveness of the DCS in aiding watershed manage- 
ment functions as compared to other conventional means. The study 
demonstrates the degree of utility of the Landsat data collection 
system with respect to reservoir regulation. A direct downlink was 
installed at Waltham as a cooperative effort between NASA and the 
Corps of Engineers. Areas that were explored in the management of 
the downlink were computer programming; data reduction, storage and 
retrievel; maintenance, downtime, service and costs; and daily 
operating procedures and personnel requirements. Comparisons of 
costs and operating procedures were to be drawn among alternative 
data collection systems such as GOES and ground-based radio. 

3. 2 Evaluation of Imagery 

The second principal objective was to develop methods of 
using Landsat imagery information to assist in the planning, manage 
ment and operation of NED water resource projects. To this end, 
several tasks were undertaken. Imagery was analyzed to evaluate 
mapping accuracy of the areal extent of snow. A preliminary rela- 
tionship was determined between the water equivalent of a snowpack 
and radiance measured by Landsat. Imagery was also used to delin- 
eate wetlands and floodwaters in New England. Along with the two 
main objectives, it was possible to develop sensors and interfaces 
for obtaining environmental real time data by the Landsat DCS. 

4.0 LANDSAT DATA NETWORK AT NED 

4 . 1 Initial Steps 

The present DCS study into operational uses of Landsat 
was an extension of the 1972-1974 feasibility investigation, and 
plans for the follow-on investigation were underway in 1973. It 
was' recommended that surplus equipment, compatible with our needs, 
be used to construct a regional downlink at Waltham . NASA per- 
sonnel provided technical support in planning and integrating 
components; and a contract was issued to Scientif ic-Atlanta in 
Atlanta, Georgia, for refurbishing the surplus equipment and fur- 
nishing several new components. NED prepared the site for the 
antenna, furnished the computer equipment for controlling the sys- 
tem, and was responsible for software development. 


9 


Installation 




The system was installed and operating in late 1975; and 
after testing and shakedown, NED examined all facets related to 
operating and maintaining a downlink. The most significant topics 
examined include: 

a. DCP development 

b. personnel requirements 

c. malfunctions, downtime and repair costs 

d. comparisons of downlink to conventional methods 

e. computer programs for orbit prediction and tracking 

f. data reduction and storage 

g. daily operating procedures 

Each of these topics is covered in the sections which 
follow; and an in-depth view of all computer programs, data files, 
special devices, and operating procedures is given in Appendix A, 

4.3 System Design and Operation 


The Landsat tracking system at NED integrates a set of 
about 20 programs or subroutines (software) , about ten disc data 
files, and several pieces of equipment (hardware) . The hardware 
components are listed in table 2 and depicted in figures 1 and 2. 
Software and data files which were designed at NED are described 
in detail in Appendix A. 

4.4 Network Description 

The overall data collection system (DCS) comprises the 
satellite, one or more ground receive stations, and many remote, 
automatic data collection platforms, which are equipped by users 
to sample local environmental conditions (NASA, 1976a). Up to 
26 platforms were deployed by NED at any time during the course of 
this investigation. 

4.5 Deployment of DCP f s 


4.5.1 


Site Selection 


NED T s data collection platforms were deployed to serve 
various functions while testing the feasibility of satellite data 
collection systems. Sites were established for: 


- field test purposes 

- expansion to areas outside the AHRRN coverage 

- demonstrations 


a . 
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TABLE 2 


COMPONENTS OF LANDSAT TRACKING STATION AT NED 


Item 

Number 

Manufacturer 

Antenna , Parabolic , 1 5 1 

custom made 

Scientific-Atlanta 

RF Feed and Waveguide Assembly 

custom made 

ii 

Tracking Pedestal 

Model 3203 

it 

Manual Command Unit 

Model 3732 

ii 

Servo Control Unit 

Model 361 5 A 

ii 

DC Amplifier 

Model 3641 

it 

SCR Power Amplifier 

Model 3635 

?f 

Digital Comparator 

Model 1848 

n 

Punched Tape Programmer 

Model 3823 

it 

Tape Reader 

Model RRS6300 

Rem ex 

Parametric Amplifier 

Model SCP-290 

Scientific Communi- 
cations, Inc. 

Air Dryer 

Model 550 

Puregas 

Digital Synchro Display 

Model 1842 

Scientific-Atlanta 

Receiver, Basic Unit 

Model 410A 

M 

UHF RF Tuner 

Model 423 

it 

PM Demodulator, Wide Angle 

Model 444A 

if 

IF Filter 

Series 430 

»t 

Spectrum Display Unit 

Series 450 

ii 

Signal Generator, S-band Test 

Model 7100 

Microdyne Corp. 

Time Code Generator 

Model 9100A 

Datum, Inc. 

Demodulator /Decoder , Data 

custom made 

NASA, Wallops Station 

Terminal, Cathode Ray Tube 

Model 4014-1 

Tektronix 

Hard Copy Unit 

Model 4631 

it 

Graphics Tablet 

Model 4953 

tt 

Minicomputer 

Model 1220 

Data General Corp. 

Disc Unit, Magnetic 

Model 4057 

ft 
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The primary criterion for DCP deployment was the need for 
continuous real time data collection in sites that supported NED f s 
water control mission. These comprised the majority of our DCP’s 
and were usually interfaced with stream or precipitation gages (see 
sitelist, figure 3). In some cases sites were chosen for their 
remoteness, climatic harshness, or on the other hand, for their 
accessibility. DCP f s were co-located at existing AHRRN stations, 
such as Plymouth and Goffs Falls, New Hampshire and Hartford, 
Connecticut to compare the hydrologic data* 

Several DCP’s extended or enhanced our coverage of New 
England, such as Cornish, North Anson, West Enfield, Fort Kent, 
Ninemile Bridge, and Allagash Falls in Maine; and Cranston and 
Forestdale, Rhode Island. DCP r s at these sites provided more com- 
prehensive real time data from New England rivers than was prac- 
tical with other methods at our disposal. 

It should be noted here that DCP’s were not placed with 
the strategic aim of integration with the Merrimack Hydrologic 
Model which has been under development by NED and the Hydrologic 
Engineering Center in Davis, California. In the early assessment 
stage of the investigation this integration was considered to 
be beyond our reach with regard to the time, manpower, and funds 
available. Implementation of the HEC model is being pursued using 
the AHERN. 


A final category of DCP site selection includes those set 
up as demonstrations and tests. These include: 

a. DCP’s interfaced with snow pillows in the winters of 
1975-1976 and 1976-1977. 

b. DCP’s interfaced with meteorological stations in 
North Dakota by CRREL. 

c. A DCP interfaced with a thermocouple chain in Maine 

by CRREL. 

d. A demonstration DCP at NED to show sensors, inter- 
faces, platforms, antennas, and typical cabling. 

e. A DCP at Manchester, Connecticut, Nature Center main- 
tained by the Hartford office of the U.S. Geological Survey, as a 
public demonstration. 
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LANDSAT-8 - DCP INFORMATION SHEE T 
U.S. ARMY CORPS OF ENGINEERS, NEU ENGLAND DIVISION 


18 OCTOBER 1577 


DCP 

NO. 

STATION NAME 

PARA- 
METER CS 

>* 

LAT 

LONG 

7912 

ST. JOHN RIUER AT NINEHILE BRIDGE, PIE . 

RS UES 

46 

42 

99 

69 

42 

59 

7273 

ST. JOHN RIUER AT FORT KENT. ME. 

RS 

47 

15 

27 

68 

35 

35 

7394 

SHIELDS BROOK NEAR ST a POMPHILE, P.G. 

RS 

46 

56 

39 

69 

38 

58 

7119 

ALLAGASH FALLS* ME. AT 

PENOBSCOT RIUER AT WEST ENFIELD, ME. 

,GT, UES 

46 

57 

95 

69 

11 

43 

79? 1 

RS 

45 

14 

12 

68 

38 

56 

7278 

CARABAS5ETT RIUER NEAR NORTH ANSON, ME. 

RS 

44 

52 

99 

69 

57 

29 

7356 

SACO RIUER AT CORNISH, ME. 

RS 

43 

48 

35 

79 

<6 

53 

7127 

SOUTH MOUNTAIN, N.H. 

PEMIGEUASSET RIUER AT PLVMOUTH, N.H. 

P 

42 

58 

59 

71 

35 

21 

7291 

RS 

43 

45 

33 

71 

41 

19 

7921 

MERRIMACK RIUER NEAR GOFFS FALLS, N.H. 

RS 

42 

56 

54 

71 

27 

52 

7246 

UACHUSETT MOUNTAIN, MA. 

P 

42 

29 

24 

71 

53 

15 

6963 

IPSUICH RIUER NEAR IPSUICH, MA. (1) 

RS 

42 

39 

35 

79 

53 

39 

7271 

NORTH NASHUA RIUER AT FITCHBURG, MA. 

RS 

42 

34 

34 

71 

47 

19 

7197 

NED, WALTHAM, MA. (LABARGE) 

P 

42 

23 

46 

71 

12 

56 

7325 

NED, WALTHAM, MA. 

UES 

42 

23 

46 

71 

12 

56 

7345 

PAUTUXET RIUER AT CRANSTON, R*I, 

RS 

41 

45 

93 

71 

26 

44 

7254 

CONNECTICUT RIUER AT HARTFORD, CT. 

RS 

41 

46 

19 

72 

49 

94 

7355 

CONNECTICUT RIUER NEAR MIDDLETOWN, CT. 
PORTER BROOK NEAR MANCHESTER, CT. (2) 

RS 

41 

33 

49 

72 

36 

45 

7296 

RS 

41 

45 

55 

78 30 

12 


6316,7042,7101 (3 )RL,AT,GST,GT.UP 

7103-7105,7107-7120.7125 ASSIGNED TO CRREL'S C'DCP'S 
7147,7220,7142,7207,7171,7233,7242,7335,7010 SPARES 


P - PRECIPITATION AT 

UES - UATER EQUIVALENT GST 

OF SNOUPACK 

RS - RIVER STAGE GT 

RL - RESERVOIR LEVEL UP 

UQ - UATER QUALITV PV 

(TEMPERATURE, T 

CONDUCTIVITY, 

PH AND DISSOLVED 
OXYGEN ) 


AIR TEMPERATURE ( S ) 
GROUND SURFACE 
TEMPERATURE 
GROUND TEMPERATURES > 
UIND PASSAGE 
PARAMETERS VARIABLE 
TEST SET 


U) DCP OUNED BY U.S. GEOLOGICAL SURVEY, BOSTON. MA. 

(2) DCP ON LOAN TO U.S. GEOLOGICAL SURVEY, HARTFORD. CT.- 
-ON DEMONSTRATION AT THE MANCHESTER NATURE CENTER 

(3) DCP ON LOAN TO U.S. ARMY COLD REGIONS RESEARCH AND 

ENGINEERING LAB, HANOUER, N.H. 

Figure 3. Data Collection Platform Sit.liat 
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f. Water quality monitors interfaced with DCP's at 
Fitchburg, Chicopee, West Springfield and Webster, Massachusetts. 

4.5.2 Interfacing Sensors and DCP's 

The initial translation of river stage to a transmittable 
signal is often performed by a float well or manometer-type gage. 

In a float well configuration, the level of a float on the water 
surface is followed by a taut vertical tape which runs over a pulley 
on a shaft. The shaft is coupled to an analog to digital recorder 
(ADR) which translates the shaft rotation to a digitized signal. 
Either a "telekit" or a semiconductor memory unit latches and holds 
the latest value for the DCP to read. 

A manometer-type gage functions by sensing the amount of 
(nitrogen) gas pressure needed to force bubbles down through a 
tube into the water body to the datum level. A moveable mercury 
well is moved by a motor to balance the pressure of the gas. The 
motor shaft motion provides the analog of water stage. From this 
point the transmission path is similar to the float well config- 
uration. 


Rainfall is usually measured with weighing or tipping 
bucket rain gages that are equipped with ADR's. A gage records 
accumulated rainfall in the range from 0 to 99.99 inches. The 
data pathway to a DCP is similar to that of river stage. Rainfall 
rates may be computed from successive values and elapsed time 
between readings. 

Landsat DCP's allow sensor input of 64 parallel bits, 64 
serial digital bits, eight analog inputs, or combinations of these 
formats. Various brands of DCP's offer differing amounts of versa- 
tility and control of external devices. Users must configure their 
sensor equipment and cabling to conform to the data formats 
accepted by the DCP. A scheme using binary coded decimal (BCD) 
bit assignments was adopted in most applications at NED. By this 
method a four-digit decimal number representing the measured 
parameter is encoded in 16 bits divided into four groups of four 
bits each. It is convenient to handle such a group in its octal 
(base-8) form. Sixteen bits of information can be represented and 
conveniently handled as two octal triplets. For example, the 
octal triplets 237 373 represent the 16-bit pattern 10 Oil 111 11 
111 Oil. 

Each position of this pattern can be assigned the following 
values, respectively: 

1 2 4 8 10 20 40 80 .01 .02 .04 .08 .1 .2 .4 .8 
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During encoding at the DCP, the bit values are ones - 
complemented so that a 0 in the binary form stands for a "high” or 
presence of the decimal value, and 1 stands for "low" or absence 
of the decimal value. By these conventions the sample bit pattern 
above represents 2 + 4 + .2, or the decimal number 6,2. 

Several of these schemes or bit assignments were used at 
NED at various times and with various types of sensors and DCP f s. 

The problem of remembering which scheme is in use at each field 
location is handled by our computer using the program P3 and the 
disc file "INDX" (see Appendix A). Remembering over a period of 
months and years requires tables of station numbers, dates, and 
parameter types. 

Whether a binary 1 stands for "high" or M low" depends on 
the user's sensor equipment and whether it performs simple switch 
closure or has another sort of input driver. The value or signifi- 
cance of each bit is determined by users and is incorporated into 
wiring and pin assignments. Aay scheme adopted should be flexible 
enough to fit as many applications as possible. This will mini- 
mize the number of schemes and decrease the likelihood of misin- 
terpretations of valid data. 

4.6 NASA Support Data 

Several forms of support data were supplied by NASA during 
this investigation. DCS dats* was provided by teletype within a few 
hours and by punched card and computer printout within a week of 
acquisition. Photographic prints and transparencies of Landsat 
scenes of New England were also furnished. The scenes were processed 
by EROS Data Center of the USGS in South Dakota and were delivered 
usually from 4 to 8 weeks after acquisition by NASA. Image quality 
was high in most cases. At first, scenes were collected from all 
of New England, but in mid-1975 coverage was reduced to a study 
area in northern Maine, thereby increasing the number of passes 
that could be handled with available funding. Cloud cover toler- 
ance was also reduced to 30 percent or less. 

4 . 7 Personnel Requirements for Ground Receive Station 


Approximately one person-year was required for planning, 
site preparations, programming, and acquisition of equipment that 
was not furnished by NASA or the antenna contractor, Mr. Timothy 
Buckelew was the principal programmer in the DCS investigation. 

His background includes five years of high level language program- 
ming, mathematical training through calculus and analytical geometry 
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and exposure to sciences. Assembly level programming was needed 
and was learned on the job. At most times during program develop- 
ment, one or two college engineering students provided valuable 
assistance. 

Once the ground receive station components were fully 
integrated and operating according to specifications, the personnel 
requirements for it diminished to one half hour of attention per 
day to perform tasks such as: 

a. printing out data collected over night 

b. checking the time code generator against standard time 

c. entering orbital elements each week 

d. inserting or removing DCP directory entries 

e. retrieving sets of data as needed 

In addition, the service of a part-time technician was 
required for deployment and maintenance of a system of 30 DCP f s. 
This job includes antenna installation, wiring, interfacing, DCP 
initialization and checkout, battery maintenance, trouble-shooting, 
DCP replacement, documentation, etc. DCP T s that were faulty were 
shipped out for repair at a NASA-affiliated laboratory. One full- 
time technician could manage the entire system including downlink 
and DCP network provided there were access to support groups for 
service when malfunctions occur. 

4.8 Reliability of Ground Receive Station 

A study of DCP reception versus peak elevation angle to 
the satellite from Waltham was conducted over a 19-day period in 
April 1977. Reception is good for satellite passes having peak 
elevations down to 14°, For NED purposes a completely useful pass 
is one in which at least one message arrives per DCP. Usefulness 
diminishes quickly from 14° to 12° passes; and over half the net- 
work becomes silent for 10° to 5° passes. The last DCP’s to report 
are the three experimental platforms in North Dakota. 

The loss of messages as the satellite progresses to the 
west is caused by several factors: antenna gain patterns on DCP’s 

and the satellite; increase in range; and possibly the fact that 
the downlink antenna at Waltham begins to be obscured by buildings 
when elevation angles are less than 13°. The significant finding 
of this examination is that the DCS is highly useful for passes 
having peak elevations as low as 12°. 

The ideal configuration for a regional downlink is similar 
to the system tested during this experiment . The sizes of compo- 
nents used were ones available as surplus equipment. A 15-foot 
dish was installed at NED, but from discussions with other people 
in this field it was clear that a 12-foot dish would have been 
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adequate. Most other subsystems in the downlink are nearly ideal 
because the hardware and software were tailored to fit needs of 
this specific purpose. 

The weakest link in NED f s system was the S-band receiver. 
We lost many days of service when this component failed. It is 
recommended that two be purchased, if possible, for system re- 
dundancy. Margins of estimates in the RF link analysis may have 
been excessive, but it is hard to say precisely which items were 
overly conservative. Hardware present but not needed includes an 
extra position indicator (which came as surplus equipment) ; a 
spectrum display unit in the receiver,; and a programmable paper 
tape reader for input of azimuth and elevation angles. The tape 
reader would not be necessary if the pedestal were interfaced 
directly to the computer or the pedestal had auto tracking capa- 
bility. 


The performance of the ground receive system conformed 
to NASA’s specifications and far exceeded NED’s original expecta- 
tions. Bata reception quality (according to system specifications) 
had to be at least 90 percent error-free when the satellite is 
above elevation 25°. This degree of performance has been consis- 
tently demonstrated by the system. Moreover when examined in the 
context of collection of hydrologic data, the performance of the 
system including a set of DCP T s transmitting at 3 -minute intervals 
is much better than that which is indicated above. For hydrologic 
events, rates of parameter change are typically slow, so a high 
proportion of messages from one DCP in one satellite pass are 
identical. Receipt of several identical messages does not increase 
the amount of intelligence received; thus, if at least one message 
comes in during a pass, it still represents a high degree of success. 
(Upward and downward trends of a parameter can be calculated from 
comparison of individual DCP’s messages from successive passes.) 

DCS data furnished from_NASA by t elet ype w as us ed as a 
backup to our ground receive station and as a standard of compari- 
son for the station’s performance. NED’s receipt of DCP messages 
ranged from 80 to 90 percent that of NASA. This performance is 
considered excellent. Four messages per DCP per pass is the prac- 
tical maximum for NASA for high satellite passes. Receipt of five 
messages is rare. In the calculation of performance as a percent, 
there is a discrete downward jump as one message per DCP per pass 
is lost, causing a strong tendency toward 75 percent, or three out 
of four. Our 15-foot antenna may miss a DCP’s message, transmitted 
when the satellite is low, which NASA’s high gain 40 or 60-foot 
dishes will receive; but the performance of the 15-foot dish for 
higher segments of passes compares very well to larger dishes . 
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Putting the 80-90 percent performance in other words, slightly 
over one half the time we got the same number of messages as NASA 
and the rest of the time we missed one message. It was rare for 
NASA to receive a DCP that we failed to get at least once per 
pass, and vice versa . 

The sources of errors experienced were not checked, be- 
cause more sophisticated test equipment would have been needed. 
Reliance was placed on the decoder/demodulator/buffer subsystem to 
flag errors detected during convolutional decoding and not to flag 
valid messages as erroneous. After that, data were reduced and 
checked by the minicomputer for reasonableness (that is, to fall^ 
within limits). 'Finally data were scanned by the users for reason- 
ableness or the presence of discontinuities; an example can be seen 
in figure 4. The abrupt increase shown in this figure was judged 
to be valid because of the likelihood of "Ice runs 11 in April; and 
hourly values recorded by the USGS showed a smooth, rapid increase. 
Users should suspect nonvalid data if the discontinuity varies 
consistently by one bit value in the "binary coded decimal" coding 
scheme. 


A longitudinal study of the performance history of all our 
General Electric .DCP T s placed in the field from 1972-1976 shows 
that they are very reliable. Available data compiled from tech- 
nician’s records include the duty cycles (sometimes intermittent) 
of 23 DCP’s over a span of approximately four years for a total 
service history of 70.1 DCP-years. During this time several events 
could have forced a visit to the DGP site: initial installation, 

mechanical or electronic failure, preventive maintenance, battery 
change, or vandalism. For this study, final removal was not de- 
fined as a forced visit since removal was not urgent iii most cases. 
Over the 70.1 DCP-years there were about 95 forced visits for an 
average of 1.4 forced visits per DCP-year. In other words, slightly 
under nine months elapsed between two forced visits to the same DCP. 
Our early experience (or inexperience) with DCP’s is included in 
this analysis, so it may be assumed that for us or another exper- 
ienced group to start now with the knowledge already gained would 
involve far fewer DCP failures and on the average assure nine 
months or more between forced visits. 

4. 9 Malfunctions, Downtime, and Repair Costs 

After many months of operation of the ground receive sta- 
tion, it became possible to assess the amount of downtime experienced 
and the nature of malfunctions serious enough to cause the whole 
system to be useless. During most of the early part of this in- 
vestigation, the Landsat was not rigorously tracked every day for 
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reasons other than malfunctions, so there is no valid method of 
quantifying downtime for that period. However, during an all-out 
endurance test in May and June of 1977, near 100 percent dependa- 
bility was attained in ground receive station performance for 60 
days. 

Downtime was experienced occasionally during the investi- 
gation, and it must be viewed in the context of the mission of the 
Water Control Branch: when operationally useful data was needed, 

effort had to be channelled into our operational system, the AHRRN. 
At those times less attention was paid to rapidly correcting mal- 
functions in the satellite ground receive station which was con- 
sidered experimental. Some of the main causes of malfunctions are 
described below. 

1. The radio receiver crystal oscillator failed repeatedly. 

Solution: entire receiver was returned to factory or field service 

organization for repair*; several days at a time were lost. 

2. An RF cable broke or vibrated loose on several occasions 

in the tracking pedestal. Solution: field service technicians 

made on-site repairs. 

3. Computer malfunctions were caused by power interrup- 

tions or surges. This type sometimes led to disc file loss, a time - 
consuming mishap. If the computer aborted early on a weekend, this 
led to loss of whole days of data. Solution: run the computer 

with the power failure protection option. However, if this is done 
on our system, only a small disc is on-line. For matters of con- 
venience, it was usually decided to run without the power failure 
protection. 

4. There were also computer malfunctions associated with 

a long-standing and obscure electrical problem. Occasionally this 
problem was acute enough to disrupt operations for days. It caused 
erroneous reads and writes of the disc files, and most often caused 
false azimuth and elevation commands to be sent to the tracking 
pedestal, thereby causing it to slew wildly and blow fuses. Solu- 
tion: a field service technician found an intermittent open cir- 

cuit in the main power line to the computer. 

5. Our minicomputer is apparently sensitive to high tem- 
peratures. Solution: keep the computer room temperature at about 

18° C. 


^Between December 1976 and July 1977 malfunctions in the antenna and 
receiver subsystems occasioned repair costs of $2,445. 
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It should also be noted that most of the components of the 
system were completely free from malfunctions, including the follow- 
ing: 


- antenna reflector 

- feed 

- waveguide 

- air drying and pumping system 

- par amp 

- all standing cables 

- data demodulator 

- digital input/output interface 

- time code generator 

- S-band test signal generator and feed 

- software (after initial debugging) 


The weakest points in the system were the S-band receiver 
and power instability nriii r^pr^ er. Th e latter problem can be 
solved by power protection devices, and the former can be solved by 
component redundancy. 


4.10 


Cost Comparisons of Data Collection Systems 


Experience with data collection technology during and 
prior to the investigation enables us to estimate the costs of 
various data collection systems. Estimates given below were derived 
from conversations and correspondence with manufacturers and other 
Government agencies. For a fuller economic analysis forecasting 
market conditions and costs, a report by Ecosystems International, 
Incorporated (NASA, 1977), is recommended. Our observations are 
consistent with that report. Estimates have been made for only 
those media that are both available to us and reliable during storms 
that reach hurricane intensity. Those restrictions exclude meteor 
burst and telephone communication, leaving line-of -sight radio and 
satellites . 


4.10.1 Line of Sight Radio (LOS) 


NED ! s Automatic Hydrologic Radio Reporting Network has 
been in use for over seven years. It consists of an automated data 
collection system superimposed over an older voice radio network. 
Experience with the procurement and operation of the AHREN suggests 
that a system can be installed for about $20,000 per reporting 
station. A 40-station network would therefore cost $800,000. 
Maintenance is currently $65,000 per year on just data collection 
hardware, exclusive of voice radio equipment; and the cost is in- 
creasing about 10 percent per year. Maintenance over 10 years 
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(the standard life time adopted for Government -owned equipment) 
would amount to over $1,035,000. Total cost for installation of 
a system of 40 stations and maintenance for 10 years would be over 
$1,795,000. 

4.10.2 Landsat Downlink 


A ground receive station capable of tracking a polar 
orbiting, Landsat-type satellite requires less hardware than the 
radio network. A 12-foot parabolic dish would be sufficient, and 
either an autotracking or step-tracking pedestal could be employed. 
The former type pedestal would simplify the configuration and cut 
costs. Costs may be summarized as follows: 


12-foot antenna with tracking pedestal $100,000 
computer and receiving equipment 100,000 

40 DCP's with interfaces and procurement 

overhead 160,000 

annual maintenance 

computer 8,000 

pedestal 2,000 

40 DCP’s 4,000 

Sub-total 14,000 

10-year maintenance total, 

assuming 10% increase per year 210,000 


Total, including procurement and 

10-year maintenance $570,000 


The computer mentioned above consists of minicomputer with 
32 to 64 K memory, simple disc storage, 1/0, interfaces, etc.. 
About $40,000 would be spent for these items. Receiving equipment 
comprising preamplifier, receiver, decoder, cabling, etc., costs 
an additional $60,000. 

DCP’s currently can be obtained for about $3,000 each, 
but $4,000 is a realistic figure to allow for interfacing, pro- 
curement overhead, batteries, and cabling. A price reduction 
approaching 50 percent can be expected during the next five to 10 
years, .parallelling the price reductions in the semiconductor in- 
dustry. The use of ordinary microprocessors will generalize DCP 
design and cut expense of research and development on tailored 
DCP’s. 


ORIGINAL PAGBJS 
OF POOR QUALITY 
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4.10.3 GOES Downlink 

The equipment for a GOES downlink differs from that of 
Landsat in that a larger dish is needed (18 feet is desirable) , but 
simpler aiming devices can be employed since GOES hovers around one 
point in the sky. If more than one channel is to be received, 
■multiplexing and buffering equipment is needed to route simultan- 
eous messages into one computer. Estimates for components and 10 
years of maintenance may be summarized as follows: 


18-foot antenna with pedestal $ 50,000 
computer and receiving equipment 100,000 
40 DCP f s 160,000 
maintenance 210,000 


Total, including procurement and 

10-year maintenance $520,000 


Costs of DCP’s and maintenance are assumed to be the same 
as those encountered with a Landsat station. 

4.10.4 Significance of Cost Comparisons 

In summary, GOES and Landsat downlinks are approximately 
equal in cost and both are cheaper than LOS radio for comparable 
networks. The cost advantage at this time depends on the existence 
of Government funded satellites (NASA, 1977) . As long as such 
satellites exist, an agency gains a strong advantage by using 
them instead of procuring a whole system. 

5.0 RESULTS AND DISCUSSION OF SATELLITE DATA COLLECTION 

The following items are considered significant results of 
this investigation in the use of the Landsat DCS: 

1. The reliability of data collection platforms, Landsat 
itself, and NED's downlink were judged to be high. This finding 
has led to our continued involvement with satellites as a data 
collection medium. 

2. Component redundancy is an important consideration in 
any operational flood control data collection system. Our exper- 
ience has shown that with many components in series (as they are in 
the DCS) a failure in one component causes an entire breakdown. 

3. A significant outcome of installing the DCS downlink 
was a potential increase in ability to function in the flood control 
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mission with no additional manpower. This potential became apparent 
during the storms which struck New England during this investigation. 
For example, during April and May of 1975 our DCP f s on the Saint 
John River in northern Maine provided hydrologic information for 
relay to flood forecasters at the New Brunswick Power Commission 
office in Fredericton, New Brunswick, Canada. Prior to Land sat, 
part of this area could not be monitored in real time because of 
the high costs involved. In August 1976, the remainder of hurri- 
cane "Belle" travelled through Vermont, New Hampshire, northern 
Maine, and Canada's Maritime Provinces dumping 3 inches of rain in 
many areas. In Canada and Maine local storms dropped up to 2 inches 
during the following week, causing the Saint John River to reach 
near flood stages at Fort Kent, Maine. During this event DCS data 
from the region was received from three of our DCP's and from one 
belonging to Canada 1 s Department of the Environment. The resulting 
high runoff after these storms was studied in connection with the 
proposed Dickey-Lincoln School dams to be built in that area, and 
it was found that for this watershed creditable flood hydrographs 
could be generated from DCS data. Corps personnel who have seen the 
system in operation have deemed the potential sufficient to warrant 
further investment in the field of satellite data collection. 

4. Service was required in the field for malfunctioning 
DCP f s or battery replacement on the average of every nine months, 
and the method that evolved was to visit them as needed, rather than 
to make rounds on schedule. A set of strong arguments favoring 

our adopted method is found in recent research (NASA, 1977). 

5. A single Landsat is suitable for data that is needed 
daily or twice daily. Slowly changing or integrated parameters 
such as wind passage and snow accumulation may be monitored with 
ease; but the system is not suitable for rapidly changing or vari- 
able data such as peak rainfall or insolation. In general it is 
not adequate for collecting hydrologic data in small or flashy 
watersheds. 

6. Costs of satellite data collection systems are less 
than comparable line-of-sight radio systems as long as the invest- 
ment in the satellite is borne by a single agency. It is recognized 
that commercial satellite may be economically favorable within a 
few years (NASA, 1977). 

7. Experience with Landsat and exposure to GOES lead to 
the conclusion that for reservoir regulation an ideal data collec- 
tion system would draw features from both satellites. Our needs 
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would be best suited by a geosynchronous satellite having random 
reporting, the timing of which could be governed by the DCP. 
Alternatively, a geosynchronous satellite operated to allow hourly 
reports or emergency reports would be acceptable. 

8. In August 1977 we invited an electronics team from 

the National Space Technology Laboratory, Bay St. Louis, Mississippi, 
to perform tests on the 15-foot antenna at NED, to determine its 
adequacy to receive GOES data. As a result of these tests and 
our own favorable evaluation of GOES for use in reservoir control, 
we are adding GOES acquisition capability to the ground receive 
station. lie intend to retain Landsat tracking capability as a 
backup in case the entire GOES system fails. Appropriate elec- 
tronic switch-over devices are being built into the combined sta- 
tion. Receiving and computer equipment is on order, and installa- 
tion is scheduled for May 1978. 

9. NED is strongly considering further use of general 
purpose microprocessors for DCP control, data handling at the re- 
mote site, threshold detection, encoding and modulation, along the 
lines of the work already performed at the University of Tennessee 
(NASA, 1976b). The capabilities of microprocessors as control and 
monitoring devices have been adequately demonstrated recently, and 
cost savings and greater flexibility are anticipated with their 
use. 

6.0 BACKGROUND ON IMAGERY STUDY 

6.1 Project Sequence 

The U.S. Army Engineer Division, New England and the Cold 
Regions Research and Engineering Laboratory have been involved in 
the Landsat DCS and imagery analysis since the launch of ERTS-1, 
now known as Landsat-1, in July 1972 (Cooper, et al , 1975). During 
the Landsat-1 experiment CRREL participated in the DCS studies by 
developing sensor interfaces for the Landsat data collection plat- 
forms (DCP f s) and evaluating performance of DCP installations. 

During the last two years (1975-1977) of the Laudsat-2 program 
CRREL was involved in the digital processing of the Landsat com- 
puter compatible tapes (CCT T s) and in sensor interface development 
for the DCP * s . 

6. 2 Hydrologic Parameters 

A coordinating committee comprising personnel from NED, 

CRREL and the University of Connecticut selected the following 
hydrologic parameters that would affect reservoir operation and 
management: 
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Snow cover (areal extent, water equivalent correlation) 

Soil moisture regimes 
Wetlands delineation 
Slope/ topography 
Ice cover 
Flooded areas 

The hydrologic parameters selected for detailed Landsat 
imagery analysis using the Landsat digital data were snow cover and 
delineation of wetlands and flooded areas, 

7.0 DIGITAL PROCESSING OF THE LANDSAT CCT'S 

7 . 1 Description of the Landsat CCT’s 

Landsat-1 and Landsat-2 circle the Earth in a 572-mi 
(920 km) near-polar orbit once every 103 minutes, each completing 
approximately 14 orbits per day. The multispectral scanner (MSS) 
on each satellite continuously scans perpendicularly to the space- 
craft’s direction with an oscillating mirror (NASA, 1976a). Six 
lines are scanned simultaneously in each of four spectral bands 
for each mirror sweep, and radiation is sensed simultaneously by 
an array of six detectors in each of four spectral bands from 0.5 
to 1.1 ym (NASA 1976a). The spectral information is transmitted 
in digital form to ground stations. During image data processing 
at the NASA Goddard Space Flight Center a black and white photo- 
graph can be produced depicting an area approximately 115 mi (185 
km) on a side for the following spectral regions: MSS band 4 (0.5- 

0.6 ym) , MSS band 5 (0.6-0. 7 ym) , MSS band 6 (0.7-0. 8 ym) and MSS 
band 7 (0.8-1. 1 ym) . This information is also available in digital 
form on a computer compatible tape, or CCT. 

The standard Landsat CCT was computer processed to pro- 
duce a geometrically corrected tape with observations transformed 
to a UTM (Universal Transverse Mercator) projection. This geo- 
metrically corrected CCT comprises 2,432 scan lines with each scan 
line covering 3,200 pixels* (Ungar, 1977). Differing levels of 
radiant energy for each pixel within the scene are registered on a 
scale from 0 to 127 (minimum to maximum) for bands 4, 5 and 6, and 
0 to 63 for band 7 (Thomas, 1975). 

7.2 Description of the Computer Analysis Algorithm 

The geometric correction of the digital data and the com- 
puter classification algorithms used in the analyses were developed 
at the NASA Goddard Institute for Space Studies (GISS) (Ungar, 1977). 


*Picture element, representing an area on the ground having dimen- 
sions of 61 x 76 meters. 





The geometric correction provides for a 1:24,000 scale computer 
printout which permits accurate location of test sites. The com- 
puter classification algorithms developed for analysis of the digital 
data allow for both components of the data (each of the four wave- 
length bands and associated energy value for each pixel) to be 
evaluated when classifying the Landsat data into various categories. 
In addition, atmospheric corrections were applied to the Landsat 
digital data (Ungar, 1977). 

The Landsat MSS observation (pixel) may be thought of as 
a point in a four-dimensional "color” space, where the values along 
each axis represent the radiant energy received by the satellite 
in one of the four bands (figure 5a) . Observations lying in a 
similar direction from the origin in this four-dimensional color 
space are said to be similar in color regardless of their total 
radiant energy. The distance of an observation from the origin is 
a measure of the total radiance associated with that data point. 

The algorithm is primarily designed to combine observations that 
are similar in color into the same classification category. There 
are provisions for evaluating brightness differences and for weigh- 
ing these differences in with the color discrimination when con- 
structing the classification categories. 

Discrimination based solely on color is obtained when the 
difference in direction between the color vectors (observations) 
is examined. If the angle between the observations is smaller than 
some user-defined aviterion y the vectors are considered to be lying 
in the same direction and*, therefore, the observations are placed 
in the same category. 

There are two modes in which this classification scheme may 
be employed, supervised and unsupervised. In the supervised mode 
the user specifies a signature (the energy distribution in four 
Landsat bands) . If an observation lies within a solid angle smaller 
than the user-defined criterion, Stnax* it is said to belong to the 
category represented by the multispectral signature (figure 5b) . 
Therefore, all vectors lying within a cone of angle, about the 

signature representing category X belong to category X. 

In the unsupervised mode the color vector corresponding to 
the first observation is compared with all subsequent observations. 

If color vector 1 is similar In direction to color vector 2 (i.e., 

<$9 JS ^max) > observation 2 is placed in the same category as the 
first observation (figure 5c). In a similar fashion observations’ 
subsequent to observation 2 are compared to the second observation 
and so on right up to the last observation. If in the process of 
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a. A color vector in a four-dimensional space. 


Band 3 



Band 2 

b. Supervised mode . The user-defined criterion, 6 max , defines 
category X about the specified signature If. Any color vector that 
lies within this cone belongs to category X. This is illustrated 
for three bands, however, all four Landsat bands are used in the 
computer classification algorithm. 


Band 3 



Band 2 

c. Unsupervised mode . B 3 is similar in direction to B 3 (60 6 max ) 
and placed in category l.^B^ is similar in direction to §2 and placed 
in category 2. However, B 4 is also similar in direction to §3 (cate- 
gory 1). Therefore, category 1 is merged with category 2. 

Figure 5. The concept of the four-dimensional "color" space used in 
computer classification algorithm. 




constructing categories, a member is found which belongs to a pre- 
vious category, the new category is chawed (or linked) to the 
original classification category forming one joint category (Ungar, 
1977) . In effect the unsupervised classification will form several 
categories based on a criterion specifying maximum color difference 
permissible between members of the same category. 

In addition to discrimination based solely on color, the 
GISS algorithm provides the capability of weighting total radiance 
differences into the discriminant equation for classification. The 
percent difference in brightness between two observations is com- 
puted. The calculated normalized difference is then combined with 
the color difference angle (expressed in steradians) by performing 
a weighted average in the RMS (root mean square) sense. This 
brightness-weighted quantity is now compared with the user-defined 
criterion (6 max ) . Thus, in the classification process, a relatively 
small weighting of brightness allows very large brightness differ- 
ences to disqualify observations that are similar in color from 
membership in the same category, thereby adding a second level of 
discrimination. Discrimination of the classification categories 
based partially on overall brightness differences plays an impor- 
tant role in the work discussed in this report. 


7 . 3 Computer Data Handling and Analysis 


A Harris COPE 1200 remote job entry terminal was utilized 
at CRREL for computer data handling and analysis. The remote ter- 
minal was used in the analysis of the DCP data cards (hexadecimal 
format) obtained from NASA and the digital processing of the Landsat 
CCT's, 


Data reduction of the DCP cards was accomplished by using 
the Infonet computer system located in Chicago, Illinois. The 
computer programs were developed at CRREL and used for analysis of 
the data from each interfaced DCP sensor. These data included 
temperature, using thermocouples, and water equivalent from a snow 
pillow. 

The digital processing of the Landsat CCThs was accomplished 
through a cooperative agreement with NASA GISS, Computer algorithms 
for the analysis of the digital data were developed at GISS (Ungar, 
1977). These algorithms were accessed using the CRREL remote entry 
terminal to the main computer facility located at GISS in New York 
City. 



8.0 


SNOW COVER ANALYSIS 


8.1 Literature Review 


Manual methods have been used to delineate the areal 
extent of snow and the mean altitude of the snowline from Landsat 
photographic data products (Barnes and Bowley, 1974; Meier, 1975a, 
1975b). However, .a quantitative measure of the water equivalent 
of the snowpack has not been obtained from Landsat photographic 
data products. Usually the areal extent of snow has been related 
indirectly to subsequent watershed runoff occurring during the 
springtime (Meier, 1975c; Anderson et al f 1974). Also, the changes 
in the areal extent of snow cover measured on Landsat imagery have 
been found to correlate with changes in water equivalent recorded 
by a snow pillow (Anderson et al, 1974). Another Landsat manual 
interpretation method has used a coded snow cover classification 
scheme to account for vegetation cover, density, aspect, elevation 
and slope to map the areal extent of snow (Katibah, 1975). 


A snow mapping experiment comparing the identification of 
six snow cover types was accomplished using three image processing 
systems - LARSYS Version 3, STANSORT-2 and General Electric Image- 
100 (Itten, 1975) . In addition, other studies have focused on 
digital analyses of Landsat data in defining various snow cover 
types (Bartolucci et al, 1975; Dallam, 1975, Luther et al, 1975; 
Alfoldi, 1976). In these studies a quantitative estimate of water 
equivalent content associated with snow cover types was not made. 
In one case it was suggested that spectral variations within the 
snowpack area could not be reliably determined because of detector 
saturation problems (Bartolucci et al 1975). 


Another study used simulated infrared Landsat color composites 
and snow course data to estimate water equivalent related to the 
snowpack (Sharp, 1975). Sampling units on the Landsat image were 
mapped to determine the areal extent of snow. An estimation of a 
snow water content index was calculated using a linear regression 
equation relating the imagery to ground truth data. 


It has been stated that remote sensing of the snow cover 
may have useful applications since the magnitude and wavelength of 
reflectance vary with snow types (Mellor, 1965). Also, the albedo 
is high for a layer of new snow and as the new snow coalesces and 
coarsens in texture the albedo falls steadily (Bergen, 1975) . In 
addition, a reduction in the spectral reflectance occurs from the 
combination of densif ication and increased particle size associated 
with aging (O’Brien and Munis, 1975). Therefore, it is believed 
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that the Landsat CCT f s may contain information that can be used to 
estimate water equivalent in a snowpack. If so, then Landsat 
digital data can be used for estimating volume of spring runoff from 
a watershed more accurately than is now possible. 

8.2 Approach 

8.2.1 Selection of Site 

The Dickey-Lincoln School Lakes Project, Maine (figure 6), 
currently being evaluated by NED for the generation of hydroelectric 
power, flood control and recreational purposes, was an idea], site 
for an analysis of snow cover utilizing the Landsat CCT's. Due to 
the remoteness and inaccessibility of the area there is not an 
adequate data collection system for evaluating the water equivalent 
of the snowpack each year. 

The climate of this region is characterized by short, cool 
summers and long, cold, windy winters. Average annual temperature 
is 39°F (3.9°C) with extremes of -40°F (-40°C) in the winter and 
99°F (37°C) in the summer (New England Division, Corps of Engineers, 
1967). The average annual precipitation is approximately 36 inches 
(91 cm) and occurs uniformly throughout the year with about 30 per- 
cent in the form of snow. The average annual snowfall is about 
100 inches (254 cm) which occurs during 8 months of the year (NED, 
1967). Average winter snow depth ranges between 20 and 40 inches 
(51 and 102 cm) with the upper limit exceeding 50 inches (127 cm). 
Water equivalent of the snowpack reaches a maximum in late March 
and usually exceeds 10 inches (25 cm) . The geology and vegetation 
of the Dickey-Lincoln School Lakes Project area was previously 
mapped and served as a data base of site characteristics in this 
study (McKim and Merry, 1975; McKim, 1975). 

Ground truth for the snow cover analysis consisted of 
selected snow courses obtained in the upper Saint John River basin 
by the U.S. Geological Survey and the Allagash Wilderness Waterway 
Agency for the winter season of 1972-1973 (U.S. Department of 
Commerce, 1973) (figure 7) . The Allagash B snow course is located 
within a mixed forest near the confluence of the Allagash and Saint 
John Rivers at an elevation of about 640 feet (195 m) msl and is 
characterized by a southeast exposure with gently sloping terrain 
(figure 8) . The Beech Ridge snow course is located within a mixed 
forest near the Front ier-Chur chill Road near Umsaskis Lake at an 
elevation of about 1,300 feet (396 m) msl and is characterized by > 
a western exposure on gently sloping terrain (figure 8) . The 
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Ninemile B snow course is located within a coniferous forest on 
the floodplain near the USGS gaging station on the Saint John River 
at an elevation of about 950 feet (290 m) msl and is characterized 
by a northwest exposure (figure 8) . 

8.2.2 Snow Course Data 

The available cloud-free Landsat CCT's of the upper Saint 
John River basin were selected for four seasons and included the 
following dates: 11 February 1973 (image ID 1203-15002), 23 July 

1973 (image ID 1365-14593), 26 November 1973 (image ID 1491-14572) 
and 19 April 1974 (image ID 1635-14541) . It was not possible to 
obtain a cloud-free Landsat scene for February 1974, which would 
have been desirable. Therefore, the 11 February 1973 CCT was 
selected for the snow cover analysis. 

A 

The snow depth and water equivalent were estimated from 
figure 7 for the three snow courses (table 3) for 11 February 1973 
(U.S. Department of Commerce, 1973; Li and Davar, 1975) • These 
data seemed reasonable when compared to local climatological data 
obtained at Caribou, Maine, located 50 airline miles (80,5 km) east- 
southeast of the Dickey-Lincoln area. A major snowstorm occurred 
on 29 January 1973 between the last snow course measurement (28 
January 1973) and the date of the CCT (11 February 1973). A total 
of 9.2 inches (23.4 cm) of snow associated with a water equivalent 
of 0.65 inches (1.6 cm) was recorded at Caribou, Maine. In addi- 
tion, on 8 February 1973 there was a minor snowfall of 1.6 inches 
(4.1 cm) with a water equivalent of 0.12 inch (0.3 cm). Based on 
these data, the estimated 9,5 inches (24.1 cm) was assumed to be 
reasonable. 


TABLE 3 

SNOW COURSE DATA 
(11 February 1973) 


Snow 



Latitude/ 

Course 

Sampling 

Snow 

Water 

Snow Course 

Longitude 

Length 

Points 

Depth 

Equivalent 



(ft) 


(in)* 

(in)* 

All ag ash B 

47°05'N/69°04 , W 

1,000 

10 

42 

9.6 

Beech Ridge 

46°36’N/69°28 , W 

- 

- 

41 

9.4 

Ninemile B 

46°42'N/69°43’W 

1,000 

10 

38 

9.5 


*1 inch = 2.54 cm 



TT 








j ■ t / r: i i 




’ X r«*nr. : r as , 


3BSSSL3: f ^ - W JT ■ tITTv*:] 


The snow courses were located on each of the CCT’s by 
generating a geometrically corrected, 16-level grayscale computer 
printout of a 320 by 256 pixel area (146.9 mi^ or 380.6 km^) at a 
scale of 1:24000. Each observation was assigned one of 16 levels 
of gray depending on its radiance value in MSS band 7. The snow 
course sites were located on the grayscale printouts using avail- 
able topographic njaps for orientation. 

The computer test site containing each snow course is 40 
by 32 pixels for a total area of 2.3 miles^ (6.0 km^) . The snow 
course was located in the center of each computer test site. The 
computer algorithm described previously was applied to extract 
information concerning the spectral characteristics of the snow 
cover /vegetation within the snow course computer test sites. 

8 . 3 Results and Discussion 

Unsupervised classifications were performed on the three 
snow course sites for the 11 February 1973 CCT for a range of <$niax 
(delmax) values between 0.02 and 0.04 with several brightness 
weightings (for example, 0.1, 0.2, 0.3) for initial classification 
of the digital data. Computer runs which produced large numbers of 
categories were selected so that several signatures could be ex- 
tracted for the pixels contained within the snow course areas. 

This allowed for an evaluation of signature variation within each 
snow course. The sun elevation angle (23°) of the scene was cor- 
rected to zenith to account for seasonal variations in irradiance. 

The three snow course computer classification printouts for 
11 February 1973 are shown in figures 9, 10, and 11 (Merry, et al, 
1977) . The location of the snow course is shown outlined on the 
computer test site. The radiance values associated with each pixel 
within the snow course are indicated by the arrows within the cate- 
gorization summary shown in figures 9, 10 and 11. In this summary 
the two left-hand columns are the category symbol and the number of 
pixels {man) within each classification category, respectively. 

The tot column is a measure of the variation between the signatures 
within a category, with the smaller numbers indicating little vari- 
ation. The normalized radiances in each band (B1 f 3 B2 f > BS f and 
B4 f ) always add up to unity and are listed for each category. Also, 
the true radiances are listed for each band ( Bl 9 B2> B3 and B4) ; 
these values sum to the total radiance (mW/cm^ sr*) listed in the 
extreme right-hand column. The delmax and the brightness weighting 
used in the unsupervised classification are shown at the bottom of 
the categorization summary. 


^Milliwatts per square centimeter per steradian. 




ALLAGASH B SNOW COURSE 

640' ELEVATION ASPECT: SOUTHEAST 


WATER EQUIV. 9.6" 
SNOW DEPTH 42" 
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Figure 9 Computer classification printout of the Allagash B snow course (ll February 1973) 





:iv ■ 


'".V; 


yiaii&ay&aay^ 




„ ", -r gfeggggggfta 


C*f£fiU«II»TIUh »U«*4B? 




PUJTJ2C0 OUTPUT UF SUM€e"tM<5) eOu. 
inORtNS «$ CATeSQNlfcl 

/MU//*4/I/*///*///*2Dy*y4/yy !♦//// /*/•,♦» 
/L«////La///MU/»/Uf3 ///*/♦♦♦/♦//// UWn 


out 

s< 

TM 




C*#YY /■■t.//-/YY?T« I 

rj.fi-TT «Mti 


//♦//»♦///♦♦/♦••/! 0 
♦/♦y/0 JU///44//4 // 

T /■■&•//. /y»YT» l.,V IV 

« T}FJ,»«TT *k»»» 

us.th i6syr, % uqobdI oli-ae.u,J!*iCS 
/nIffc)**ae7//**. .*//••*. 

<T) 14T«//U/U£yy/ h/ wjTwuIpS H 

»)«c»rntn Jrt7% t yY*Y2)uuon/, 

»0*tIU)mjYt £*uU®UO®HSTT* ♦ ««. ♦ **J 

A* £O0Au7BQ E,9U*®Iij|l*£// ♦ IW H// 0 

/ //*!.; »E/,«/»7tff!l*4 5 B 

/? «j? 

• %Q»»ACveuov« /*&ce f n 

I* ..47)77 T«f*4«I»«BJTw. • X. 

0*tu*:omt)TT < uTRujRrJui *♦ 

*!• *.f)y F*afi* /ot )&*.♦ r •/ (i 

I t«T» ♦////U*0MWT/NQU//*/ *♦ 0/ 


•J |«**Y Y*QOOy«*B//// ♦/ 


D/Y* 

1M//JO 


♦ 0 fc 

10 * 1 / 


/ 140 

^ * l*° 


► • 
► * 
► * 

» * 

) 

t 

9 

T 

Y 

► U 

► t 

► 0 

I 

9 

A 

9 


it? 

Tl 

9* 

10 

• 9 

• J 
«U 
3* 
3? 
3* 
St 
13 
12 
to 
tu 

1 


BEECH RIDGE SNOW COURSE ! 
I300‘ ELEVATION ASPECT: WEST ► . 

WATER EQUIV. 9.4" SNOW DEPTH 4l“ 



nON*** 

TBl*C 

Lllfel* **OI*5«Cfc3 / 

«*ci ANCt» 


r i1 L. 

*1 • 

04 * 

*3 * 

*4' 

TOTAL 

91 

*0 

*3 

94 

RADIANCE 

t.i 

0.13 

0,1« 

0 • ON 

0,4* 


0*74 

0.53 

0.3? 

«.ll 

4,0b 

U.9 

0,1 2 

0.1« 

U,?8 

0,04 

0.«F 

9.74 


0,70 

0,1* 

J.*J 

v.9 

8:11 

0.1** 

c.eo 

8$ 

U.4? 

3 »*2 

9.84 

1 .2 

0,11 

0,13 

0.04 

0,4* 

4.94 


u*»o 

Q.F9 

0.34 

4.44 

0,9 

o.ts 

0.14 

0,04 

U.44 

9,3* 

0,70 

0.73 

0.3* 

3,49 

0.9 

0,11 

0,l« 

0.04 

0.4? 


ojio 

U, 79 

0.39 

4.7* 

9,9* 

U.% 

0,11 

0,13 

0.04 

0 ,0* 



0.65 

0.4« 

0.30 

3.99 

9.14 

0,B 

0,1 1 

0.13 

0,09 

0.70 


0 .49 

0.77 

0,3t 

i.W 

9.40 

0,T 

t-M 

8:H 

0. 04 

0.J4 

O.J? 

4.7* 

9,40 

o,e' 

8:1! 

8:18 

8:1$ 

S:?! 

9,*1 

0.9 

0,1* 

o:?9 

i:\i 

0.04 

0.39 

5:52 

9.4* 

0,4 

0,14 

8:K 

0.07 

0.44 

9,44 


0,40 

0.3? 

4.7* 

u.? 

0,14 

0.1* 

0.0? 

0.43 

4.02 

0*49 

0,49 

o.»o 

4.47 

9.9 

0,1« 

0,13 

0.04 

«.*? 

9.87 

u,*o 

0.7* 

0.4? 

3.47 

Q .3 

u.ll 

0,13 

0,04 

0.70 

3.83 


0.45 

0.79 

0.32 

«.ll 

0,* 

0,13 

0.1« 

0.04 

U, 4 • 

*,44 

0.49 

0.4? 

0.29 

4.34 

0,4 

0,1* 

v.13 

0.04 

0.4* 

9.11 


o.»o 

0.4* 

0,2* 

4.93 

0.4 

U,13 

0.19 

0.04 

0.43 

4.02; 


0,40 

0.91 

0,3* 

5.9* 

0.7 

U, 19 

U, 14 
1.0* 

0.07 

0,42 

4.9* 


0,44 

G.a? 

4.11 

0.9 

0,17 

0,17 

0,5* 

0.09 

0,71 



0I4O 

0.29 

5.99 

».« 

0,1 

U.13 

0.79 

0,19 

0.07 

0,*5 

5,76 


C ,89 

0,40 

4.7* 


NOBmai izco RADIANCES * 
TRUE RADIANCES 



NUM. 

mu 

Bl* 

81 

82 * 
82 

8 J 1 
93 

84' 

8* 

TOTAL 

radiance 


M 

7 

o.« 

O.U 

0.40 

0.13 

0,80 

o.o? 

3,4 1 

0,44 

4.7* 

9.8* 


4 

0.3 

0.11 

0.47 

0.13 

0.83 

0.04 

0.43 

0.4* 

4,4? 

4.49 


4 

0.9 

0.14 
r • *5 

C, if 
u,*l 

O.U? 
0, *1 

0,43 

4.?* 

9.81 


4 

0.1 

0.13 

0.49 

0,1 « 
0,40 

0,07 

0.48 

0,49 

i:?« 

3.?4 


9 

0,4 

o.n 

0,70 

o.ll 

0.71 

0,09 

n.l* 

U, 72 

4,48 

8.21 


« 

0.9 

o.l* 

3.44 

0,19 

1,00 

0,04 

0.43 

0.49 

4.0* 

8.87 


a 

0.4 

0.1* 

O.fC 

o.l* 

0.7* 

C .04 
0.40 

0,44 

l!lk 

9.10 


• 

0.4 

0.19 

0.89 

0.15 

0,04 

C.O? 

0,* 1 

0,43 

l|3! 

9.64 


4 

0.2 

0.19 

0,80 

0.19 

U.V3 

C.O? 

0.83 

0.43 

4.92 

6.19 


4 


0.12 

0.4* 

•5:1! 

0.05 

0.14 

0.?1 

4,4* 

6.86 


4 

0.2 

o.n 

U.49 

8:15 

0.05 

Q.2T 

0.72 

*!n 

9.72 


4 

0.4 

0.15 

0.80 

0.1* 

0,78 

0*04 

0.39 

0,49 

3.39 

9.*0 


4 

0.4 

0.19 

1.02 

0. 19 

1. C2 

0,07 

0.4* 

0.4* 

4.4* 

7.01 


4 

0.* 

0,1* 

i!o* 

0.14 

0.88 

0 , ON 
0,44 

0,49 

*.*? 

?.l* 


3 

0.4 

0.1* 

0.80 

0.1* 

0.*8 

0,0? 

0.3* 

0.43 

4.99 

9.42 


J 

U.4 

0.13 

0,83 

0.14 

1.01 

0.07 

0.0 

0.42 

3.82 

6.31 


3 

0.3 

0.12 

C.88 

0,14 

1.03 

0.04 

0,*4 

0,48 
4, 8N 

?.!* 


3 

0.2 

0,14 

0.89 

0,12 

o.ei 

0.00 

0.-0 

0,04 

*.*? 

4.73 


4 

0.2 

0.13 

0.40 

o.n 

0,70 

0,05 

0.32 

0.71 

4,48 

6,3b 


4 

0.3 

0.14 

l.IT 

0.19 

1.24 

0,0? 

0.9? 

i:ll 

8.«3 


3 

0,3 

0.17 

0.80 

0,12 

0.8Q 

0.05 

0,3? 

0,72 

5,91 

7.77 


3 

0.3 

0.11 

0,70 

0,11 

0.75 

C.09 

0.3* 

0,?3 

4.84 

6.63 


MO, F0t*T8 CLA88IFU0 ■ 10)0 NO, *CInT* UMCLA99IEJ1D • 

UClMAI • O.OJtO BRIGHTNESS WEIGHING -} 0,iOO 


100.00 »e 
ino.uo pi 
too, oo Pk 

100.00 Rfe 


CENT OF AND 
CENT OF BAND 
CfNT OF band 
cent of BAND 


1 U9E0 

I u S r ‘ 

3 Uli 

4 Utl 


III) In 
I feO IN 
IfO IN 
€0 IN 


ALulOO CALC. 
AlblOO CALC. 
ALBEDO CalC. 
albeoo CALC. 


Figure 10 Computer classification printout of the Beech Ridge snow course (ll February 1973) 
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Figtne 11 Computer classification printout of the Ninemile B snow course (ll February 1973). 





The Allagash B snow course computer classification printout 
is shown in figure 9. The total radiance of pixels contained in the 
snow course area varied from 6.93 to 7.74 mW/cm2 sr (categories: 

E T & P 2 6), which corresponded to a water equivalent of 9.6 inches 
(24.4 cm) obtained from the snow course data. An important obser- 
vation was that the radiance for MSS band 7 was consistently 2.99 
or 3.18 inW/cm^ sr, a difference of only one energy level. 


The Beech Ridge snow course computer classification print- 
out is shown in figure 10. The total radiance of pixels contained 
in this snow course varied from 5.34 to 6.54 mW/cm^ sr (categories: 
* . + Q W U I 0 G) and corresponded to a water equivalent of 9.4 
inches (23.9 cm) obtained from the snow course data. The Ninemile 
B snow course computer classification printout is shown in figure 
11. The total radiance for the pixels contained in the site varied 
from 5.45 to 6.87 mW/cm^ sr (categories: / * W C) and corresponded 
to a water equivalent of 9.5 inches (24.1 cm) obtained from the 
snow course data. 


The radiance varied from 5.34 to 7.74 mW/cm^ sr over the 
three sampled snow course areas, which corresponded to a water 
equivalent value of approximately 9.5 inches (24.1 cm) (Merry, et 
al, 1977). The range in radiance values from 5.34 to 7.74 raW/cm^ 
sr may be attributed to variations in vegetative cover, slope and 
aspect which occurred among these snow course areas# 


The greatest radiance occurred in cleared areas such as 
fields and river channels. As an example, the snow cover on the 
Saint John River (figure 9 and 11 showed the highest radiance, 
which ranged from 20.23 to as high as 29.22 mW/cm 2 sr. These high 
radiance values occurred in areas where there was a minimal vege- 
tative cover. The important factors to be considered in the snow 
cover mapping and assessment of the water equivalent analysis based 
on radiance values obtained from the Landsat CCT’s in the Dickey- 
Lincoln School Lakes area are probably the vegetative cover, slope, 
aspect, geomorpliic position and, to a lesser degree, elevation. 

A preliminary analysis of the 11 February 1973 CCT using 
the GISS computer algorithm showed that the radiance of snow cover/ 
vegetation varied from approximately 20 mW/cm 2 sr in non-vegetated 
areas to less than 4 mW/cm 2 sr for densely covered forested areas. 
Comparison of the digital data from three snow courses in the Dickey- 
Lincoln School Lakes area to the radiance value of the snowpack at 
these sites indicated that the radiance of the pixels contained in 
the snow courses varied from 5.34 to 7.74 mW/cm 2 sr, with the 
average radiance value being 6.4 + 0.6 mW/cm 2 sr. The water equiv- - 
alent of the snowpack for this range of radiance values was approxi- 
mately 9.5 inches (24.1 cm) of water. 


42 






Average mill t i spectral signature*; were extracted for the 
three snow course sites from the four-hand energy values for 
11 February 1973 and 23 July 1973 to study seasonal variations in 
radiance. The multispectral signatures from these dates are shown 
in table 4. 

TABLE 4 

AVERAGE MULTISPECTRAL SIGNATURES (mW/cm 2 sr) 

FOR THE THREE SNOW COURSES 

Total 


Snow Course 

Date 


MSS 4 

MSS 5 

MSS 6 

MSS 7 

Radiance 

Allagash B 

ii 

Feb 

73 

1.741 

1.189 

1.209 

3.084 

7.223 


23 

Jul 

73 

0.518 

0.247 

0.603 

1.760 

3.128 

Beech Ridge 

11 

Feb 

73 

1.053 

0.600 

0.720 

1.927 

4.300 


23 

Jul 

73 

0.507 

0.233 

0.561 

1.705 

3.006 

Ninemile B 

11 

Feb 

73 

1.428 

0.534 

0.975 

4.033 

6.970 


23 

Jul 

73 

0.550 

0.282 

0.591 

1.879 

3.302 


The areas where the three snow courses are located showed 
about the same radiance value (3 mW/cm^ sr) for the month of July. 
This low value would be anticipated due to the absence of snow. 

The minor differences observed in the four MSS bands can be attrib- 
uted to the difference in vegetative cover. The Allagash B and 
Beech F-idge sites have a mixed forest cover and their multispectral 
signatures are similar. The Ninemile B site is in a coniferous 
forest and the four-band multispectral signatures are slightly dif- 
ferent from those of the other two sites. 

Figure 12 shows the computer classification of snow cover/ 
vegetation classes for the 11 February scene of a selected area 
(35.2 mi^ or 90 km^) near the confluence of the Saint John and 
Allagash Rivers with the correlative USGS topographic map. The 
multispectral signatures (table 5) for a supervised classification 
were derived from table 4 for the three snow course sites during 
February and another multispectral signature was derived from the 
four-band energy values for the Allagash and Saint John River chan- 
nels. The four snow cover /vegetation classes are shown in table 6. 
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TABLE 5 


INPUT MULTISPECTRAL SIGNATURES (mW/cm 2 sr) 
FOR THE FOUR SNOW COVER/ VEGETATION CLASSES 


Symbol 


Multispectral Signatures 


on Map 

Delmax 

MSS 4 

MSS 5 

MSS 6 

MSS 7 

Total Radiance 

W2 

i 

0.25 

1.74 

1.19 

1.21 

3.08 

7.22 

0.3 

2 

0.25 

1.05 

0.60 

0.72 

1.93 

4.30 

0.3 

3 

0.50 

1.43 

0.53 

0.98 

4.03 

6.97 

0.3 

4 

0.25 

5.22 

4.23 

3.43 

7.52 

20.40 

0.3 


Symbol 
on Map 


TABLE 6 

SNOW COVER/VEGETATION CLASSES MAPPED 
FROM THE 11 FEBRUARY 1973 CCT 


Snow Cover /Vegetation Characteristics 


Mixed forest, 640 ft (195 m) msl elevation, 
southeast exposure, gently sloping. 

Mixed forest, 1,300 ft (396 m) msl eleva- 
tion, western exposure, gently sloping. 

Coniferous forest, 950 ft (290 m) msl ele- 
vation, northwest exposure, level. 

Open nonvegetated areas, lowest elevations, 
level. 


Water 

Equivalent 

(in)* 


*1 inch ~ 2.54 cm. 


Patterns of snow radiance values can be observed on the 
computer classification printout (figure 12) which suggests the 
interrelationship of vegetation, slope and aspect. The symbols 1 
and 2 are predominant, indicating mixed forest at various eleva- 
tions. The symbol 3 occurs in isolated areas such as hilltops and 
along the river channels . The symbol 4 occurs as expected along 
the Allagash, Saint John and part of the Little Black Rivers. 

Also, there were a number of unclassified pixels (the dashes) which 
can be seen in figure 12. This was as anticipated due to other 
snow cover/vegetation classes that were not defined during this 
exercise. It is suggested that these undefined snow cover /vegeta- 
tion classes would be for areas that have water equivalent values 
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either greater or less than 9*5 inches (24 ,1 cm) or other types of 
site characteristics, 

9.0 INLANDS MAPPING 


9. 1 Literature Review 


Landsat MSS imagery has also been used to delineate the 
extent of wetlands. The wetlands maps produced have generally re- 
sulted from a visual interpretation of the Landsat photographic 
data products. Mapping accuracies between 70 and 85 percent have 
been achieved using Landsat MSS photography (Seevers, et al, 1975; 
Anderson, et al, 1973a, 1973b; Higer, et al, 1975; Rehder and 
Quattrochi 1976). 


Accurate inventories of wetlands larger than ten acres 
were made in Nebraska for four categories: open water, sub irri- 

gated meadows, marshes, and seasonally flooded basins. The 
inventories were made by using imagery from two seasons and an 
electronic image-enhancing system. Positive print enlargements 
of MSS bands 5 and 7 at a scale of 1:250000 (acquired in the spring) 
as well as band 7 (acquired in late summer) were used to delineate 
all wetlands (Seevers, et al, 1975). Electronic enhancement of 
MSS band 6 (acquired in the fall) was used as an aid to differen- 
tiate marshes. 


A wetlands map of Wisconsin was prepared at a scale of 
1:500000 using Landsat MSS bands 5 and 7 and an additive color 
viewer as a data analysis system (Frazier, et al, 1975). Wetland 
areas in this investigation were defined as those which had 
enough water during June to adversely affect the infrared reflec- 
tance of growing plants. These included areas with wetland cover 
types (marsh, sedge meadow, shrub-carr and lowland forest) and 
poorly drained agricultural cropland areas. The primary criteria 
for delineation of wetlands were the reduced infrared reflectance 
of broad leaf plants growing in wet areas, the dark red tone of 
spruce bogs, and the black color of organic soil areas observed 
when using the additive color viewer (Frazier, et al, 1975). 

Significant changes occurred in the size of wetlands be^ 
cause of seasonal fluctuations in vegetation characteristics and 
precipitation (Rehder and Quattrochi, 1976) . The dynamic charac- 
teristics of the wetlands were not attributed exclusively to 
seasonal factors, as significant changes in wetland morphology 
were found to occur within the MSS bands for the same date. 


The following features were determined from Landsat MSS 
band 5 and 7 imagery enlarged to a scale of 1:250000 for test 
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sites located in Maryland and Georgia: (a) upper wetlands boundary, 

(b) drainage patterns within the wetlands, (c) plant communities 
such as Spartina altemiflora , , Spartina patens 3 Juneus roemerianus , 
(d) drainage ditches associated with agriculture, and (e) lagoons 
for waterside home development (Anderson, et al, 1973a, 1973b). 
Mapping at a scale of 1:250000 was adequate for the general deline- 
ation of large marshes and for rather gross wetland species asso- 
ciations. 


In addition, digital processing of Landsat MSS imagery 
has been used* to map wetlands (Anderson, et al, 1973b; Cartmill, 

1973; Flores, et al, 1973; Klemas, et al, 1973). Mapping accuracies 
ranging generally from 78 to 99 percent have been achieved using 
digital processing techniques. 

Seven categories of marsh vegetation and marsh features 
were identified at an approximate scale of 1:20000 (Anderson, et al, 
1973b). These categories included: water, sandy mud flat, organic 

mud flat, sparsely vegetated, spoil, Iva frutescens , Spartina patens 3 
and Spartina alterni flora. 

Seventeen to thirty spectrally homogeneous land use classes 
were defined in the Texas coastal zone using two clustering algo- 
rithms available from the NASA Johnson Space Center (Flores, et al, 
1973). Many classes were identified as being homogeneous features 
such as water masses, salt marsh, beaches, pine, hardwoods, and 
exposed soil or construction materials, with most classes identi- 
fied as mixtures. 

Eight vegetation and land use discrimination classes were 
selected to map and inventory the significant ecological communi- 
ties in the coastal zone of Delaware (Klemas, et al, 1973) . These 
classes were: Phragmites communis (giant reed grass), Spartina 

altemiflora (salt marsh cord grass), Spartina patens (salt marsh 
hay), shallow water and exposed mud, deep water (>2 m) , forest, 
agriculture, and exposed sand and concrete. The Spartina alterni - 
flora was discriminated with an accuracy of 94 to 100 percent, the 
Phragmites communis showed a classification accuracy of 83 percent, 
but the discrimination of Spartina patens was only 52 percent. 

9.2 Approach 

The site selected for the wetlands mapping analysis was a 
124-km^ (48-mi^) area of the Merrimack River estuary (figure 6). 

This area contained the largest variety of land use and vegetation 
classification units to be found in the Merrimack River basin. In 
addition, the Merrimack River basin was a primary test site for the 
NED-CRREL Skylab Earth Resources Experiment Package (EREP) project 
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(McKim, et al, 1975b). Ground truth data in the form of land use 
maps prepared from satellite and aircraft photographic data products 
were available for this site. 

9.3 Results and Discussion 


The CCT's were obtained of the Merrimack River estuary for 
6 July 1976 (image ID 5444-14082). A grayscale computer printout 
of the Merrimack River estuary was obtained for the purpose of lo- 
cating potential training sites for wetlands. In addition, a 
1:24000 NASA RB-57/RC-8 photograph and a USGS topographic map of 
the study area were available for ground truth comparison to the 
computer printout. An overlay was prepared from the photograph 
showing the delineation of the water and the extent of wetlands. 

Two training sites were located on the north side of the estuary 
for use in the wetlands analysis (figure 13) . 

A n ground truth" computer algorithm developed at GISS was 
used to determine an average multispectral signature for the wet- 
lands and the water classification categories. This computer 
program allows one to "tag" certain pixels of a category into a 
32 by 40 array; these specified pixels in the array are then used 
in the computation of an average multispectral signature. 

The wetlands overlay that was prepared from a photointer- 
pretation analysis of the 1:24000 photograph was placed over the 
grayscale computer printout (scale 1:24000) of the Merrimack River 
estuary. The pixels that were within the boundaries of the de- 
lineated wetlands and the water classification categories were 
tagged as the pixels to be used in the specified array for the 
ground truth computer program. The ground truth computer program 
determined the average multispectral signatures for the two spec- 
ified categories (table 7). 

TABLE 7 

AVERAGE MULTISPECTRAL SIGNATURES (mW/cm 2 sr) 

FOR THE WETLANDS AND WATER CATEGORIES 


Category 

MSS 4 

MSS 5 

MSS 6 

MSS 7 

Wetlands 

0.584 

0.326 

0.581 

1.510 

Water 

0.551 

0.237 

0.140 

0.131 


The multispectral training signatures (table 7) were used 
in a supervised classification of the northern portion of the 
Merrimack River estuary. Various values of delmax were used until 
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an optimum computer classification map of the wetland areas was 
obtained. The computer classification map is shown in figure 14. 
The symbol 1 indicates a wetlands unit, the symbol 2 indicates a 
water unit, and the dashed lines (-) indicate unclassified pixels. 
The wetlands/water photointerpretation map overlay was used in 
comparing the accuracy of the computer map of wetlands. A classi- 
fication accuracy of 75 percent was obtained for the wetlands unit, 
taking into account the misclassif ied and the unclassified wet- 
lands pixels. 

The reason for the 74 percent accuracy may be the vari- 
ability of wetland species, because it was assumed that all the 
wetlands contained the same vegetation type. If there was a large 
variability in species, there would be different multispectral 
responses. Also, changes in moisture conditions and tidal fluc- 
tuations would contribute to multispectral variations. This would 
prevent the wetland areas from being classified in one broad cate- 
gory . 


10.0 MAPPING OF FLOODED AREAS 


10. 1 Literature Review 


Landsat MSS data have been used for flood observations 
because of the relatively high resolution, cartographic fidelity 
and the near infrared sensors (Rango, 1975). Flood area measure- 
ments for areas of 100 km2 (39 mi^) or more have been made with 
less than 5 percent error (Rango and Salomonson, 1974) . 


Flood-prone areas have been shown to have multispectral 
signatures indicating categories of natural vegetation, soil 
characteristics and cultural features which are different from the 
signatures of surrounding non-flood-prone areas (Range and Anderson, 
1974) . These differences have developed over a period of time in 
response to increased flooding frequency which enabled the signatures 
to be distinguished from the non-flood-prone areas (Rango, 1975). 

The areas subject to flooding along the Mississippi River were 
identified by observation of various floodplain indicators such as 
natural and artificial levee systems, soil differences, agricul- 
tural pattern and vegetation differences, upland boundaries, back- 
water areas and special flood alleviation measures in urban areas. 



Landsat imagery has also been used to trace the details 
of inundation and drainage of flood areas and deltaic lowlands 
(Burgy , 1973). An overall lightening of tone observed on MSS band 
7 imagery of the Andrus Island flood area in California was attrib- 
uted to an increase of bottom reflection with the lowering of the 
water level. 
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Input signature table 


Sym 

Num 

Delmax 

81 

82 

83 

B4 

Albedo 

W2 

1 

1172 

0.150 

0.58 

0.33 

0.58 

1.51 

3.00 

0.20 

2 

2565 

0.200 

0.55 

0 24 

0.14 

0.13 

1.06 

0.20 



Actual signature table 






81 

J32 

B3_ 

84 

Albedo 





0,74 

0.42 

0.69 

1.77 

3.62 





0.60 

0.25 

0.15 

0.12 

1*12 
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No. points classified s 3737 No. points unclassified * 1383 


Figure 14 Computer classification printout of wetlands for the 
Merrimack River estuary (image ID 544U-l4082). 
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Flood inundation mapping has been accomplished on MSS band 
7 imagery enlarged to a scale of 1:250000 acquired one week to ten 
days after a flood (Hallberg, et al, 1973; Rango and Salomonson, 

1974; Morrison and Cooley, 1973; Schwertz, et al, 1976). The inun- 
dated areas showed sharply reduced infrared reflectance on MSS band 
7 because of standing surface water, excessive soil moisture and 
stressed vegetation (Hallbert, et al, 1973). These data compared 
favorably to flood extent mapping accomplished on low altitude air- 
craft photography. Also, areas affected by severe sand and gravel 
erosion and sediment deposition were detected on MSS band 5 
(Morrison and Cooley, 1973). 

Color enhancement techniques were used to produce a variety 
of multispectral color composites at a scale of 1:250000 of flood- 
ing along the Mississippi River in the spring of 1973 (Deutsch, et 
al, 1973; Deutsch and Ruggles, 1974). Two color composites of MSS 
bands 6 and 7 were enlarged and registered to 1:250000 scale topo- 
graphic maps and used as the data base for preparation of flood 
image maps., Specially filtered three-color composites of MSS bands 
5, 6 and 7 and 4, 5 and 7 were used to aid in the data interpreta- 
tion. In addition, two-color temporal composites of pre-flood and 
post-flood MSS band 7 images were used in interpretation. These 
indicated that flooding caused changes in surface reflectance char- 
acteristics, making it possible to map the flooded areas after the 
flood waters had receded. Also, Landsat MSS data have been digit- 
ally processed to produce water distribution maps and map overlays 
chat show the areal extent of flooding during the 1973 flood of the 
lower Mississippi River (Williamson 1974, 1975). 

Landsat MSS band 7 digital data were superimposed on a 
digital image display and manipulation system (IDAMS) developed at 
Goddard Space Flight Center (Rango and Anderson, 1974) . This en- 
abled a quantitative change detection analysis for determining 
normal surface water area and areas susceptible to flooding. Also, 
the General Electric Multispectral Information Extraction System 
(GEMS) was used to classify and measure water areas according to 
differences in reflectance resulting from physical differences in 
depth and/or sediment load (Rango and Anderson, 1974) . Therefore, 
in general, preliminary digital Landsat flood and flood-prone 
area maps have been produced at a scale of 1:62500; however, most 
mapping has been done on a regional basis at a scale of 1:250000. 

10.2 Approach 

Franklin Falls reservoir, New Hampshire (figure 6) was 
selected for the flooded areas mapping analysis based on a particular 
storm. During the last four days of June 1973, a strong, moist 
tropical airflow in conjunction with a stationary frontal system 
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resulted in moderate to heavy rain over much of New England, For 
example, the total rainfall was 5.1 inches (13,0 cm) for the three- 
day period ending at 0800 hours 1 July 1973 in the Franklin Falls 
reservoir. In the northern portions of the Merrimack River basin 
this storm caused the largest summer flood on record. Sixty-six 
percent of the storage capacity at Franklin Falls was utilized in 
controlling the flood waters. Large areas were inundated for 
periods of one to two weeks (McKim, et al, 1975a). This flood was 
unusual because of its magnitude, the extremely high concentration 
of suspended sediment in the flood waters, and the fact that it 
occurred at the height of the growing season. 

On 6 July 1973 a Landsat pass occurred over the New 
Hampshire area at peak flood conditions, hue to partial cloud cover 
the entire surface area of the flood waters could not be accurately 
delineated; however, the satellite imagery did provide a look at 
most of the peak flood conditions in the Franklin Falls reservoir. 

10.3 Results and Discussion 

The areal extent of water was best displayed in the near 
infrared band, MSS 7. Therefore, MSS band 7 grayscale printouts at 
a scale of 1:24000 were obtained of the Franklin Falls reservoir 
area for the dates of 27 October 1972 (image ID 1096-15065) at 
low-water stage and 6 July 1973 (image ID 1348-15064) during flood 
stage. The MSS band 7 grayscale overprint symbols representing 
energy intensity levels were used in differentiating the reservoir 
water from the land. 

Figure 15 shows the Franklin Falls reservoir area on 
27 October 1972. The extent of water is delineated on the gray- 
scale printout by a solid line. The dotted line indicates the out- 
line of the Pemigewasset River (elevation ranging from 320 to 360 
feet (98 to 110 m)) obtained from TJSGS topographic maps. Both out- 
lines show very good agreement in total area of water depicted at 
the low-water reservoir stage. 



Figure 16 shows the Franklin Falls area on 6 July 1973. 
The extent of water is delineated by a solid line on the grayscale 
printout and the dashed line indicates the extent of water in the 
northern part of the reservoir area, which had to be estimated 
because of cloud cover. The dotted line shows the maximum inun- 
dation level of 369.5 feet (112.6 m) for 6 July 1973 delineated 
from USGS topographic maps. Again, the agreement of the areal ex- 
tent of water between the computer printout and the ground truth 
data is extremely good. 
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The acres of water were quantified for the 27 October and 
6 July Landsat scenes. Table 8 shows the total number of pixels 
and acreage determined to be water within the reservoir area for 
both dates. It shows approximately 60 percent more water in the 
Franklin Falls reservoir on 6 July 1973 than on 27 October 1972, 
which was as expected (McKim, et al, 1975a). 

TABLE 8 

AREAL EXTENT OF WATER WITHIN THE FRANKLIN FALLS 
RESERVOIR AREA ON 27 OCTOBER 1972 AND 6 JULY 1973 


Date 


Pixels 

(#) 


Areal Extent 
(acres)** 


27 Oct 1972 790 

6 Jul 1973 1,333* 


1,011.2 

1,706.2 


*Estimated 
**1 acre = 4,046 m 2 


11 • 0 DCP SENSOR INTERFACE DEVELOPMENT 

11.1 Snow Pillows 


During the summer of 1975 two snow pillow transducer sys- 
tems were interfaced to a General Electric DCP and installed at 
Ninemile Bridge on the Saint John River and at Michaud Farms on 
the Allagash River. A circuit diagram of the interface system is 
shown in figure 17. 


The computer program used for decoding the data is shown 
in Appendix A. A graph of the water equivalent data from these two 
sites during the 1975-1976 winter season is shown in figure 18. 

The sudden increase in water equivalent around April for the Nine- 
mile site cannot be explained and is probably spurious. 

The Bournes transducers used in the interface were tested 
under controlled temperature and pressure conditions during the 
summer of 1976. The resulting temperature calibration curve for 
the Bournes transducers indicated that the system became erratic 
below 0°C (32°F). Therefore, a CRREL in-house study on the relia- 
bility of a number of transducers was initiated and an Endevco 
transducer was used to replace the Bournes transducer in the snow 
pillow interface. 
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1. 


I 






I 1 


Pressure 

Transducer 

CJ.C. Type 4000 

5kn 


Battery 6 volts 


'Ch 6 | 


Black & White or Green & Yellow 




I l 

DCP 

Note 1 ; On some OCR's, data gate is not wired 
to pin 13, J4; in which case connect 
data gate lead to pin 8. J3. 


a. INTERCONNECTION DIAGRAM — Snow pillow transducer to interface to DCP 
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b. CIRCUIT DIAGRAM - interface 


Figure 17 Circuit diagram of the snow pillow interface used in the 1975-76 winter season, 




















The snow pillow transducer system using the Endevco model 
was used during the 1976-1977 winter season. A circuit diagram of 
this interface is shown in figure 19. Two snow pillow interface 
installations were emplaced, one at the Allagash Falls location in 
northern Maine and one at NED, Waltham, Massachusetts. However, 
incorrect data were obtained from these two systems during the 
1976-1977 season. The problems encountered were inadvertent break- 
age of the transducers during handling and shipping, and unex- 
plained, questionable data telemetered by the transmitter within 
the DCP. 

11.2 Water Quality Monito r 

A DCP was installed on the Saint John River in northern 
Maine at the Dickey Bridge, approximately one mile upstream of the 
confluence of the Saint John and Allagash Rivers. A Martek water 
quality monitor interfaced to the DCP transmitted the following 
water quality information: pH, dissolved oxygen, river stage, 

water temperature and conductivity. The sensor interface develop- 
ment and the computer program used to decode the data were accom- 
plished during the Landsat-1 experiment (Cooper, et al, 1975; 

McKim, et al, 1975c). 

The water quality data from this exercise are shown in 
figure 20. The dissolved oxygen probe operated intermittently and 
the river stage measured less than 2 feet (62 cm); therefore, these 
data were not included in figure 20. The pH between 10-12 August 
probably did not fluctuate as indicated on the graph. The water 
quality data compared favorably with onsite analysis of these para- 
meters during the first week of operation. The water quality in- 
formation will serve as part of the baseline data for the upper 
Saint John River. 

11.3 Thermocouple Interface 

A site for testing a thermocouple interface to monitor 
air and ground temperatures was located at Sugarloaf Mountain, 

Maine (45°02 f 56 ft N 70°23 I 21 |, W) . This was the first time the ther- 
mocouple interface was tested under field conditions. The data are 
presently being evaluated for accuracy. The emplacement and inter- 
facing techniques developed during this field experiment will be 
used for installation of thermocouples in Alaska. When validated 
this system could be used in reservoirs to monitor water tempera- 
ture at various depths on a daily basis. 

There is evidence of permafrost on the upper 1,000 feet 
(305 m) of Sugarloaf Mountain. The summit of the mountain is 
veneered with active, turf -banked terraces which have moved down- 
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slope during the past five years at a rate of approximately 2 
inches/year (5.1 cm/yr) (Borns, 1075 ). When the gondola ski lift 
was installed in 1967, the subsurface was frozen in August and 
the water in the drill holes had to be melted before the tower 
foundations could be installed. In addition, the two towers at 
the top of the mountain moved in the time period from 1967-1976 
and had to be realigned. The deep pit at the summit that houses 
the counter-weight used in tensing the gondola cable contains 
water that is frozen year-round. The objectives of this study 
were to evaluate the reliability of the thermocouple interface 
under field conditions, to monitor the ground temperature measure- 
ment to a depth of 100 feet (30 m) , and determine if permafrost 
conditions occur year-round at this latitude and elevation. 

A complete circuit description of the thermocouple inter- 
face is illustrated in Appendix B. The thermocouple interface 
accepts inputs from 112 copper-constantan thermocouples. The 
inputs are arranged in 16 "banks" of seven inputs each. One bank 
of inputs is recorded with each update of the DCP. A total of 
28 thermocouples are being monitored at Sugarloaf Mountain with 
seven readings being transmitted in one Landsat message (table 9) . 

TABLE 9 

ARRANGEMENT OF THE 28 THERMOCOUPLES IN FOUR BANKS 


Bank 

Thermocouple No. 

Depth 
(ft) * 

i 

1,5,9,13,17,21,25 

0,8,24,40,56,72,88 

ii 

2,6,10,14,18,22,26 

0,12,28,44,60,76,92 

hi 

3,7,11,15,19,23,27 

0,16,32,48,64,80,96 

IV 

4,8,12,16,20,24,28 

4,20,36,52,68,84,100 


*1 ft = .3 m 


The temperature measurement range of the thermocouple 
interface unit is -34 to +32°C (-29 to 90°F) . The resolution of 
a DCP data word is + 0. 25°C (+0.5°F) or lOyv and the copper- 
constantan wire is guaranteed to be accurate within 0.5 to 0.75°C 
(0.9 to 1.35°F). Therefore, the accuracy of the temperature data 
is. + 0.5°C (0.9°F). The results from this experiment will be re- 
ported at a later date. 

11.4 Tensiometer/Transducer System Interface 

A soil moisture tensiometer/transducer system has been 
successfully interfaced to the Landsat data collection system 
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(McKim, et al, 1975a), The interface system enables moisture 
tension and soil volumetric moisture content data to be obtained 
in near-real time. 

The instrument is currently being tested under simulated 
field conditions at CRREL. Preliminary results indicate that 
tension as low as 10 cm (4 inches) and as high as 900 cm (354 in) 
of water can be obtained. 

Typical data for a soil that has a bulk density of 1.37 
g/cm 3 with a specific gravity of 2,63 and a volume of voids of 
47.9 percent are shown in table 10, In the initial series of tests, 
tensiometer values could be obtained that ranged from about 10 to 
150 cm (4 to 59 in) of water. Previously it has been very diffi- 
cult to get accurate and reliable numbers on field moisture ten- 
sion less than 300 cm (118 in). It is suggested that this method 
will not only give reliable data for this range of values, but also 
supply the information in near-real time. 

The precision and accuracy of the data are being evaluated. 
The results will be reported at a later time after the first field 
operations have been completed. 

12.0 CONCLUSIONS ON IMAGERY AND SENSORS 


Preliminary analysis of the Landsat digital data using 
the GISS computer algorithms for the 11 February 1973 scene showed 
that the radiance of the snow cover /vegetation varied from approxi- 
mately 20 mW/cm 3 sr in non-vegetated areas to less than 4 mW/cm^ 
sr for densely covered forested areas. Comparison of the digital 
data from three snow course sites in the Dickey-Lincoln School 
Lakes area with the radiance value of the snowpack at these sites 
indicated that the total radiance of the pixels contained in the 
snow courses varied from 5,34 to 7.74 mW/cm^ sr with the average 
radiance value being 6,4 + 0.6 mW/cm^ sr. The water equivalent 
of the snowpack for this range of radiance values was approximately 
9.5 (24.1 cm) of water. Since data from only three snow courses 
were available, the correlation between radiance values and water 
equivalent of the snowpack still needs to be drawn. However, if 
the relationship holds with more extensive ground truth data in 
the Dickey-Lincoln area, it is anticipated that extrapolation of 
radiance values for snow cover/vegetation in large areas of the 
watershed may prove useful in estimating water equivalent. This 
method, after tests for accuracy, would find use in estimating the 
volume of spring runoff. 
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TABLE 10 


TYPICAL DATA FROM THE 
TENS IOMETER/ TRANSDUCER INTERFACE SYSTEM 

USA CRREL Test Set 
Hanover, New Hampshire 
DCP 7110 


Calculations 


Date 

Time 

Tensiometer 

% Water 

% Voids filled 
with Air 



(cm of water) 

(volume) 


2/22 

1752 

-10.7 

40.0 

17.0 

2/23 

622 

-10.7 

40.0 

17.0 

2/23 

1801 

-22.5 

40.0 

17.0 

2/24 

628 

-22.5 

40.0 

17.0 

2/24 

1804 

-22.5 

40.0 

17.0 

2/25 

633 

-28.3 

39.5 

17.5 

2/25 

814 

-22.5 

40.0 

17.0 

2/26 

638 

-28.3 

39.5 

17.5 

2/26 

820 

-28.3 

39.5 

17.5 

2/27 

644 

-34.2 

38.5 

19.7 

2/27 

829 

-34.2 

38.5 

19.7 

3/1 

1832 

-40.1 

37.8 

21.1 

3/2 

702 

-40.1 

37.8 

21.1 

3/2 

844 

-40.1 

37.8 

21. r 

3/3 

1703 

-51.9 

36.0 

24.8 

3/3 

2030 

-57.8 

35.4 

26.1 

3/4 

856 

-57.8 : 

35.4 

26.1 

3/6 

1721 

-81.3 

30.2 

37.0 

3/6 

1859 

-81.3 

30.2 

37.0 

3/7 

911 

-87.2 

29.2 

39.0 

3/8 

735 

-93.1 

28.8 

39.9 

3/8 

924 

-99.0 

27.9 

41.8 

3/9 

741 

-110.7 

26.6 

44.5 

3/9 

924 

-110.7 

64 

26.6 

44.5 
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The different changes in radiance for the 23 July 1973 
(3 mW/cm2 sr) and 11 February 1973 (5.34 - 7.74 mW/em^ sr) Landsat 
scenes for the three snow course sites were attributed to the 
presence of snow and vegetation changes in the 11 February scene. 
Multispectral training signatures for four snow cover /vegetation 
classes were derived from the four-band energy values for the 
11 February 1973 Landsat scene. These signatures were applied to 
the digital data and four snow cover/vegetation classes were 
mapped. 


Landsat digital data were also used to delineate flood 
waters in the Franklin Falls reservoir area, New Hampshire, for the 
6 July 1973 scene. Low-water reservoir and flood water stages were 
mapped from grayscale printouts of MSS band 7 for 27 October 1972 
and 6 July 1973, respectively. Comparison with ground truth in- 
formation indicated very good agreement with the Landsat digital 
data. Approximately 60 percent more water was observed on the 
6 July 1973 scene than on 27 October 1972, which was as expected. 

In addition, Landsat! digital data were used in a wetlands 
analysis of the Merrimack River estuary for 6 July 1576. A multi- 
spectra! signature was developed for a x^etlands category from two 
sites on the north side of the Merrimack River estuary. Wetlands 
were mapped with an accuracy of 75 percent compared to ground truth 
information. 

Snow pillow transducer systems for measuring the water 
equivalent of the snowpack in northern Maine were interfaced and 
field tested. Little valid data was transmitted from the field, 
because problems with the transducers were encountered during both 
winter test periods. 


A thermocouple system was successfully interfaced and field 
tested at Sugar loaf Mountain, Maine. Temperature data from the 
surface to a depth of 100 feet (30 m) were transmitted through the 
Landsat DCS. The emplacement and interfacing techniques developed 
during this experiment will be used for installation of thermo- 
couples at Alaskan sites. In addition, this system could be used 
in reservoirs to monitor water temperature on a daily basis. 

A soil moisture tensiometer/transducer system was success- 
fully interfaced to the Landsat DCS. Laboratory results indicated 
that tension as low as 10 cm (4 in) and as high as 900 cm (354 in) 
of water can be obtained. 'Presently the system is being tested 
under field conditions. 




A water quality monitor interfaced to the DCS was also 
field tested in northern Maine* The water quality data compared 
favorably to onsite chemical analysis and will be used as baseline 
information for the Dickey-Lincoln School Lakes Project, 

13.0 RECOMMENDATIONS 


The following recommendations are made based on the result s 
of the Landsat-2 investigation. 

1. Further research work should be aimed at relating a 
wide range of ground truth sets to multispectral signature sets. 

It is clear that a method of quantifying the water equivalent of 
snow cannot be perfected with a small number of samples. Landsat 
CCT's should be furnished on a continuing basis for research in 
evaluating the snow cover of a watershed, especially during crit- 
ical spring runoff periods (March through May) . Nine day coverage 
provided by the Landsat satellites is not adequate for the opera- 
tional needs of the NED Reservoir Control Center, 

3. When a mapping accuracy of 75 percent Is required, it 
is recommended that the Landsat digital data be used in regional 
evaluation of wetlands. 

4. It is clear from this investigation that the Landsat 
Data Collection System is suitable for monitoring events for which 
dependable sensors are available. With more evenly timed coverage 
over each day, a system of polar-orbiting satellites would be 
operationally useful for collection of hydrologic data such as 
river stage and precipitation. 

14 . 0 FUTURE PLANS 

The snow cover analysis work will be continued by CRREL 
with the FY 78 work unit entitled, r, Snow Cover Analysis in New 
England Using Landsat Digital Data”. Site selection based on vege- 
tation, slope, aspect and elevation will be accomplished in the 
Dickey-Lincoln School Lakes project area and other selected water- 
sheds. Ground truth measurements of snow depth and water equiva- 
lent at the selected sites will be taken in conjunction with times 
of the Landsat imagery acquisition. A meteorological station will 
be installed to obtain data on local climatic conditions. If 
available, cloud-free Landsat CCT’s will be acquired for the 1977- 
1978 winter season and analyzed. 
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Cloud-free Landsat CCT’s will also he obtained over the 
Sleepers River watershed in Danville, Vermont for the 1972-1977 
winter seasons. Detailed measurements of the snow cover are 
available at this site from December 1968 to present. This water 
shed was chosen because it is hydrologically representative of 
most of the glaciated upland areas in New England and is exten- 
sively instrumented. 
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GLOSSARY 


ADR - Analog to digital recorder. Typically a Fisher-Porter or 
Leupold-Stevens recorder, often equipped with a telekit 
for interfacing to DCP f s. 

AHRRN - Automatic Hydrologic Radio Reporting Network, NED’s 
ground-based radio system. 

Albedo - Refers to the fraction of incident radiation which is 

reflected by a material, summed over all wavelengths for 
a sunlit surface. 

ASCII - Abbreviation for USA Standard Code for Information Inter- 
change. An eight-bit character code. 

Azimuth - Horizontal angle measured clockwise from north. 

BCD - Binary Coded Decimal. 

CCT — Computer compatible tape. 

DCP - Data Collection Platform. Field installation used for 
sensing parameters, encoding data, and transmitting 
data to satellite. 

DCS - Data Collection System. Refers to configuration including 
field stations, medium of transmission, and central 
station. 

Epoch - A particular instant. 

Julian Date - A count of days and fraction of days since a par- 
ticular reference instant. 

Micrometer - 10~^ meters; used in wavelength measure. Symbol: ym. 

Multitasking - In a computer, several program tasks competing for 
devices and the central processor on a priority or 
queued basis. 

NORADC - North American Air Defense Command. 

Octal - Refers to a number system having 8 as a base. 

Pixel - Picture element, the unit area of Landsat scenes. 
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Real Time Clock - Device in the Data General NOVA computer that 

consists of a crystal controlled clock and associated DG 
system software that are used (1) to keep track of date 
and time of day and (2) to provide for low resolution 
timing. 

Sidereal Time - The relationship between an observer's meridian 
and some inertial coordinate system, for example, one 
based on the constellation Aries. 

Stage - Water level. 

Time Code Generator - A device which keeps standard time and 
outputs it in a form suitable for input to another 
device. , 

Tracking - Keeping the antenna pointed at the satellite, and in 
conjunction with that, logging any incoming data. 
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APPENDIX A 


PROGRAMS AND OPERATING PROCEDURES USED IN GROUND RECEIVE STATION 

Preparation for the ground receive station began over two 
years before the pedestal installation in October 1975. After firm 
commitments to construct the downlink were made in 1973, site 
preparation, hardware fabrication, and computer programming began 
late in 1974. A program was written by NED to predict satellite 
passes and a rudimentary version of a program to aim the antenna 
and simultaneously log incoming data was readied by June 1975. 

All programs (figure A-l) were written in-house, because previous 
experience had shown that complex programs written by outside con- 
sultants require extensive modification when installed and that 
in-house personnel may not know the groundwork and assumptions 
that go into them. 

PREDICT 

The program that predicts satellite passes over NED, "PREDICT" 
was the first one undertaken. While it was known that the problem 
of predicting the time qf a satellite’s rise over one’s horizon 
and its path across the sky had been solved before, its solution 
was not available in a form that was compatible with NED’s mini- 
computer, a Data General Nova 1220 with 32 thousand 16-bit words 
of core storage. (The Nova had been selected several years before 
for use with the Automatic Hydrologic Radio Reporting Network.) 

A prediction program was written in Fortran, because it is readily 
understood by in-house programmers and lends itself to later mod- 
ifications. At that early stage of programming , assistance was 
provided by individuals at NASA Goddard Space Flight Center to 
quickly assemble the most useful algorithms and background mater- 
ial for writing PREDICT. 

The task of debugging the combined routines that comprised 
PREDICT in the early stages was facilitated by the process of 
"units analysis", i.e., the terms in the right side of each line 
of Fortran were assigned appropriate engineering units, and sim- 
plified by cancellation wherever possible. The simplified ex- 
pression of units was then associated with the single variable 
on the left side of the Fortran statement, and those units were 
employed in later statements when the same variable occurred. By 
checking the units on paper for the entire program (which ran to 
about 300 lines) down to the final result, the likelihood of 
immediately getting correct numerical answers was increased greatly. 
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PREDICT starts with inertial orbital elements and calculates 
station range vectors (azimuth and elevation) at selected time 
increments. Further program steps select only those vectors that 
are above the horizon. The program is generalized for any sta- 
tion and for any satellite in a nonequatorial orbit. For an aid 
to understanding the discussion of PREDICT see the flowchart in 
figure A-2 and the program listing. 

After the current date and time are read by PREDICT, they are 
converted to sidereal time by the subroutine, TCALC. Sidereal 
time is an angular relationship between the Greenwich prime meri- 
dian and an inertial X-axis which points toward the first point 
of the constellation Aries (Escobal, 1965). This angle is denoted 
by THETA (0), local sidereal time. Knowing the east longitude 
(Ae) of an observer's station and 0g (Greenwich Sidereal time), 

0 can be easily determined by the relationship 0 = 0g - Ae, where 
0 < 0 < 2tt (Escobal). 

To find the sidereal time, the Julian Date must also be 
calculated. The Julian Date is continuing count of each day 
elapsed since some arbitrarily selected epoch. (The epoch selected 
for Landsat orbital predictions in TCALC is January 1, 4713 B.C.) 
Each Julian day is measured from noon to neon; hence, it is an 
integer 12 hours after every midnight (Escobal, 1965). After the 
Julian Date and sidereal time are found, TSINCE, the number of days 
since the most recent NORADC EPOCH is calculated. TSINCE is then 
used to determine the unit vector pointing toward the satellite 
(see SGP of this text). This unit vector is converted to azimuth 
and elevation angles at the observer's station (degrees clockwise 
from north and degrees above the horizon, respectively - see SRV 
of this text) . These two values are then written on the disc 
along with the time at which they occur in the file named PTAE. 

The _Time of interest is then inc remented seconds or minutes by 
the routine INCT, and the next values of azimuth and elevation 
are determined, etc. In this incremental fashion, the computer 
is able to predict the path of the satellite across the sky at 
the observer's site. 

The file name PTAE stands for Paper Tape Azimuth .Elevation. 

On rare occasions, for diagnostic purposes, it may be necessary 
to transfer the file PTAE to the paper tape punch to generate a 
paper tape that is suitable for input to the paper tape reader on 
the pedestal control equipment. PTAE is a disc file which com- 
prises a time, a time increment, many pairs of azimuth and eleva- 
tion angles, and a special file terminator. A sample of PTAE is 
shown in figure A-3. 
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type ptae 

281.979629630 

10.0 

A 95.S9E 0.22 J 

A 94.9SE 0.61* 

A 94.0SE 1.001 

A 93.09E 1.391 

A 92.10E 1.771 

A 91.09E 2.161 

A 90.06E 2.541 

A 89.00E 2.921 

A 87.91E 3.291 

A 86.80E 3.661 

A S5.66E 4.021 

A S4.49E 4.381 

A 83.29E 4.731 

A 82.06E 5.081 

A 80.81E 5.411 

A 79.53E 5.741 


JULIAN DATE OF SATELLITE 
RISE TIME MEASURED FROM 

BEGINNING OF CALENDAR YEAR 

o DT. TIME INCREMENT (SECS) 

BETWEEN THE FOLLOWING 
ANGLE PAIRS 


PAIRS OF AZIMUTH AND 

ELEVATION ANGLES 


A 57.43E 8.981 
A S5.82E 9.05! 
A 54.21E 9.10! 
A 52.59E 9.12! 
A S0.92E 9.13! 


A 16.83E 3.81! 

A 15.72E 3.45! 

A 14.64E 3.08! 

A 13.59E 2.70! 

A 12.56E 2.32t 

A 11.56E 1.94! 

A 10.58E 1.56! 

A 9.63E 1.17! 

A 8.71E 0.78! 

A 7.80E 0.39! 

A 6.93E 0.00! 

A 0.00E 90.00! -=> 


PEAK ELEVATION 
FOR THIS PASS 


PAIR OF STOW ANGLES: 
SENDS ANTENNA TO 
STOW POSITION 


ORIGINAL PA.GEIS 

OF POOR QUALITY' 


Figure A-3. Format of Film "PTAE*. 


OEDIT BCDAZEL 


A 

0 

1 

£ 

3 

4 

5 

6 
7 

10 


28 
- 1 . 
>•97 
'96 
>29 
'63 
>0 
/ 

/I . 


11 /0<015> 
0 

12 /0002 

13 /Q00225 

14 /00O042 

15 /0000 

16 /00017T 

17 /0002E5 

20 /O00045 

21 /0000O0 

22 /000160 

23 /00O225 

24 /00OO51 

25 /00O000 

26 /000141 

27 /00022S 

30 /000O63 

31 /000O00 

32 / 000 122 

33 /0O0225 

34 /0O0067 

35 /0O0000 

36 /0OO103 

37 ' 000225 
40 /O0O101 
HOME 



LOW ORDER AZIMUTH: 

.8 .08 . 0 ! 

211 = 010 001001 - .89 
8 2 10 

HIGH ORDER AZIMUTH: 
80 10 4 I 

255 = 010 010 101 = 95 
8 2 10 

LOW ORDER ELEVATION: 

.2 .02 

042 ~ 000100010 = .22 
8 2 tO 

HIGH ORDER ELEVATION 

000 - 0 
8 10 


I 


95.09 


0.22 


Figure A-4. Format of File "BCDAZEL*. 
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In practice, PREDICT calculates azimuths and elevations at 
10-second increments, so the angle pairs in PTAE are pointing 
angles for instants ten seconds apart. If these angles were fed 
to the pedestal, the antenna would jump quickly to the next posi- 
tion every ten seconds. The progress of the satellite is uniform, 
and it has been found that one-second incrementing of dish posi- 
tion is sufficiently small for continuous radio lock. Therefore, 
the program INTERPl is executed right after PREDICT to interpolate 
ten angle pairs for every one pair in PTAE. Furthermore, INTERPL 
recodes the angles from ASCII characters to a binary coded deci- 
mal (BCD) format suitable for the electronic interface enroute to 
the pedestal command equipment. The new angles are stored in a 
binary file called BCDAZEL (figure A-4) and the number of pairs of 
angles in BCDAZEL is stored in the file NANGLES . 

The recoding in the program INTERPl is done by bit-mapping. 

A set of special interfaces is used to route the BCD angle informal 
tion from the Nova to the Scientific Atlanta 1848 Digital Compara- 
tor (see figure 1). Input to the comparator is in the form of 
two 18-bit BCD words representing azimuth and elevation. However, 
the Nova can output only 16-bit words, by way of the 4065 Digital 
Interface. For this reason, it was necessary to concatenate two 
pairs of 16-bit words into two 18-bit words in the sequence shown 
in figures A- 5 and A-6. 

TCALC 


TCALC is a time handling subroutine called by PREDICT which 
calculates the following three variables: 

1. XJD - Number of Julian days - it is used to find 
TSINCE and THETA. 


2. TSINCE - Number of days since most recent NORADC 
EPOCH. This is used by SGP to find the unit vector pointing 
toward the satellite (RDOT) . 


3. THETA, 0_ - Sidereal time (measured in radians). The 
angle between a line from the center of the earth L& the first 
point of the constellation Aries (T) , and the plant of observer’s 
meridian. . » 


ObMrvtr'* mtridlan 




Ctntar of aarth 
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NED, Waltham 
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18 -BIT SPECIAL 16 - BIT 

WORDS INTERFACES WORDS IN 

OUTPUT | I THE NOV/3 






START 

PULSE 


START 

PULSE 


LOW 

ORDER 

AZIMUTH 


START 

PULSE 


START 

PULSE 


HIGH 

ORDER 

AZIMUTH 


AZIMUTH 


LOW 

ORDER 

ELEVATION 


HIGH 

ORDER 

ELEVATION 



TIME (/v/ m s) 


Figure A-5. Order of Transmission of Angular 
Information from Minicomputer. 
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THETA is used by SRV in determining Azimuth (A) and elevation 

(H). 

The following calculation (based on Escobal, 1965) derives 
sidereal time (THETA) from the following instant expressed as a 
date and Greenwich mean time: 

August 23, 1975, at 10 hours, 15 minutes, 0 seconds. 
Number of hours, minutes, seconds expressed in minutes: DT = 

615 minutes. 

XJD = 2442648.5 DTHDT - .25068447 

TU = (XJD-2415020) /36525 = .7564271047227926 
THETA G0 = 99.6909833 + (36000.7689) (TU) + (.00038708) 
(TU)2 = 331.6485916015490 

THETA G = THETAG0 + (DT) (DTHDT) = 125.8195406515490° 
THETA * (THETA G + LAMBDAE) (2tt/360) = 

.9530184529430946 radians. 


SGP AND SRV 


SGP is a Fortran subroutine called by PREDICT embodying a 
truncated simplified general £erturbation theory for use in the 
determination of Landsat pointing elements, SGP computes osculat- 
ing position, velocity and mean classical elements. SGP is a 
first order analytical integration of the equations of motion in- 
cluding perturbations caused by the first two zonal harmonics of 
the geopotential. The zonal harmonic constants account for the 
effects of the noncircularity of the meridian cross sections of 
the earth. The perturbations caused by these harmonics are inde- 
pendent of the longitude of the satellite, SGP is based on the 
orbital elements a, A xn , Ay n , i, ft, and L which are well defined 
for all elliptic orbits except those that are nearly equatorial. 

For equatorial satellites, the elements A xn and Ay n are ill- 
defined because of the indeterminacy of the node angle 0 . to 
which they are referred. The SGP mathematical model is adequate 
to handle a majority of routine cataloguing. Accuracy is said 
to be better than one part in 10^. 

SRV ("Slant Range Vector”) is a Fortran subroutine of PREDICT 
which transforms the orthogonal vectors and the time angle, 

THETA, from subroutine SGP into an azimuth/ elevation coordinate 
system with the observer’s station as the origin. Files of azimuth 
and elevation angles in this coordinate system describe the path 
of Landsat over a particular station during some interval. 
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TRACKING THE SATELLITE: PROGRAM 11 TRACK* 1 


After PREDICT calculates the satellite’s sky path and INTERP1 
prepares a file of pointing angles, control goes automatically to 
the program TRACK which will perform six main functions. It is a 
complex multitasking program in which internal control shifts ac- 
cording to time as counted down by the Real Time Disc Operating 
System (RDOS) and according to the random arrival of DCP messages. 

TRACK carries out the following main tasks: 


1. Schedules itself by looking at the starting time of 
the upcoming pass. This time is the first number stored in the 
disc file BCDAZEL. 


2. Orients the antenna 1-1/4 minutes before the satellite 


rises, 


3. Starts repositioning the antenna second by second 
beginning at the instant the satellite rises. 

4. Logs any data that arrive by way of the antenna/ 
receiver/decoder pathway (see figure 1). 

5. Accepts corrections from the terminal and console 
switches to advance or retard the tracking antenna by a certain 
number of seconds to Improve antenna position. 

6. Restores the antenna to the stow (upright) position 
when the last angle pair in file BCDAZEL has been sent. 

7. Transfers the DCP data that have come in from a 
core buffer to a temporary disc file called ,, SDF M (Satellite Data 
File) . 

8. Finally chains to a program called QD4 which will 
decode field data from binary to an octal format similar to one 
used by NASA at Goddard . 

Once PREDICT and INTER? 1 have been executed for an upcoming 
pass, TRACK can be run at any time up to one minute 30 seconds 
before satellite rise time. Execution of TRACK after that causes 
problems which are signalled by the message "TOO LATE 1 ' being printed 
at the terminal. One then must quickly reset the system clock to 
an instant one minute and 30 seconds prior to the satellite rise 
time*; execute TRACK; and six seconds* after the computer types 
"l", enter positive corrections that stand for numbers of seconds 
to enable TRACK to catch up with real time. (If this method has 
to be used, note that incoming data will be incorrectly time- 
tagged.) 


*A typical command would be STOP 14:30:10: TRACK 
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Corrections (in seconds) are entered by means of the NOVA's 
console switches numbered 0-15, and the Tektronix terminal. The 
switches have binary significance in descending order from left to 
right, and bit 0 indicates the sign, + or 


BIT; 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 0N+ 

_± ± 32 16 8 4 2 1 OFF! 


Thus, a correction of +20 would require that only switches 
11 and 13 be up, whereas the correction -20 would have bits 0, 

11 and 13 up. 

To effect the correction set in the switches, strike SPACE 
BAR on the terminal. Positive corrections will cause the antenna 
to jump ahead, and negative corrections will cause it to pause. 

For DCS data to be usable, it must be timely and reduced. 

For this investigation timeliness was guaranteed by use of the 
direct downlink with no delays except insignificant propagation 
time between a DCP's report and our reception of that report. 

The matter of data reduction required further analysis and pro- 
gramming. It was necessary to screen out erroneous messages and 
any unwanted ones which were transmitted by DCP’s which belong 
to other agencies. (In some cases other users’ messages were 
gathered by our system for cooperative use.) Whenever two or 
more successive messages were identical they were compressed into 
one disc entry with the time of only the first transmission and the 
number of identical messages included. At this point the messages 
were stored in ASCII code in a main file on our largest storage 
medium, a twenty surface disc pack. Their format is a mixed 
decimal and octal representation that is a compromise between 
maximum density and understandability . For final reduction to 
readily comprehensible symbols, two or three programs are executed 
to retrieve data and either decode it into legible lines having 
engineering units or rapidly plot it as in the case of hydrographs. 
The overall handling scheme is shown in figure A-7 . 

Data gathered during a satellite pass are reduced and auto- 
matically displayed within two minutes after the pass. Old data 
are retrieved, decoded, and tabulated or plotted in less than two 
minutes also. A sample of plotted and tabulated data for a flood 
hydrograph of an actual flood which occurred at Fort Kent, Maine, 
in April of 1977 is shown in figure A-8. 
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230030 

1S3000 

160000 

140000 

120003 

100000 

80000 

60000 

40000 


77 -42i S-37 7273 36550 
77 4211018 7273 36903 
77 42112 1 7273 37033 
77 4212012 7273 43160 
77 4212156 7273 43640 
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7273 FORT KENT, PIE . 


125 


Figure A-8. Samples of Data Reduction by Plotting 
and Tabulation of a Flood Hydrograph. 
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OUTPUTTING DATA: QD4 AND P3 


QD4 and P3 are programs that condition the raw field data 
received by the ground receive antenna for disc storage or for 
easy reading. 


The temporary disc file "SDF n produced by TRACK immediately 
after a satellite pass becomes the input to QD4. Output from QD4 
goes to a temporary file f, LS2DAT M and a permanent file "STORAGE". 
Formats of these files are shown in figures A-9 and A-10. 

The arrival time of each DCP message is recorded by LSI by 
storing a seconds counter with each message. This number of 
seconds is accumulated from the beginning of each pass, and QD4 
calculates a message arrival time by adding the number of elapsed 
seconds to the starting time. Coordinated Universal Time is 
employed to avoid problems that arise with standard versus day- 
light savings time. 


Legible output of DCS data is obtained by executing P3. 

Input to P3 is from the temporary disc file "LS2DAT". The pro- 
gram examines each message for the platform ID number, and looks 
up the ID in a table (INDX", figure A-10) which contains indices 
concerning the parameters being measured* These indices then 
direct program control to appropriate subroutines for calculating 
floating point decimal numbers and attaching labels. The kinds 
of parameters handled by P3 are shown in the site list, figure 
3. 

STORAGE NEEDS 


The amount of disc storage needed for long term retention 
of Landsat DCP data is small and inexpensive. For parameters that 
are slow to change, such as river stage or precipitation, data 
may be compressed merely by counting successive identical messages 
as one message. Concatenations of message components would per- 
mit further compression of data on the order of 7:1. This has not 
been done at NED. We have found that a compression of 3:1 is 
attainable with no sacrifice of legibility of storage files and 
no highly sophisticated techniques. Hence, for a system consist- 
ing of one DCP reporting hydrologic data randomly at three minute 
intervals and one satellite having an orbit like Landsat 1 s, the 
following calculation is typical: 

ORIGINAL PAGE IS 
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HOUR. 


MINUTE. 


SECOND 


ERROR CODE. 


PLATFORM I.D. 


DATA 


TYPE STORAGE 

.610 1132230 7273 37 65377377377377377377 
610 1132321 7021377377267323377376 57126 
6,10 1132353 7246377163377377377377377377 
610 1 1 32442E7254 67366345313232 64153227 
610 .1 132544 7356 77237377377377377377377 
610 1132746 7142377377 0 0 0 20 0200 

610 11328 3E720 12771 4 1377377 33226215 31 
610 1132847 7273 37 65377377377377377377 
610 1 132853 7356 77237377377377377377377 
610 l‘l5 2 2 7142377377 0 0 0 20 0200 

610 115 232 7207337375377377377377377377 
610 115 , 245 606 3277121377377377377377377 
610 life 253 6504 14341 23 14 0 0 14 6 

610 115 330 7127377 76377377377377377377 
610 115 3 50 7345 77341377377377377377377 
"610 1 15 V 1 7124377377377377377377377377 
610 115 412 7021377377267323377376 57126 
610 115 415 7010160 40 0 0 0 0 0 0 

610 115 428 7 254327137377377377377377377 
610 115 5 ?' ‘7246377163377377377377377377 
610 115 526 7220277223377377377377377377 
610 115 528 7142377377 0 0 0 20 0200 

610 115 539 7272 77265377377377377377377 
610 115 546 6504 14341 23 14 0 0 14 6 

610 115 555 7775337375377377377377377377 
610 115 634 ,7127 377 76377377377377377377 
6ld 1 l5 648 7010160 40 0 0 0 0 0 0 

610 115 7 6 7273 37155377377377377377377 
610 115 7 8 7201277341377377377377377377 
610 115 744 7124 00000000 

610 115 745 7254377377377377377377377377 
610 115 752 7356 77277377377377377377377 
610 115 757 7171 77173377377377377377377 
610 115 810 7246377163377377377377377377 
610 115 822 6063277121377377377377377377 
610 115 833 7220277223377377377377377377 
610 115 354 7142377377 0 0 0 20 0200 


Figure A-9. Formats of File ■STORAGE" or * LS2DAT " . 





NUMBER OF LINES IN THIS FILE 


PLA T FORM 
I.D. NUMBER 



INDEX USED IN ‘P 3' TO DIRECT 
PROGRAM CONTROL 


2 
a 
a 
a 

3 
2 
2 
6 


7071 
7272 
7356 
7127 
7201 
7207 
7331 
7214 1 
7246 3 
6063 
7106 
7142 
7021 
701010 
7220 2 
7345 
7254 
7242 
7206 
6504 
R 


2 

2 

4 

4 


2 

2 

2 

2 

2 


NINEMILE BP.. HE — 
FORT KENT, HE. 

(JEST ENFIELD. HE. 
NORTH ANSON, HE. 
CORNISH, HE. 

SOUTH NT., N.H. 
PLYHOUTH, N.H. 

GOFFS FALLS, N.H. 
CRREL, HANOUER, N.H. 
CRREL, HANOUER, N.H. 
UACHUSETT (IT., HA. 
IPSUICH, HA. 
FITCHBURG, HA. 
CHICOPEE, HA. 
UESTFIELD, HA. 
NED.UALTHAH, HA. 
FOREST DALE, R.I. 
CRANSTON, R.I . 
HARTFORD, CT. 
HIDDLETOUN, CT. 
HANCHESTER, CT. 

ST. FRANCIS R., N.B. 


PLATFORM LOCATION 


Figure A-10. Formet of File * INDX 
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Since these disc packs sell for as little as $60, the cost 
of the storage medium for an entire year of data from one DCP is 
no more than 15c?. 


OPERATING PROCEDURES 


The automatic tracking system for receiving Landsat data at 
New England Division consists of a 15-foot dish antenna, a track- 
ing pedestal, pedestal control equipment, and a Data General Nova 
minicomputer with various accessories and interfaces. The rela- 
tionship of all these parts is shown in the subsystems diagram, 
figure 1. Little knowledge is required to place the system in 
operation, because the programs have been written to guide the 
operator through the various options using either the Tektronix 
or Teletype terminal (figure A-13) . Full control of the system 
to handle unusual situations requires some knowledge of the pro- 
grams and various information files that are kept on disc. The 
operator may have to power up the Data General Nova minicomputer 
and execute programs to track the Landsat and store incoming data. 
Having the correct time of day and executing the tracking pro- 
grams, the Nova should continue until interrupted by the operator. 
The programs are cyclical, and if one of them is int erupted, it 
may be restarted later; that is, the operator may re-enter the 
cycle at several points. 

The simplest procedure for tracking is as follows (refer to 
figure A-l) : 

1. Power up all equipment. 

2. Execute the program PREDICT at least five minutes 
before an expected satellite rise time. This is accomplished by 
typing "PREDICT" followed by a carriage return. 

REAL TIME CLOCK 

To track Landsat accurately the Nova T s real time clock must 
be set to Coordinated Universal Time (UTC) . Accuracy of one- 
fourth second is sufficient. Two methods may be used - either 
manual or automatic. 
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Automatic Method . Execute the program CLOCK . Within two 
minutes the computer will signal completion by a message, and it 
will progress to the next program. If problems arise, use the 
manual method. 

In the automatic method the computer uses an interface which 
samples (at 550 Hz) the serial code output of a time code genera- 
tor. The time code generator (made by DATUM, Inc.) outputs a 
code which contains day, hour, minute, and second information 
within the frame of one minute. 

Theory of Operation . The pulses in this code consist of 0 
to 5 volt shifts for varying lengths of time. A positive shift 
lasting .2 seconds signifies a logical "0", a shift for .5 
seconds signifies a ,, 1 M , and a shift for .8 seconds signifies a 
special mark pulse. The computer samples this shifting voltage 
a 550 Hz; thus for a logical ,, 0” the voltage would be "high” for 
.2 x 550 or about 110 times; for a "1” voltage would be high for 
.5 x 550 or about 175 times, and for a mark pulse the voltage is 
high for .8 x 550 or about 440 times. The string of 0's, 1 1 s and 
mark pulses are interpreted in a tabular fashion, made possible 
because their positions within one minute determine their values. 

Manual Method . The Nova provides for its clock to be set 
by the teletype command STOP hh mm ss , where hh, mm, and ss 
stand for hour, minute and second. 

Dial up the FTS number 8-323-4245 to get the National Bureau 
of Standards audio time signal. When you have found out what time 
mark will be coming up soon (e.g., the next minute), use the STOD 
command to prepare to enter that upcoming time, and hit RETURN 
exactly when the time marker is heard. 

Setting the Time Code Generator . The NBS voice time signal 
may be used for the external reference (as suggested above for 
setting the Nova's real time clock manually). After becoming 
familiar with the controls and adjustments on the time code gen- 
erator, the following procedure may be used to set it into opera- 
tion: 

Step 1. Set the POWER switch to the ON position. 

Step 2. Press the STOP pushbutton. The visual time 
display should not be updating. 
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Step 3. Set the time of day by the following procedure: 

Set the number desired in the unit seconds position on 
the PRESET thumbwheel switch; then press the SET button located 
directly under the unit seconds display; then proceed to preset the 
time desired in the tens of seconds position, pressing the tens of 
seconds SET button. Once again set the time desired for unit 
minutes in the PRESET switch, pressing the SET button under unit 
minutes. Follow this procedure proceeding from the least signifi- 
cant digit of the days display. Note that it is generally good 
procedure to preset a time as much as a minute ahead of the actual 
time. This is to give the operator time to set up the controls 
and be ready to start the generator on time. 

Step 4. Observe the external time reference, and at the 
instant it coincides with the preset time in the time code genera- 
tor display, press the START button. The time code generator 
should now start updating at one pulse per second. 

ENTERING ORBITAL ELEMENTS 


To track Landsat, the computer must predict when the satel- 
lite will rise over the horizon and compute the correct azimuth 
and elevation angles for the tracking pedestal. To predict those 
times and angles, a numerical description of the Landsat's orbit 
is entered into the computer through either the teletype or CRT 
terminal. This orbital information is contained in the element set 
provided by the North American Air Defense Command, Ent AFB, 
Colorado*. Eight of the elements (shown in figure A-ll) are im- 
portant to the NED system and their meanings, formats, and units 
are as follows: 

1. EPOCH - An arbitrarily chosen recent instant ex- 
pressed as a Julian date, at which the rest of the element set 
was determined. XXX.XXXXXXXX (DAYS) 

2. NDOT0** - First derivative of mean motion, + or - 

.XXXXXXXX 

3. 10 - Inclination. XX.XXXX (DEGREES) 


*Questions about NORADC elements can be addressed to: 

SPACE DEFENSE CENTER 
(Cheyenne Mountain Office) 

ENT Air Force Base, Colorado 

**"0" stands for zero; "0" is the 15th character of the alphabet. 


A- 20 




FT" 


ri 


T 


ARMY ENGRS WAL 


GRIFFISS ROME 
5-9 

710-324-6949 VIA 3l5-337“6275 MSG NBR 050857 
R 050857Z DEC 75 

FM SPACE DEFENSE CENTER ENT AFB COLO 
TO USA ENDE WALTHAM MA 
BT 

UNCLAS SDC-0 F05085I 0819 DEC 75 

NEDED-W/ATTN COOPER 


t. EPOCH 


2. N DOT 0 




1 0761 5U 75 4 A 75 338.62508676 -.00000340 00000-0 -23329‘3 0 01 14 

2 07615 99.0522 35. 15 44 001 3684 158. 7570 201 .41 14 13.95215112 4404 



4. NODE0 


7. M0 


5. E0 


6. OMEG0 




Figure A- 1 1 . Format of Orbital Elements provided by NORADC. 
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4. NODK0 - Right Ascension of the Ascending Node. 
XX.XXXX (DEGREES) 


5. E0 - Eccentricity. . XXXXXXX (NO UNITS) 

Notice that the decimal point is not printed on NORADC 
message, but must be supplied to Nova system. 

6. OMEG0 - Argument of Perigee. XXX.XXXX (DEGREES) 

7. M0 - Mean Anomaly. XXX.XXXX (DEGREES) 

8. N0 - Mean Motion. XX.XXXXXXXX (REVS/DA\J 

The orbital element set is entered into the system by execut- 
ing a program called "REWRITE", which stands for "Element Writer". 
REWRITE is an interactive program which guides the operator in 
entering the numbers correctly. Because the numbers have many 
digits, it is easy to mistype them on the teletype keyboard; there- 
fore, REWRITE echoes each number as it is entered and allows re- 
vision of that one number. If no correction is needed, the opera- 
tor types Y_ after "OK?" and enters the next number. If a correc- 
tion is needed, the operator types N_ and retypes the same number. 

An example of the operator /computer dialogue for the element set 
of figure A-ll is given in figure A-12 . 

DCS RECEIVE SITE DECODING EQUIPMENT 

Background . DCS Receive Site Decoding Equipment (or "decoder") 
consists of one rack-mounted cabinet containing RF and logic 
modules that demodulate, decode, and buffer incoming messages for 
presentation to a Data Genei~al 1220 minicomputer. The decoder was 
built by Robert Snyder of NASA, Wallops Island, Virginia. His 
documentation is available upon request. The functions performed 
by the decoder are shown in figure A-14, and the relationship of 
the decoder to the rest of the system is shown in figure 1. 

In the field, data collection platforms collect, encode, and 
transmit sensor data to Land sat in a manner shown In the following 
block diagram. 
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ElUPITE 

ENTER EPOCH, NDOTO, IO,NODE0,E0,OnEG0,n0,N0 
DO YOU NEED FURTHER EXPLANATION? YES OR NO N 
ANSWER YES ( V ) OR NO(N) TO OK? 

192.64706325 


192.647063250 

OK?Y 

0.0 

0.000000000 

OK?Y 

98.9660 

98.966000000 

OK?Y 

246.7513 

246.751300000 

OK?Y 

.0013733 

0.001373300 

OK?Y 

8.7956 

8.795600000 

OK?Y 

351.3396 

351.339600000 

OK?Y 

13.95255730 

13.952557300 

OK?Y 


R 

TYPE ELEMENTS 


0 . 1926470632500000E 3 

0 . 0000000000000000E 0 

0 . 9896600000000000E 2 

0.2467513000000000E 3 

0 . 1373300000000000E >2 
0 . 8795600000000000E 1 

0.3513396000000000E 3 

0 . 1395255730000000E 2 

77 
7 

12 

R 

Figure A-12. Dialogue between Operetor and Computer 
during Execution of 'ELURITE*. 
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LOCAL USER TERMINAL FUNCTIONAL DIAGRAM 
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Before each transmission, each DCP message is encoded using 
a rate one half, length five, convolutional code, to produce a 
190-bit message output (NASA, 1976a) . The message is double the 
format shown in the following table. 


Bits 



11 5 

Preamble 


16-17 

Synchronization 


18-27 

Platform ID ( 


28-35 

Data Word Nnmbor 1 j 


3643 

Data Word Number 2 1 


44-51 

Data Word Number 3 

Sensor 

52-56 

Data Word Number 4 \ 

1 through 8 

60-67 

Data Word Number 5 / 

occupy 

63-75 

Data Word Number 6 \ 

64 bits 

76-03 

Data Word Number 7 ] 


84-91 

Data Word Number 8 / 


92-9S 

Encoder run-out bits 



The Landsat spacecraft acts as a simple relay receiving, 
frequency translating, and retransmitting the message butsts. No 
onboard recording, processing, or decoding of the data is performed. 
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The output of the DCS receiver is applied to the premodulation 
processor where the DCS data are put on a subcarrier of the S-band 
equipment. This equipment retransmits the DCP messages to receive 
sites on the ground. 

At the NED Ground Receive Site, the composite S-band signal 
is received, and the DCS subcarrier is extracted and Inputted to 
the decoder which is the subject of this section. The decoder 
receives a 1.024 MHz subcarrier from the S-band receiver and 
applies the necessary operations to recover and decode the specific 
DCP data. The extracted data are formatted and inputted to the 
computer system by means of digital I/O interface and the computer’s 
interrupt system. 

In addition to the data decoding equipment described above, 
there is one other device in the decoder cabinet which performs a 
separate function. It is a special computer /tracking pedestal in- 
terface which reformats azimuth and elevation command words 
from the computer for presentation to the satellite tracking 
command equipment. 

The relationship of the decoder to the rest of the system is 
shown on figure A-14. Also shown on this figure is the data 
stream along the pathway from the receiver to the functions "ESK 
demod", "Bit sync", "Convolutional Decoder", and "Computer Inter- 
face". Greater detail of this pathway and the interactions to 
carry out these functions is shown in figure A-15. Schematic 
diagrams showing individual components are available from NED. 
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DCS DECODER FUNCTIONAL DIAGRAM 
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r cpruc 

C RNC 19 JULY 197? 


c 

C ROUTINE TO RETRIEUE DCS DATA FROM OCTAL 
C STORAGE FILE BV STATION OR DCP NUMBER 
C 

C INPUTS OCTAL DATA FILE 

C STATION OR DCP NUMBER 

C STARTING DATE 

C ENDING DATE 

C 

C OUTPUT! OCTAL DATA FILE CONTAINING 

C RECORDS FOR STATION BETWEEN 

C SPECIFIED DATES 

C 

C SUBROUTINES* OPNFLS.UDII.UICHDCP 
C 

C FILES* 'STADCP' 

C 

DIMENSION IACI€>,IP(4),IE<4).IFILE(5) 

CALL NEUPAG 

CALL 0PNFIS(2,41,3,-1 ) 

If URITECID.il ) 

11 FORMAT C 1 X, * STARTING DATE* ’,2) 

READC 11, 12 )IP(2),IP(3),IP( 1 ) 

IPY-IPC1) ,2 DIGIT VEAR 

IPC 1 )-NOD(IP(l ), 10) »1 DIGIT VEAR 

IP(4)-D {DUMMY HOUR 

12 FORMAT! 12, IX, 12. IX, 12) 

URITEC1D.13) 

13 FORMATC IX, ‘ENDING DATE* *,2) 

READC 1 1, 12 )IE<2), IEC3 ), IEC 1 ) 

IE(4)-0 

IEC3 )«IE(3H1 | MAKE INCLUSIVE 

IEC 1 )-MOD( IEC1 ), ID) iSANE AS ABOUE 

ACCEPT ‘STATIONCD FOR DCP)*‘,NSTA 
IFCNSTA.EO.D )GO TO 2D 




CALL UDII(IP(2),IP(3),IPY,NDAV) {DAYS SINCE 1 JAN 
CALL UICHDCP(NSTA,NDAY, IPD, JPRN) 

GO TO 25 

29 ACCEPT *PID» *, IPD 

25 KK*D 

KS-® 

DO 8D 1-1,3 

KC-1 {SEEK COUNTER 

30 K-KCXC 10*X(3-I ) )+KS {SEEK BLOCKS 

C TYPE ‘K-‘,K 

CALL FSEEKC2.K) 

3S READ C 2 , 4D, END -7D ) I A 

49 F0RMAT(11,SI2,A1. 14,813) 

DO 45 J-1,3 

IFC IAC J ) .EQ.9 )GO TO 3S 
45 CONTINUE 

DO 50 L-1,4 

C TYPE 1 IA, IPs * ,IACL),IPCL ) 

C PAST STARTING DATE ? 

IFCIACL).EO.IPCL))CO TO S0 
IFCIA(L).LT.IP(l))GO TO 60 



I ON THE NOSE 
I NO 


> YES 


GO TO 70 
50 CONTINUE 

GO TC 90 
60 KC-KC+l 

GO TC 30 

70 KK-(KC-1)*C 10*1 (3-1)) 

ICS-KS+KK 
80 CONTINUE 

C OUTPUT DATA 

90 IF(IA(8).EQ.IPD )URITE( 10, 196)IA 

IF(2A(8).EO.IPD)URITEC3,1D0)IA 
100 FORMATC1X, II, SI2.A1. 14.813) 

READ (2,40, END-208 ) I A 
C CHECK FOR ENDING DATE 
ICT-0 

DO 119 LL-1,3 

IFCIA(LL).GE.IECLL) )ICT-ICT+1 
IFCICT .EG.3JGO TO 200 
110 CONTINUE 

GO TO 90 

C EXIT OR CONTINUE 
200 URITEC 10,202) 

202 FORMAT (IX. ‘MORE? ‘.2) 

CALL NOECHCl.IAMS) 
IFCIAMS.E0.«9)TVPE ‘YES* 

IFC IANS. EG. 19 >00 TO 10 
TYPE ‘NO* 

CALL EXIT 
END 


jCO TO OUTPUT 
{ANOTHER BLOCK 

{BACK UP A BLOCK 
{DECREASE BLOCK SIZE 


TYPE OPNFLS 


SUBROUTINE TO OPEN INPUT AND OUTPUT FILES 

ICH-0. ..NO INPUT FILE 

IBK>0. FOPEH INPUT FILE 

JCH-0 NO OUTPUT FILE 

JBK>0. ..... ..FOPEN OUTPUT FILE 

JBK<0 .APPEND OUTPUT FILE 

SUBROUTINE OPNFLSC ICH. IBK.JCH, JBK ) 
DIMENSION IFILEC5 ), JFlLECS 
IF(ICH.EO.0)GO TO 25 
URITEC 10, 10) 

FORMATC IX. ‘INPUT FILENAME* *.2) 

READ(11,20)IFILE 

FORMAT (»A2) 

IF(IBK.EO.0)CALL OPENCICH.IFILE.l.IER) 
IF(IBK.GT.0)CALL FOPEMCICH.IFILE.IBK) 

IFC JCH.EO. 0)00 TO 40 
URITEC10.30) 

FORMATC IX. ‘OUTPUT FILENAME* *.Z) 

READC 11,20 ) JFILE 

IFC JBK.EQ.t )CALL OPENC JCH. JFILE.3.IER) 
IFCJBK.GT.0)CALL FOPENCJCH. JFILE, JBK) 

IF ( JBK . LT . 0 >CALL APPEND ( JCH, JFIli. 3, IER ) 

RETURN 

END 


52 

Si 

Cl s? 
£ !> 


T'/PE TRACKLOAD 

RLDR 4/K 16/C TRACK T1 T2 T3 UG CHIC ANT ITTO/L * 

TASKCALL fnt.lb fort. lb rrent 

R 

TYPE TRACK 

C TRACK CALLED IV UHOSNEXT OR EXEC'D BY ITSELF 
C WAITS FOR SATELLITE RISE TIRE - WIN. 

C ORIENTS ANTENNA, TRACKS SATELLITE, LOGS DATA C 

C ACCEPTS CORRECTIONS, DUMPS DATA AT END OF C 

C PASS, STOWS ANTENNA IN UPRIGHT POSITION 

C TURNS ON AND OFF CIRCUITRY IN COMMAND SI 

C EQUIPMENT, CHAINS TO 0D3 

C TDI 14 NOV 75 C 

COMPILER DOUBLE PRECISION 

DIMENSION IDO), IT (3), IT2C11 ),IT3dl ),ITSdt ) 51 

COMNON/KBLK/IT1 (11) 

EXTERNAL T1,T2»T3 
C0MM0H/KEY/KEY1,K£Y2,KEY3,KEY4 

OUR 

C KEY- '5 ARE USED TO PASS MESSAGES AS FOtlOUSi 
C KEY1 T1 TO LSI 

C KEY2 T1 TO UG 

C UG BACK TO T1 

C KEYS T5 TO Tl 

C KEY4 Tl TO UG TIRE COUNTER 

CONMON/IBLK/DAYSINMOCIR) 

DATA DAYSINM0/3l..8*..31.,3l..31.,3t.,31.,31.,3S.,31., 
•31. ,31./ M 

C TYPE 'HIT CTRL F* 

WIN- 1. 25/1441. 

C CALL DFILWC'SDF'.IER) 

C CALL CFILU( a SDF a ,8»IER) 

CALL OPEN (5, 'BCDAZEL* , 1, IER ) 

CALL APPEN0C15, a TL0G a ,3,IER> 

CALL OPEM( 12, 'STARTANGLE', 1,IER ) 

CALL 0PEN(7,'NANCLES',1,IER) 

CUROTE-I. 

CURTIRE-I. 

KEV2-I 

KEY3-I 

KEY4-I 

CALL TIMECIT.IER) 

IF ( IER. HE. 1 )TYPE ‘TIRERR'.IER 

CURTIHE-DFLOATCITd ) WDFL0ATCITC8) >/SR. )+(DFL0AT(IT(3) )/ 
I36M.) 

CALL DATE! ID, IER) 

IFCIER.NE.t )TYPE *DE',IER 
II-IDC 1 )-l 
DO 6 J-1,11 

6 CURDTE-CURDTE+DAYSINMOU) 

CURDTE-CURDTE+DFLOATC ID (2 ) ) 


XFCCJRDTE.LE. ("‘SlNCE-ZNIN ) )GO TO 51 jTASK SCHEDULING 
TYPE 'xxxxxxxx* TOO LA T E xxxxxxxxx' .TOO LATE FOR CU 
RRENT PASS 

CALL CHAINt ‘PREDICT. SUMER) 

IF ( IER.NE. 1 )TVPE *PREDICTRD‘ , IER 
CALL EXIT 

C TASKER Tl SENDS ANGLES, T2 ORIENTS, T3 GATHERS DATA 

C T5 ACCEPTS CORRECTIONS 


READ(?.51)IT1(2) 


TINES TO EXECUTE 


NOTE THAT IT1(2) IS THE LOC THAT IS NODIFIED BY TS 


FORMAT (IB) 

XTINE-(TSINCE-IDINTCTSINCE) )I24. 

IT1(4)-XTIHE (STARTING HOUR 

IT1 ( 5 ) • < XTIHE-DFLOAT ( IT1 ( 4 ) ) )S36M . 


j SECOND WITHIN H 


IF(<ID<3>/4t4.EQ.ID<3)).AHD.CID(D.GT.8T>CURDTE-CURDTE+l. CALL FOTASK(DUN,TS.IT5,IER,-l ) 
CURDTE-CURDTE+CURTINE/24. (CURRENT JUL TINE SINCE 1 RAN IF ( IER.NE. 1 )TY PE 'FOB*, IER 

READ (5, ll.END-lM )TSINCE,C< IONE-1 

F0RNAT(F13.9,F6. 1 ) CALL RECCKEYl.IONE) ^ ^ 

IF(CURDTE.GT.364. )TYPE*RESET EP.YR.tFEB IN CLOCK, WTI I, GTCALC CALL FDELY(12i) jUAIT FOR ANT TO RCH STOW POS. 
SEE P.18 IN ESCOBAL' 11*1 


ITl(8)-3 

IT1(7)-DT 

ITKin-lH 

CALL FOTASK(Di*.Tl,ITl,IER,-l) 

IFCIER.NE.DTYPE *F01 , ,IER 

ININ-IT1(S)/8R 

ISEC-IT1(S)-IRINBSR4.1 

URITE(18,R8 )IT1(4 ),ININ, ISEC 

URITE(15,R8)1T1(4),IRIN.I$CC 

FORMATC* NEXT PASS AT ■«.•» M2, * i a ,I2> 

IT2(2)-1 (ORIENT ANTENNA ONCE 

XTINE-CC TSINCC-WIN )-IDSNT ( TSINCE-WIN ) >814 . 

IT2< 4 )• IDINT ( XTINE ) 

IT2(S )■ (XTIRE-DFLOAT( IT2C4 ) ) >83SM. 

IT2<€)-3 

ITf(7)-S 

IT2(li)-2ll 

CALL F0TASK(DIR5,T2,IT2,IER,-1 ) 
IFCIER.NE.DTYPE 'FG2MER 
IT3<8)-1 
IT3(4)-IT1(4) 

IT3(5)-IT1(5) 

IT3(6)*t 
IT3(7)-S 
IT3C11 )«3SS 

CALL FQTASKCDUJ.T3.IT3.IER,-D 
IFCIER.NE.DTYPE a FQ3 a ,IER 
IT5(2)-l 
IT*(4)-IT1(4) 

ITSC5)-IT1CS> 

ITS (6)- II 
IT5C7)-I 
ITSdl )-5H 

CALL FOTASKt DUN.TS, IT5, IER,-1 ) 
IFCIER.NE.DTYPE a FQS a ,IER 
IONE-1 

CALL REC(KEY1,10NE) _ ^ 


l* 


f 




CALL ANTCH.tI.ri.il. II) 
CALL CLOSE ( 5, ItR ) 


,TURN OFF CIRCUITRY 


3> 

I 

s> 


CALL CLOSE!?, IER) 

CALL CLOSE! 12, IER) 

URITE! lt,?t) 

URITE<15.?t) 

70 FORMAT!/) 

CALL CHAIN! *004. SUMER) 

IF ( XER.NE. X )TYPE 'CHERMER 
100 CALL EXIT 

CALL CLOSE! 15, XER) 

END 
R 

TYPE Ti 

C TDI 9 DEC 75 

C SENDS ANGLE INFORMATION TO 4M5 DIGITAL I/O 
C WARD DURING SATELLITE PASS 
TASK TI 

COMMOH/KEY/KEYl , KEY2 , KEY3, KE Y4 
CORnOH/KBLK/ITKll) 

IM 

KEY4-KEY4+1 
CALL CHK!KEY3) 

IF(KEY3.GE.t)G0 TO IS 
KEY3-KEY3+1 
IT1(2)-IT1<8H1 
GO TO 25 

1A READ BINARV!S)K1 

IF<K1.EQ.-1>00 TO 3t j£NO OF FILE? 

READ BINARV(S)K2.K3.K4 jNO, GET 3 MORE NOS. 
CALL AHT(K1,K2.K3»K4,I§) iSCND THEM TO ANT 
IF!KEY3.LE.t)G0 TO 25 {ADJUST IF NECESSARY 
KEY3-KEY3-1 
IT1(2)-ITK2)-1 
GO TO It 
CALL KILL 

CALL XMT(KEV2,l,S5t) 

IONE-t 

CALL REC(KEY2,I0NE) jUAIT FOR UG TO FINISH 
CALL XMT(KEYI,I,S5t) {TELL LSI EOF HAS BEEN 
CALL KILL 
TYPE ‘XHTERRI* 

CALL EXIT 
It-1 

TYPE * ERRTN* 

CALL ANT ( It, It, It, It, It) 
URITE(l5,7t)KI,K2JC3,K4 
FORMAT! IX, 'ERRTN FROM RD BINARY' ,40I«) 

CALL EXIT 
END 
A 

TYPE T2 
C TDI 8/7/7S 

C SENDS ANGLE INFORMATION TO 4tfi5 I/O BOARD 
C DURING ORIENTATION BEFORE THE SATELLITE PASS 
TASK T2 
DIMENSION JC4) 

CALL APPEND ( IS , ■ TLOG 1 , 3, IER ) 


C 

It 


IS 

2t 


ltt 

iet 

21 


25 

3t 


13t 


25 

3t 


St 

6t 


7t 


jTELL UG TO DUMP TO DISC 


ia-e 

READ BINARY! 12 ) J ;GE T FIRST ANGLE PAIR 

FROM -startangle* 

IP0INT»1 

GO TO ltt i GO SEND PAIR 

CALL FDELY!2t) »UAIT FOR ANTENNA TO GET THERE 

DO 2t 1-1,21 tSEND ANTENNA CCU TWICE 

READ BINARY!12)J 

IPOINT-2 

GO TO ltt 

CALL FDELY!21) 

CONTINUE 
GO TO 21 

CALI ANT! J(1 ),J!2),J(3),J!4),It) 

URITE ! 15, 12t) 

WRITE! It, 12t) 

FORMAT! • .*,2) 

GO TO (It, 15,25.35 ),IPOIMT 

UUSWiiitf" " KTO "‘ “ °" e£ 

IPOINT-3 
GO TO ltt 
CALL FDCLY(21) 

CONTINUE 
READ BINARV(5)J 


IPOINT-4 
00 TO ltt 
URITE! It, 13t.) 
URITE! 15, 13t) 
FORMAT!* !*) 

CALL KILL 

CALL CLOSE! 15. IER) 

END 


REACHED 

TYPE T3 

TASK T3 

CALLS UG !UORDG£TTER) 
CALL UG 
CALL KILL 
END 

R 


{FIRST LOOK ANGLE, 
THEM CONTROL GOES TO 


\ rl 


u 





> 

I 

u> 

u> 


O O 

§| 

tr* 

£> ^ 
fg 


TYP£ QD4 

CnnttttnnntnnntiM it t i 1 * 1 1 1 1 i 1 1 t * i t : n 1 1 1 1 
C 


TDI 13 APPIL 1977 


PPOGPPfl TO CONVERT LANDSAT RSOS TO FORMAT INCLUDING 
TINE. DCP NUMBER, AND t OCTAL TRIPLETS. 

THIS PROVIDES INPUT TO B P B OR B P3 B 

PROGRAN ACCEPTS NSCS FROfl AU DCP'S AND SCREENS OUT 

THOSE WHICH ARE NOT NED'S 


C 
C 
C 
C 
C 
C 

C 
C 
C 
C 
C 

Ct > 1 1 1 1 1 1 1 1 1 u t i n t i 1 1 i < 1 1 i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

DIMENSION I A ( € ) , XI ! 12 ) , I DATE ( 3 ) , IA2( 16 ) » 1D2 (16 ) 
DINENSION ID0NE(ttl99) 

INTEGER DAVSINRO 
COHRON/ 1 BLK/DAY5XNR0 ( 12 ) 

DATA DAVSINRO/31,26,32,36.31, 30,31. 31, 30,32,30,31/ 

CALL APPEND! 15, ‘TLOC‘,3, IER ) 
ifcier.ne.dtyae •qlermer 
CALL OPEH<S t ‘SDF‘,i.IER) 

CALL DFILU(*LS2DAT*,XER) 

CALL CFILW(‘LS*DAT‘.2.IER> 

CALL 0PEN<6,*LS2DAT*,3,IER) 

CALL APPEND <7. ■rrORAQE\3, IER ) 

CALL APPEMXS, *EST0RAQE*,3,IER> 


IHR-IhR+1 

IF UHR.L". 24)00 TO S 

IRET-2 

00 TO 50 

CONTINUE 

if(ia( 1 ).eq.o)co to ioe 

I COUNT -7 COUNT +1 

DO 20 1-3,6 

ID-IAC I ).AND.37?K 

IC-ISHFT(IA(I ),-•) 

IC-IC.AND.377K 

J-IM-l 

II(J)-IC 

J-J+l 

I1(J)*ID 

CONTINUE 


10 


■OPENERR 1 


IFUER.NE.DTYPE 
WRITE < 10, 10) 

WRITE (1S, 10) 

FORMAT!//) 

NEDCOUNT-O 
I COUNT -0 
ICT-0 

CALL DATE ( I DATE, IER ) 

IYR-IDATE(3)-IDATE(3>/10*10 
IF(IDATE(3 )/4S4 .EG. IDATEC3 ) . AND. IN0NTH.E0.2) 

•TYPE ‘WATCH LEAP DAY!* 

GET STARTING TIRE 

wec- isec+so*i Ain*" 5 ironth, iday , inr, ihin, isec 

G® CONDO COUNTER 

ISC- ISC+ISEC-ICTS 3600 

IFUSC.LT. 3000 >00 TO 4 

ICT-ICT+1 

ISC- ISC-ICTS3000 

IHR-IHR+1 

IFUHR.LT. 24)00 TO 4 
IRET-1 
GO TO SO 
IRIN-USC/60) 

ISC-ISC-SOIIRIN 
IF UR IN. LT. 60)00 TO S 
IRIN-0 


ICHK-IA(2).AND.« 

IA(2)-IA(2).AND.7777K j STRIP EUERTHING 9UT THE DCP 0 

XF(1A(2).E0.S063K.0R.IA(2).E0. 6604000 TO 2S 
IF(IA(2).LT. 7000000 TO 1 
IF UA(2 ) .EQ. 70C7K )00 TO 1 
IF(1A(2).EQ. 7514000 TO 1 
IFUA(2).E0. 7346000 TO 1 

25 IER-S224 j ASSURE NO ERROR, OUTPUT A BLANK 

IFUCHK.QT.l ) IER- 17006 jIF ERR IS FLAGGED OUTPUT AN 
•E* 

NEDCOUNT-NEDCOUNT+1 

IFUH0NTH.E0.DIYR-7 

WRITE (6, 36 ) I YR. IRONTH, IDAY, I HR, IHIN, ISC, IER, IA(2). U» 
(K >,K-6, 12) 

X WRITE! 10, 35 )IYR, IRONTH, IDAY, I HR, IRXN, ISC, IER, IA(2 > . < I 

B(K),K-5,12> 

WRITE! 15,35 )IYR, IRONTH, IDAY, I HR, ININ, ISC, IER, XR!2 ), ! I 
0(0,05.11) 

CONTINUE 

FORMAT! IX, II, 512, Al, 014,1013) 

00 TO 1 


36 

50 


65 


IHR- IHR-24 
IDAY-IDAY*1 

IF UOAY.LT. DAYS IHRO (IRONTH ) )00 TO 55 
IDAY-1 

IRONTH-IRONTH+1 
IFURONTH.LT. 13)00 TO 55 
IRONTH- 1 
IYR-IYRU 

IF!IYR.0E.10)IYR-IYR-10 
00 TO (4,6), IRCT 


CiMtltttlltttlltllllttlllttttttttttttttlttltlttt 

c 


GO *0 1 


> 

1 

OJ 

• C - 


C ROUT I ME TO COMPRESS IDENTICAL LANDSAT MSGS 
C INTO ONE RECORD IN •STORAGE' WITH THE NUMIER 
C Of MESSAGES IN COLUMN 13. 

C 

C 

C i : > : : > t : 1 1 1 i i > 1 1 1 1 1 1 1 1 1 1 > > i « i 1 1 1 1 n 1 1 n t 

l®6 DO IRS I ••,193 

IRS I DONE ( I )-6 

CALL CLOSECS.IER) 

CALL CL0SE(6.IER> 

CALL F0PENC6.'L$3DAT-,41) 

IREC-R 

ISREC-® 

111 CALL FSEEK(S.IREC) 

113 READ ( 6, U5,END*3RR )IA2 

US FORMAT! Il.SI3.Al. 014. 8013) 

IREC-IRECU 

IF! IDONE! ISREC ).EQ. 1 )G0 TO 111 
ISREC a IREC 

IF! IA3! 7). NE. 17898)00 TO 181 ilTSM-'E* 
URITE(8,3S)IA8 j STASH ERRONEOUS MSGS IN 
IRR-IREC-1 
IDONECIRR)*! 

CO TO 113 
lit ISREC-IREC 

CO TO 113 

131 IA3(7)-l 

135 REA0(SU1S,END-8SR)II3 
IREC-IREC+1 
IRR-IREC-1 

IF! IDONECIRR ). EC. 1 ICO TO 18$ 

IFUI3(7).NE. 1789000 TO 136 

WRITE! 1.35)113 

IRR-IREC-1 

IDONE! IRR )>1 

CO TO 135 

136 IF! 1*3(8 ).NE.IA3!8))C0 TO 135 

DO 135 1-9.16 

IF! 1*3(1 ) .NE. IA3! I ) )G0 TO 135 
135 CONTINUE 

CO TO 14® 

14® IA3C7)-IA3(7)+1 

IRR-IREC-1 
IDONE! IRR)- 1 
GO TO 135 


'ESTORAOE* 


35® URITE(7,3S1)IA3 jPUT COMPRESSED NSCS IN 

351 FORMAT (IX, 11,513. 11.014,8013) 

IREC- ISREC 
REWIND 6 


■STORAGE* 


36® WRITE! 16,136 JICOUNT.NEDCOUNT 

URITE! 15, 136 )ICOUNT,NEDCOUNT 

136 FORMAT! • TOTAL NUMIER OF MESSAGES • \I3/* NUMIER OF 
NED 

* MESSAGES • *,I3) 

X WRITE! 10. 110 ! 

URITEC1S.116) 

116 FORMAT!/////) 

156 I F C NEDCOUNT . EQ , 6 ) GO TO 166 

CALL CHAIN! -P3.SU-,IER) 

IF(IER.NE.1)TVPE -QDER'.IER 
CALL EXIT 

16® TYPE -WAITING 1 NIN. FOR SATELLITE TO SET- 

CALL FD£LV(69® ) i WAIT FOR SATELLITE TO SET 
CALL CHAIN! -PREDICT. SU-.IER) 

END 

R 



















TYPE P3L0AD 

PLDR P3 ERTDA BlNAl BINBi DCBRY UTII UTQLY * 
FRAN LSTGE LAI SHOP 1TTO/L FORT. LB 
P 

TYPE P3 

C RMC 24 BAY 78 

C PROGRAM TO DECODE LANDSAT MESSAGES AMD PRINT 
C THEM '2UT WITH LADELS AMD ENGINEERING UNITS 


DIMENSION IDXd4,S9 ), ITEMP1 (G ), ITEMP2( 13 ) 

COMMON/ JBUC/NUM ( 24 > , SUM ( l 2 ) . X ( 8 ) , JB IN ( 9 ) . LI I N ( 3 ) 
COMMON/ JBUC/ J» K, S,F,JX, IX. NAME (13), CHK 
COMMON/ J1U/COHD, DOX, TEMP, PH 
COMMON/ KBLK/LS3< 59 ),IQD 
CALI DFILUCDATA'.IER) 

CALL CFILUCDATA' 3.IER) 

CALL APPEND<5, a DATA a ,3,IER) 

CALL APPEND CS. a TLOG a ,3,IER) 

CALL OPENC7, • INDX* ,1, ICR ) 

CALL OPENCt, a DIRC a »l,IER> 

CALL 0PEN(9, a L$89AT a ,l,IER) 

CALL 0PCN(8, a Q0IRC a .l.IER) 

CALL F0PEN(13, a LS8DIS*.408> 

IF(IER.ME.l)TYPE a OC a * ,IEM 
URITECS.81) 

UNITE (li, 81) 

FORMAT ( * PLATFORM LOCATION*, /• I.B. »•/) 

IX«N 

READ ID'S OP GAGING STATIONS 
READ (2, 17)100 
FORMAT (18) 

READ(8,lt)(LS3(IJ),IJ-l.X0D) 

READ(7,17)IDC 

READ(?,1C)( (IDX(IK,IL),IK a l,14)IL s l,IDC) 

FORMAT! 14, 18, IX, 18A8 ) 

IFLC-S 

LF6-S 

READ<i,lS,EMD-lSS)IPD {READ ID FILE 
FORMAT (14) 

IFCIPD.EQ.7187.OR.IPD.EO.784C.OR.IPD.EO.734S.OR.IPD. 
•7854 ) GO TO 85 
GO TO SS 

URITE(S.B) {SKIP LINE BETWEEN STATES 

URITEU9.8) ) a 

CONTINUE 

FORNATC/) 

DO 39 I-l.IDC 

IFdDX(l,I).EQ.IPt)QO TO 4S 

CONTINUE 

IDK-IDCd 

IFd.EQ.IDK )G0 TO ISO 
IND*IDX(2, I ) 

DO 4S J-3,18 
K-J-8 

NAME(K)-IDX(J,I) 

URITE(f,4S)IPD,NAME 






JRITE(1B,4«)1PD,NAME 
FORMAT < IX, 14, tX, 18A8, 4X,2 ) 

READ LS2DAT LOOKING FOR IPD 
READ(9,18,END*11 )ITEMPl, IPID, ITEMP8 
FORMAT (6A2, 14. 12A2) 

SORT BY ID'S INTO 'DATA' 

IF (IPID . EO . IPD )URITE(5, 19 ) ITEMP1 . IPID, ITEMP8 
IF (IPID . EO . IPD ) IFLG- IFLG* 1 
FORMAT ( 1X,6A8, 14. 12A8) 

GO TO 6 
REWIND 9 

IF(IFLG.EQ.9)G0 TO 39 
CALL CLOSE(S.IER) 

CALL 0P£N(5,*DATA",2, IER) 

DO IS IU-1.1FLG 

READ(5,2B,EMD*145)IY,IM0,IDD,IHH,MM,ISS,IERR, IPID, MUM 

F0RNATdl.SI8.Al. 14,8411) 

IFdERR.EG. 17S08.AND.IFLG.EQ.1 )C0 TO 39 

IFdERR.EQ.17S9S)G0 TO 15 

CALL UTII ( IY, INO, IDD, IHH.MM.ISS ) 

IF(LFG.NE.S)00 TO 41 
URXTE(S.49)XM0.IDD,IV 
URXTE( 19,49 )INO, IDD.IY 
FORMAT! 1H ,I8, a / a ,I8. a /\I8) 

IFG-LFO-H 
URITE(«,47)XHH,MM 
URXTEC 19,47 )XHH, MM 
F0RMAT(SX.I8, a i a ,I3, a EST a ,Z) 

00 T0(Sl.S8,S3.S4,S5.fS.C7,S8.«9,79 )IND 

URXTECC.il )NUM 

UNITE ( 19. SI )NUM 

FORMAT (BX, "TEST SET a ,lX.840Xl ) 

00 TO IS 

IF (IPID . EO . CSS4 )CALL FRAN (NUN, S) 

IF(IPID.NE.«S#4)CALL ERTDA 
IF(CHK.EQ.l. )00 TO 18S 
IFdPID.EQ. 7801 )S"S-2. 

IF( IPID. EQ.78S4. OR. IPID. EQ.73SS)S a S-19. 

IFdPID.EQ. 7873. AMD. S.GE. SO )S a S>40. 

CALL L$TGE(IPID,9,0,MANE,0A,J) 

CSM a O/DA 

URXTE(C.S4 )S,O.C$M 
UNITE ( 19,54 )*£&, CSM 

FORMAT! 13X, a STG aa ,FS.8, a FT a ,lX,F7.9, a CFS a ,lX,F4.1, a C 
GO TO IS 

CALL ERTDA _ __ 

IFCCMK.EQ.l )G0 TO 18S 
WRITE (S.SS)S 
UNITE ( 19, SS)S 

FORMAT (8X, a PRC a *,FS.8, a IN a ) 

GO TO IS 

COMD-i 

DOX-9 

TEMP-9 

PM-# 


rli*.* ... 1 ..um-uiLi,. 


1 


IfIT': 


CALL UTOLV 

ITCCMC.EQ. I. )C0 TO 185 
WHITE ( 6. SC 1COND.DOX, TEW .AH 
WHITE ! 18,58 >COHD,DOX, TEW.PH 
F0Rf»TCax»*CD-*,F6.l.aX,*D0-*.FS.3, 
«2x.*uT-*,Fs.2.ax.*PM-*,Fs.3> 

GO TO 15 

CALL $NOP<NUN t IPID, DEPTH J 

1F(CHK.EQ. 1 )GO TO 185 

WRITE!6.S?)DEPTH 

WRITE < 10,57 ) DEPTH 

FORMAT (2X, a WES-*,F5.2, *IN* ) 

GO TO 15 
WHITE! 6. 58) 

WHITE! 19,58) 

F0HMAT!3X.*CRREL*> 

GO TO IS 
CALL ERTDA 

CALL LST9E(IPID,S,Q,RAME,DA,J) 

csm-g/oa 

CALL SWOP ! NUN, IPID, DEPTH) 

IFtCHK. Eft. 1)00 TO IIS 
WRITE! 8.59 )DCPTH. S.O.CSH 
WRITE! 18.59 >D£P1Hj,0,C9N 
FOHHATiax, *WES« a ,FS.8» , IH , ,lX, , $T6“*,FS.a, 'FT*, 
81X.F7.8, *CFSMX.FS.8.*C9W* ) 

GO TO 15 
CALL ERTDA 
C0H9-8 
DOX-8 
TEW-8 


CALL UTOLV 

IF1CHK.EG.DGO TO IIS 

WRITE 18, 58 )S, COHO, POX, TEW, PH 

WRITE! 18,88 )S, COHO, DOW, TCTr, PH 

F0RHAT13X, *STO-*,FS.2,2X, *CD-*,F6.1,8X, *D0-*,F6.3, 

oax.'WT-^.Fs.a.ax, 5) 

GO TO IS 
WRITE(6.52) 

WRITE! 18,52) 

FORMAT I 2X, *DEM0. SET* > 

GO TO 15 

CALL LAB!NUN,S,ICHK) 

IFCICHK.EQ.DQO TO 12S 
GO TO 531 
WRITEI6.136) 

WRITE! 18, 138) 

F0RMATC6X, ’INVALID* ) 

CONTINUE 
GO TO 145 
WRITE!!, 168) 

WRITE! 18, 188) 

FORMAT !2X, 'NO REPORT* ) 

REWIND 5 
GO TO 4 
CONTINUE 

CALL CL09E15, IER> 


CALL DFILUf ‘DATA* , IER ) 

CALL CFILU! ‘DATA'.a.IER) 

CALL CLOSE(E, IER) 

CALL CLOSE (7, IER) 

CALL CLOSE (8, IER) 

CALL CLOSE (9, IER) 

WRITE (16,269) 

266 FORHATC ///'//) 

CALL CHAIN! "PREDICT. SU* , IER) 

CALL EXIT 

IF ( IER.NE. 1 )TVPE*CE- *,IER 
STOP 
END 
R 

TYPE ERTDA 

SUBROUTINE ERTDA 

CONHON/ JBLK/NUN ! 24 ) , LB IN ! 3 ) , JBIN ( 9 ) , X ( 8 ) . SUN ! 1 2 ) 

COnnON/JBU/J,K,S,F, JX^IX,NAHE! 12 ) .CHIC 

CONHON/ JIU/COND.DOX, TEMP, PH 

JJ-3 

S-8 

F-t. 

Jl-1 

38 J2-J1+2 

1-3 
T-8.885 
DO 18 I-J1.J2 
IFCNUNCI ).GT.T) GO TO 48 
K-Kf3 

CALL DCBRV!LB1N,NUM!I ),JJ) 

DO 18 L-1.3 
J-L-HC 

18 JBIN! J)-LBIN!L) 

Kl-E 

DO 88 1-8.9 
IF (I -8 >31,32, 31 
33 F-18.XF 

T-8.885 
Kl-5 

31 IF! JBIN! I) >88,25,28 

25 T1-F*2.W!HC1) 

T-T+Tl 

IF!F*18.-T ) 48, 28,26 

26 S-S-Tl 

88 CONTINUE 

Jl-4 

F-.81 

IF! J2-4 >38,48,48 
RETURN 
48 CHK-1. 

RETURN 

END 


1 




TYPE tINAl 
C D C 1/1/75 

SUBROUTINE BINA1 

CONNON/ J»LK/NUN< 24 ) , SUN C 12 ) , X ( S ) , JBIN ( 9 ) . L1IN ( 3 ) 
CONNON/ JBU/J.K.S.F.JXdX.NANEUD). CMC 
CONNON/ JBUC/COND, DOX, TCTP , PH 
K-3 

DO 1* 1-1.3 
IX-3*JX-<3-I>*« 

IFCNUNCIX>-7>45.45,55 
45 IC-ICO 

CALL DCBRY(L1IN,NUN<IX),3) 

DO 1* L-1.3 
J-L-MC 

JBIN(J)-LBIML) 

IB CONTINUE 

55 RETURN 

ENO 
R 


J- JJ-I+1 

iFaBiN(jnaB.2B.i5 
15 CONTINUE 

£• DO 25 I-J.JJ 

25 LBIN( I )-• 

LBIN(J)-1 
3* CONTINUE 

5 return 

END 
R 

TYPE UTII 
C RPC 13 JAN 76 

C CONCERTS ZULU TIE? TO EASTERN STANDARD TINE CEST) 
SUBROUTINE UTIKIY.IHO.IDD.IHH.NH.ISS) 

INTEGER DAVSINHO 
C0NN0N/IBLK/DAYSINN0(12) 

DATA DAYSINN0/31,29,31,3B,31,3B,31,31,3B,31,3B,31/ 
IHH-IKH-5 (CHANGE THIS IN THE FALL AND SPRING 


BINB1 

SUBROUTINE BINB1 

CONNON/ JBUC/NUN(a4>,SUHC18),X(S).JBIN<9>,l»XN(3> 

CONNON/ JBUC/J.K.S.F, J X.IX . NAME (12). CMC 

CONNON/JBU/CONfi.DOX.TCTP.PH 

DO 45 1-2.9 

R-I-t 

IF( JB1N(I ) >40.10.40 
IF <1-0)20,30,30 
X<N)-8W(6-I) 

GO TO 45 
X<n>-2*S(»-I) 

GO TO 46 

X(N)-B 

CONTINUE 

n-o 

DO SB N-1.2 

R-2BJX-C8-N) O O. 

SUN(N)-0 
DO SB 1-1,4 
LX-4BN-(4-L) 

SUN(N)-SUN(N)+X(LX) QZ 

XCLX>-0 SjP- 

CONTINUE 5° t* 

IF (SIRICA )-10 ) 00, 70,70 © h 

*c-a c5 

CONTINUE >2: 

RETURN H £ 

END H ® 



DCBRY 
C 1/1/75 

SUBROUTINE DCBRV(LBIN,NUN.JJ> 
DINE NS I ON LBINO) 

DO 1R I-l.JJ 
LBINC I )-• 

IF(NUN)5,5,35 
DO 3B ICK-l.NUN 
DO 15 I-l.JJ 


C FOLLOWING STEP IS TYPICAL TRANSITION NEA6URE 

C IF C IRO.EQ.ll.OR. ( IDAY.EQ.3B.AND.1HR.GE.8B) 

C S.OR. IDAY.E0.31 JIHH-IHH-l (RAKING A 5-HOUR DIFFERENCE 
C TEMPORARY STEP 

IFUSS.0E.3B)WWtHl 
IF(NN.GE.S0)IHH-IHH+1 
IFCNR.QE.BOMN-Nft-OO 
IF ( IHH.LT.0)IDD-IDD-1 
IF(IHH.LT.0)IHH-IHH+84 
IFdDC.LT.l HNO-INO-l 
IFdNO.LT.l )IN0>IN04ia 
IFdDO.LT.l )IDD*DAV$INNO(INO) 

IY-IY47B 

RETURN 

END 

TYPE UTOLY 

SUBROUTINE UTOLY 

CONNON/ JIU/NUP<24 ),SUN( 12 ),X(8 ).JBIN(B).UIN(3) 
CONNON/ J1UC/J.K.S.F.JX. IX, NONE (12). CMC 
CONNON/ J1U/C0ND, DOX, TEHP , PH 
DO 25 JX-l.B 
CALL BINA1 
CALL IINB1 
86 CONTINUE 

IF(K-6)lB,3t,3B 
30 CONO-B 

DOX-B 
PH-B 

COND-t .5*( 100*SUN<4 )*10*SUN< 1 )+SUN(2) > 

DOX-B.BOt IBB*SUN( 11 H10*SUN( 12 J+SUNd ) ) 
PN-0.014t(100*SUN<5)d0*SUN(6HSUN(?>) 
TENP-0.04t(100tSUNtSH10*SUN<9HSUN< IB) ) 

RETURN 
IB CH9C-1. 

RETURN 

END 


45 


le 


88 

55 


FRAN 

SUBROUTINE FRAN (NUN, SUN) 

DIMENSION NUPC84 ), JBIN(9),LBIN(3) 

$un-8 

DO 80 JX-1.2 
K— 3 

DO 18 1*1.3 

rx*3*jx-(3-i ) 

IF(NUNCIX)-7)4S,4S.5S 

K-K+3 

CALL DCBRY(LBIN,NUM( IX ),3 ) 

DO It L- i , 3 
J-l+K 

JIINCJJ-LBINU) 

CONTINUE 
DO 88 IW- 1.3 
Jl*2 
J2-5 

IF( IW.EQ.2 )Jl-6 
IFtIU.E0.2)J3*S 
iE*e-iu-e»ux-i ) 

DO 88 IC-J1.J2 
IV-IC-2 

IF(IC.G£.6)IV-IC-€ 

U-IV 

XE-IE 

sun-sun^( j>in( ic )>(2<xu >t( lttxxc > ) 

CONTINUE 

RETURN 

END 


TYPE LSTGE 

C SUBROUTINE TO CALCULATE DISCHARGES AND VALIDITY FROM 
C VALUES OF STAGE AND STATION NUMBER 


SUBROUTINE LSTGE (ISN, ST ACE, Q, NAME, DA, J) 
DIMENSION DISCH(44),HANEd8) 
C0NN0N/KBU/LS3(58),IQD 


33 


65 

C 

78 

S3 

153 

788 

68 

C 

34 

385 

58 


DC 68 — 1.43 
S-STAGE -T-NXO 
R-CS*G)/Q 
2*DISCH(N) 

P-DISCH(N*1)-Z 

IF(P)24,60,65 

IF(S)78,?0,68 

INTERPOLATE ON RATING TABLE 
IFt 2)53,53,780 
IF(R)51, 153,153 
Q* Pi CRM 1.5) 

GO TO 58 
0*2+P*R*.885 
GO TO 58 
CONTINUE 

EXTRAPOLATE ON RATING TABLE 
Q-Z+R8(2-DISCH(N-1 ) >4.885 
J-l 

RETURN 

END 


DO 85 N-l.IQD 

IFt LS3 ( N ) . EO . ISN )G0 TO 8 

85 CONTINUE LET 

GO TO 385 

8 NN*N-1 38 

CALL FSEEK(13,NN> 18 

C T* INITIAL STAGE IN TABLE; Q-STAOE INCREMENT IN 2ND LINE 

READ (13, 3 )NAME. DA , T, 0 . DISCH 

3 F0RNAT(6X,lRA2,24X.F6.8.12X,F6.R.SX/C8X,9Fi.R)) 

C J*0 FOR VALID OR J-l FOR NON UALID STAGE 
j*e 

IF(0)58, 58,21 iFOR DUMMY FILE, INCRMN T (0) WILL BE 8 
21 IF CSTAGE-T >51.22.22 

51 0— .885 

GO TO 285 


R 

TYPE LAB 

SUBROUTIIC LAB (NUN. SUP, I CHK ) 

C ROUTINE TO DECODE DATA FROM LA6ARQE CONUERTIBLE DATA 
C COLLECTION PLATFORMS. 

C EACH OF FOUR 16-BIT WORDS IS CODED BY THE DCP AS FOLLOWS > 
C .88 .84 .82 .81 .8 .4 .2 .1 8 4 2 1 88 48 28 18 
C FROM LEFT TO RIOHT. 

C INPUT TO THIS ROUTINE IS IN NUN (34), 24 SINGLE DIGITS 

C REPRESENTING DIGITS IN OCTAL TRIPLETS 
C THE 24 DIGITS MUST 8E GROUPED HERE INTO 3'S 
C EACH TRIPLET IS EVALUATED IN OCTAL l OCT(l) AND OCT (2) 

C SIT POSITIONS OF TtCSE OCTAL WORDS ARE TESTED AND 
C EVALUATED FOR BCD SIGNIFICANCE 

INTEGER 0CT(2),HUN<24) 

ICHK-8 
DO 18 1*1,2 
OCT(I )*8 

DO 38 J-l.3 
JJ*(I-1)83+J 

OCT(I )*OCT( I HNUM( JJ)8(*88(3-J) ) 


CONTINUE 

SUM-8. 

DO 48 1-1,2 


CONTINUE 


jMAKE AH OCTAL TRIP 
;FROM 3 SINGLE DIGITS 


; OCT (11 OR OCT (2 ) 

DO 68 J-l. 2 i LEFT OR RIGHT 

DO 68 IC-1,4 ;BITS, L TO R 


WC-12-48JHC ;BIT MAPPING FROM 7 TO 8 
N-OCT(I) 

CALCULATE THE DECIMAL AND BCD VALUES 

SUM-SUM- ( ITEST ( N , KK ) >8 ( 1 8 . K ( ( (I - 1 )*2+ J >-3 ) > * 








'2tT(4-iO> 

60 CONTINUE 

5# CONTINUE 

40 CONTINUE 

IF(SUf1.LT*t,t.(».SUri.GT.99.W)ICHlC-I 

RETURN 

END 

R 

TYPE SNOP 

SUBROUTINE SNOPCNUTI, IPID, DEPTH) 

DIMENSION NUN<24) 

C0ff10N/LlLK/ID$C5> 

COmON/LlU/COEFFS ( 5 ) 

CO WON/LlUC/XOfTSTTS C 5 ) 

DATA IDS/7147, 7385, 3»f/ 

DATA C0EFFS/.45S, «§81,3S#./ 

DATA X0FFSET$/-.45,-.173,3*#./ 

DO 3i 1-1,5 

IF(IPID.EO«IDS(I ) )GO TO 35 
3# COtfTINUE 

35 CONTINUE 

DS*64*NUN( IS HI*NUW< 17 )+NUN< 18 ) 
D7*S4*NUmi9H#SNUma#HHUm81) 

IFCDS* HE. #.#)CO TO 4# ) AMO ID DXUIS20N BY t 

DEPTH*#,# 

GO TO 5# 

4# DEPTH* (D7/D#)*C0EFF$( I HXOFFSETSM ) 

5# RETURN 

END 
R 

TYPE ANT 
•TITL ANT 

I SOFTWARE BETWEEN T1 OR T2 AND 4#S5 I/O BOARD 
• ENT ANT 

.EXTD .CPYL, .FRET 
•NREL 

1—157 ;LO AZ 

J-I+t | HI AZ 

PC-I+8 ,LO EL 

L-I+3 jHIEL 

A-I+4 

FS.-5 jPROG/STDBY SWITCH 
FS. 

ANTUSR • • CPYL 
STA 3 RTN 
LDA I #1 3 

NOU • • SZR ;TURN OFF PEDESTAL? 

JNP OFF ; YES 

LDA • ATI ; NO# SEND ANGLES 

STA # TT1 

LDA 2 CN4 

NOR: LDA e II 3 

LDA 1 #TTB 

CON # # 

AND 1 # 


con # e 

ISZ TTB 


DOAS B DUC 
INC 3 3 
INC 8 2 SZR 
JNP NOR 
BIO LDA 3 RTN 
JSR i.FRET 

OFF i LDA • 1# 
DOAS # DUC 
JNP BK 


tDATA AREA— 
DUC-42 
ATI* .+1 
151 135777 
141 133777 
131 187777 
l€i 117777 
Mi 1B#+112 

RTN*,-. 

CN4* -4 




.END ANT 

R 

TYPE UG 
•TITL UG 

j ACCEPTS NESSAQES FRON CONVOLUTIONS DECODER 
{AND IS ACTIVATED IV AN EXTERNAL INTERRUPT 
.ENT UG 

•EXTN .UIEX, .RSC, ♦ IXWT. .TASK, . AKILL, .XWT 
.EXTD .CPVL, .FRET 

.CONN KEY 4 j LABELLED CONNON AS IN LS 

.TXTR 1 

*NREL 

i DEVICE CONTROL TABLE 

IDDCTiDCT 

DCTt 

117 

ISR42__ ___ 

j INTERRUPT SERVICE ROUTINE 
{RESPONDS TO ONLY ONE INTERRUPT 
{FRON 4 ACS INTFC. LOADS S WORDS IN A 
{ROW AS FAST AS DECODER PROVIDES THEN. 
{TIMING HATCHES DECODER'S EXACTLY 

I SR 42 t NIOC DVC 
STA 2 URTN2 
STA 3 URTN3 
LDA • SYNC 
JSR 2.T0BUF 


LDA 1 CNS 


original page is 

OF POOR QUALOT 


nop- oia • Due 
JS* •.TOBUF 
JSR TNR 
INC t 1 SZR 
JNP NOR 
JNP OUT 

.TOBUFi TOBUF 

TOBUF* LDA 2 f.NBP 

STA • • 2 

INC 2 2 

STA 2 •.HIP 

JHP • 3 

TNRt LDA 2 CN1S 
INC 2 2 SZR 
JNP .-I 

JNP • 3 

OUT* LDA • tBNESS {GET TINC COUNTER 
LDA 2 f.NBP 

STA • • 2 j STORE TINE UITN NSC 
INC 2 2 
STA 2 «. 


SU» 1 
LDA 2 
3 


1 

URTN8 

URTN3 


NIOS DUC 
. UIEX 


URTN2*.-. 

URTN3I .-. 

SYNC* 12214 
CNSt -S. 

CH15t -IS. 
CN2lt#l-21Bt. 

CTRl 

.NIP i NBP 

UGsJSR B.CPYL 
STA 2 AC2 
STA 3 AC3 
SUB B • {GEN A • 
2 CN21M 


LDA 

STA 2 CTR 
LDA 2 PBUF 
STA • • 2 
INC 2 2 
ISZ CTR 
JNP .-3 
LDA • DUCN 
LDA 1 IDDCT 
. SVSTN 


j INIT BUFFER TO ALL t'S 


{DEFINE 


DIGITAL I/O BOARD TO SVSTEN 


. IDEF 
JNP t.ERT 
.SVSTN 
.GDAY 
JNP t.ERT 

2 3 »NEXT 
‘ PBUf 
NIP 
• 2 
2 

• 2 

2 

NIP 


2 

2 

1 

2 

• 


7 LINES STORE VR,N0,DAV IN SDP 


7 LINES STORE HR,H1N,SEC IN SDP 


nou 

LDA 
STA 
STA 
INC 
STA 
INC 2 
STA 2 
.SVSTN 
.GTOD 
JNP t.ERT 
NOU 2 3 {NEXT 
LDA 2 f.NBP 
STA 3 • 2 
INC 2 2 
STA 1 • 2 
INC 8 2 
STA • • 2 
INC 2 2 
STA 2 f.NBP 
LDA f PRESS 
SUB 1 1 
NIOS DUC 

.REC jUAIT HERE FOR LAST LOOK ANGLE TO BE SENT 


.XRHU 
JNP t.ERT 

W iPOIf 

SUB 1 1 {DEVICE CHARS 
.SVSTN 
.OPEN S 
JNP t.ERT 
LDA t PBUF 

NOVEL t t jBVTE POINTER 
LDA 1 C42tt {BYTE COUNT 
.BYSTN 

.URS t |DUNP THE BUFFER 

JNP t.ERT 

.BYBTN 

.CLOSE f 

JNP t.ERT 

LOA 2 AC8 

LDA 3 ACS 

LDA t ANCSS 

SUBZL 1 1 

.XNT 

JNP t.ERT 
JSR t.FRCT 

{DATA AREA 

Mlt,*. 

AC3< 


meSSi.GADD: ttV.l" i POINTS TO END ELEMENT IN LMEUED COIWON, “KEY* 
iNESSt .6ADD KEY, 3 jPOINTS TO TINE COUNTER (SECS) IN TI C 


ASDF: .♦l«a 
•TXT *SDF‘ ; SATELLITE WT A OUTPUT FILE 

C42M: 42*#. 


DUC1U42 

DUC *42 C 

NtP*t ; NOUAtLE BUFFER POINTER C 

.ERTiERT C 

CRT x . SYSTW 

.ERTM C 

PBUFt-BUF 

BUF: .ILK 21M. * BUFFER FOR 3N NSGT7UDS/NSQ C 

•END UG 3i 

R 

TVPE OFF 

1*1 

CPU ANTCI,I,I,I,I> E B 

STOP 

END 

R 

TVPE UHOSfCXT 

CONPILER DOUBLE PRECISION 
COWWi/IBU/DRVSINNO(12 ) 

DINENSION IT(3),ID(3) 

DATA DftVSItM)/31» ,BS« ,31 • ,3#« ,3l«,3t.,31*,31 «,3i*,31., 
*3#. ,31./ 4# 

11*1 

xmmzN*s./(M.x€4.) 

XAFTER*23*/(M.U4. ) 


li 


IS 


CURDATE*CURDATE+CURTinE/24. 


PRESENT JUL TINE SINCE 1 JAN (DAYS) 

XINTRUL*DFLOAT( INTRUt )/24. 

STAGETINE*CURDATE-DNOD ( CURDATE, XINTRUL 5 
I HOUR* ( STACETINE-IDINT £ STACETINE ) )*24 , POTENT I AL 

i 

' STARTING HOUR FOR STAGE 

READ STARTING TINE FROM 'UANTSTAGE' 

READ ( S. 3# ) I DUN, DSTGTWE 

THAT UAS CALCULATED PREUIOUSLY AND STORED ON DISC 
FORMAT (I5/F13. 9) 

CALL CLOSEC5.IER) 

IF ( STAGETIHE . LT.DSTGTINE JSTACETINE-DSTGTINE 
IF(STAGETINE.LE.(TSINCE-XHARCIN))GO TO 2GOj XEQ 'STAC 

TFCSTAGETIHE.CE. (TSINCE+XAFTER))GO TO 3GG »XEO ’LSI* 


CALL ORCNCS, ' UANTSTAGE*, 1, IER) 

IF ( IER.NE. 1 )TYRE ‘US ERR', IER DM 

READ(5,1#)IYH jNO a G,YE$ a l 
FORMAT! II ) 

IFUYH.EO.GJCAll CHAIN! 'TRACE. SU* , IER ) 

IFCIER.NE.l )TYRE ‘LS CNN £§• IER 

READ (5, 15)IHTRUL , HOURS DETU. INTERROGATIONS 

FORMAT (IS) 3M 

CALL ORCNCS, 'BCDAZEL',1. IER) 

R£ADC6,2G)TSINC£ 

F0R(1AT(F13.9) 

CALL CLOSE(S,IER) R 

CURDATE-G. 

CURTIME-0. 

CALL TINE! IT, IER) 

IFCIER.NE.l )TYR£ •TIMOW'.IE* 

CURTIME«DFLOATCIT(l ) >>(DFLOAT(IT<a) )/#•. )+(DFL0AT(IT(3) )/ 

I3SM y 

CALL DA TEC ID, IER ) 

IFCIER.NE.l )TYRE *DE*,IER 

II-IDCD-t 

DO S J-1,11 

CUR DATE “CUR DATE +• DAYS INMO ( J ) 

CURDATE •CURDATE+DFLOAT ( ID ( 2 ) ) 

IF C ( ID (3 )/«4 . EQ . ID C 3 ) ) . AMD . < ID C 1 ) . GT . I ) )CUROATE a CUROATE*l , 


IHOUR- (STAGETINE-IDINT! 


TER i ADD 
STAGETIME 


))S24 


CALL CLOSE (5, IER) 

CALL ORENtS.tUANTSTAQEM.IER) 

URITECS, 4# )IYM, INTRUL, IHOUR.STAQETIME 
FORMAT! IX, I 1/IS/IS/F14 .9 > 

CALL CLOSE(S.IER) 

CALL CHA1N( 'TRACE. SUMER) 

TYRE 'CHAIM TROUDLE' 

CALL EXIT 

TYRE'LEAUING 'UHOSNEXT' AND CHAINING TO 'AUTOSTAGE' 

CALL 0RENC5, 'UANTSTAGE' ,3, IER ) 

URITE C S ,4i ) IYM.INTRVL , I HOUR , STAGETIME 
CALL CLOSE(S.IER) 

CALL CHAIN! ' AUTOSTAGE. SU‘,IER) 

GO TO 2SG 

TYRE 'CHAINING TO 'LSI" 

CALL CHAINC 'TRACE. SUMER) 

GO TO 25# O O 


END 


*=0 50 





tv pe setstage 

COMPILER DOUBLE PRECISION 
CONMON/IiLK/DAYSINNO! 12 ) 

DIMENSION IT(3),IDO) 

DATA DAVSINBO/3 1 • # 28 • » 31 • # 38 « # 31 « # 38 •«31*»31«#38*#31*# 

,3# a6cEPT ‘USE DEFAULT TIMES? 1-YES.O-NOMANS 
IFdANS.EQ.DCO TO 29 

ACCEPT * INTERROGATION INTERVAL (HOURS )«‘ t INTRUt 
TV PE ’HOW MANV MINUTES AFTER 1130'S INTERROGATION?’ 
ACCEPT ’(3# MINUTES IS SUITABLE X \XOFF 
XOFF*XOFF».IOO S9g jCONUERT TO DAYS 
GO TO SB 

2# INTRVL • S 

XOFF-O. 


CALL 0PEN!5,*UAHTSTAQE’«2,IER> 

IFdER.NE.DTYPE *«S ERR’.IER 
CURDATE*0. 

CURTIME*#. 

CALL TIMEdT.IER) 

IFdER.NE.DTYPE 'T1PCRRMCR 

CURTIME*DFLOAT( IT( D H(DFLOAT( IT(2 ) )/BB. H(DFLOAT( XTO ) >/ 

888 ) 

CALL DATE (ID, I ER) 

IFdER.NE.DTYPE ’DC’.IER 
11*10(1 HI 
00 < J-l.II 

CUR DATE* CUR DATE+DAYS INMO! J ) 

CtfltlMTE • CLRDATE+DF1.MT ( ID ( 2 ) ) 

IF! (10(3 )/4*4.E0dDC3) ).ANO. (ID(D.GT.t ) >CURDATE*CURDATE-H . 
CURDATE-CURDATE+CURTIHE/I4. jTIME SINCE 1 JAN 
XINTRUL*DFL0AT(INTRVL)/24* 

STAGETIHE*CURDATETXINTRVL-DN0D<CUR9ATE, XINTRVL HXOFF 
IHOUR* (STAGETIME- ID I NT ( STAGETIME ) )*94 


C ! POTENTIAL STARTING HOUR FOR STAGE 

WRITE ( 5, 30 ) I 1, INTRVL, IHOUR, STAGETIME 
WRITE! 10, 40 )STAGETIME. INTRVL 

40 FORMAT ( IX, * NEXT NETS INTERROGATION IS AT *,FU.C/ 

BIX, ’AND EVERY '.la,' HOURS THEREAFTER’/) 

30 FORMAT! 1X,I1/IS/IS/F14.9) 

CALL CL0S£(5.IER> 

TYPE ’ CHAINING TO CLOCK* 

CALL CHAIN! ’CLOCK. SUMER) 

END 

R 

TYPE CLEARSTAGE 

10*0 

xo*o. 

CALL 0PEN(5»*MANTSTAG£*,3, IER) 

WRITEC5, 10)10, 10. IO.XO 
10 FORMAT! lX.il/lS/I6/rl4. 9) 


\ li l* V ^ ^ U 


CALL CLOSE (S, IER) 

TYPE ’CHAINING TO CLOCK. SV 
CALL CHAIN! ’CLOCK. SUMER) 
STOP 
END 



A-43 


| i ' 


type clock 
j 7/1/75 TIP'S 
{REVISED 4/14/74 

{PROGRAN TO READ COORDINATED UNIVERSAL 

j TINE FRON Tire CODE GENERATOR 

{THE SIGNAL LEAD FRON THE TINE CODE GENERATOR 

jNUST K PLUGGED INTO SOCKET P7 OF THE 

{C0NNUN1 CAT IONS RULT I PLEXOR PLUG STRIP 

• TITL CL 

.ENT CL 

•EXTH .UIEX,.REC..IXNT». TASK.. NEC,. KILL 

•TXTN 1 

.NREL 


CLiJNP .+1 
IDA 4 CUD 
LDA t ACDCT 
.SYSTR 

. IDEF , DEFINE 4G2S 
JNP .ERT 


LDA 1 .S tSTRTG ADDA OF TASK 


UK I .9 tSTWTIi nil 

LDA • DPRI jPRI*# 

.TASK 

JNP .ERT 


StSUD G • 

STA G II iINIT tS FLAG 


Ct SUD G G jINIT COUNTER 

STA • CTR 

STA • NIN 

STA • m 

STA • DAY 

STA • UT1C 


AND TIRE VALUES 


Dt LDA • AHLOC )ADDR OF NSQ LOCATION 
NIOS DUC 

•PEC iUAIT FOR NSO FRON ISR 
LDA • 1IT7 

AND • 1 Sm j DATA IN AC 1, CHECK FOR HI OR LOU 
JNP LO i LOU 

LDA • CTR (HIGH, INC THE COUNTER 
INC • 0 
STA • CTR 
JNP D 

LOi LDA 1 CTR 

LDA • C44# 

ADCZt 1 G SZC jIS CTR>»44G? 

JNP CMC# jNO 

SUDZL 1^1 >YES,CHK II TO SEE IF ONE .tS ALRDV 
SUM • 1 SZR 

JNP SET jNOT SET, THIS IS THE FIRST .8S 
JNP CHK1 ; SET, THIS IS THE SECOND .«S 
CHKSt NOV I 1 SZR j ACl HAS CTR TO SEE IF 





W SIGNAL USED TO BE HIGH 

S j NOT G.USED TO »E HI 
JNP D iG. NEVER UAS HI 
SET* SUIZL 1 1 

STA 1 II {FLAG RECEIPT OF FIRST .8S 
JNP C {GO WAIT FOR ANOTHER 

.ERTt .SVSTN 
.ERTN 

jDATA AREA 

Ka !_ a 

izRCL 

. Ai A 

♦CRSTi CRST iADDR OF COUNTER RESET 
.Si S jADLR TO START OVER 

CTR I {STORES NO. OF 4G26 TICKS THAT SIGNAL IS HIGH 

Hi.-. (LOOP COUNTER FOR PROGRESSING THRU NINUTE 
UAL«.-. 

NINt.-. {CALCULATED NINUTE 

HRi.-. ,CALC HR 

DAY I {CALC DAY 

UTICi .-. j CALC UT CORRECTION 

imi 117 {DATA CONES IN ON LINE 7 OF 4G26 BOARD 

ANTDLI NTDL {POINTER TO TIL OF SUDROUTINE HANES 

A DCDlDCD {POINTER TO TDL OF DCD VALUES 

iDLAYt.-. {FUMS WHETHER TO DELAY OR ADUANCE UT1C 


III .FLAG 1ST .SS 

AHLOCt NIOC {NSQ LOCATION FOR .XNT, .REC 

NLOCt 

AC# i .-. {SAVE ACS 
ACS I .-. {SAVE ACS 

dpri t s {Task priority 

CLKDl 24 (4G2S DEVICE CODE 

DVC-24 {4«2S DEVICE CODE 

C44GI 4GS. {LOWER UNIT TO SIGNAL .SS 

ACDCT i CDCT jADDR OF DEVICE CONTROL TADLE 

CDCTi .-. 

IDS 

ISR 

ISRt DIDC 1 DUC {4DBS ISR 
STA I ACS 
STA 3 ACS 
LDA S BIT7 

SUM S 1 SNR {BIT 7 SET? SKIP IF NO 
JNP XOUT {OK AS IS 

LDA 1 KS I NO. IT'S AG— NAKE IT A 2- 
XOUT t LDA • ANLOC 

. I XNTi SIGNAL NAIN PROO THAT I NT HAS OCCURRED 
JNP .ERT 
LDA 2 ACS 
LDA 3 ACS 

SUDZL 1 1 {FORCE RESCHEDULING 


SKIP IF NO 


U 



Mini Din ,TMLE OF no. OF DAYS IN nONTH 
Dim 3i. 

n. iXtmmnnMTCH fok leap vems i mmsmi 


ailda t Min i next a* lines calc day no. in vn 
com LDA its 
now i i szn 

JNP CONI 
SUBZL 3 3 

SUt 3 I 

jnp coni 

conn now i i snr 

JNP CONS 
LDA 0 DAY 

SUB 1 • 

STA • DAV 

N0W2LS • • SNC ; IS DAV<t? 

JNP non i NO 
JNP CONI iVCS 

nomnow t t snr ,is dav down to zero? 

JNP COM | YES 
INC » I i NO 


INC I I j NO 
JNP CON 

CONSi LDA 1 • 2 
LDA • DAV 
ADD 1 • tRESTORE 
LDA 1 MIR 

SUB 1 I 

INC I 8 | NON. NO. 


.uiex 

CMKlt SUR » • 

STA • H 

CHK2t JSR •.CAST 
IDA • H 

INC • • jN POINTS TO TAKE VALUES 
STA # N 

LDA 1 C42 j AND SUBROUTINES TO CALC TINE 
SUB • 1 SNR .DONE LOOPING THRU SUBROUTINES 
JAP S.A ,YES 

R2> LDA • ANLOC {NEXT 24 LINES DET'N S,l, OR NARK PULSE 

NIOS DVC 

.REC 

LDA • BIT? 

AND • 1 SNR jDATA IN AC1 
JHP LOU1 
LDA • CTR 
INC B • 

STA • CTR 
JHP R2 

LOUlt LDA • CTR 
HOU t • SNR ,CTR»S? 

JHP R2 iVES 

LDA 1 LLS7S j NO, CMC FOR 27S 
SUBZS 1 • SZC iIS 27SJCTR? 

JHP CHK3 iNO, .5 OR .1 

JHP LDS jtES, . 2, i RCUD 

CHK3* LDA 1 UL27S 

SUm 1 • SNC jIS 2RB>CTRT 

JHP LD1 j NO, 1 RCUD 

JHP LD2 jVES, HARK PULSE RCUD 

tNOU PASS B.l, OR 2 TO SOHE SUBROUTINE 
LMt SUB • # j GEN A • 

JHP STC 

LDlt SUBZL • • jGEN A 1 
JHP STC 
LMt SUB • • 

INCZL • • ; GEN A 2 

STC« STA • TC iHAVBC NOT NEC IF ACS PRESERVED 

LDA 2 ANTBL »GET ADR OF SUBROUTINE TABLE 

LDA 1 N jLOAD NO. OF STEPS UE HAVE GONE THRU TABLE 

ADD 1 2 jCALC ADDS UITHIN TABLE 

LDA 1 • 2 jGET SBRTN NAME IN ACL 

STA 1 NAHE ; STORE IT 

LDA 2 ABCD ,GET ADDR OF BCD TABLE 

LDA 1 N 

ADD 1 2 ; DET'N ADDR UITHIN TABLE 

LDA 1 S 2 ,G£T BCD VALUE 

STA 1 UAL jFOR CURRENT SECONO UITHIN HINUTE 

JSR SNAHE | GO TO THAT SUBROUTINE 

JHP CMK2 ,60 UAIT FOR NEXT SECOND 


LDA e HIN | YES 
ADD 1 S tACCUn 
STA S HIN 


1 I YES 
t ACCUN 


jup bk 


HRAD > LDA S TC 

HOVZR S 9 SNC 

JHP BK 

LDA • HR 

ADD 1 S 

STA S HR 

JTP BK 

DAVAD i LDA # TC 
HOVZR S S SNC 
JHP BK 
LDA S DAY 
ADD 1 # 

STA S DRV 
JHP BK 

CHK.lt LDA S TC 
SUB 1 1 

INCZL 1 1 jGEN A 2 

AND S 1 SNR , IS .t ON SCHEDULE? 

JHP t.S jNO 

JHP IK jVES 

BLANK t JHP .41 
JHP S 3 

BKt JHP t 3 , GENERAL RETURN 

jDATA AREA 

UL275t3BS. , UPPER LIHIT OF .SS 
LLSTSt 2SS. , LOWER LIHIT OF .SS 
C42* 48. i HO. OF SUBROUTINES 
JHP . ,TEHP 

TCI ,TIHE CODE S.l^OR 2 
NAHE > } VARIABLE SUBROUTINE R 


, AC TION TABLE 

AtBLI S jFIRST LOCATION IS DUHHY 


HINAD t LDA • TC 
HOVZR • S SNC ;Uf 
JHP BK ;N0 


(GET TINE CODE RECEIVED 
t TC-1? 


S&SLa rj 


*oo a i 

LDA 2 YR 

.SYST* 

.SDAY 

J*P •* AERT 


J*P F jOATA 

YRt 9. i DOC SAYS 1977 IS V* 9 

.AERT: *tRT 


CRST* SUI • • i COUNTER RESETTER 

STA • CTR 

J«P • 3 iRETURN 


£ 41 * 41 . 

C1NUW. 

C1HU3S. 


CBSStSSS. 


chh*hh. 


sum.-. 


jDATA MCA 


Dicom z 


FtLM 1 DLCOM 

.SVSTH 

•DELAY 

jrp .♦! 

LOA • S41 jAT LAST! ! UE CAR SET WE RDM. TIRE CLOCK 
LDA 1 AIN 
LOA 2 HR 
.SYSTR 
.STOD 

JAP S.AERT 

LOA • BPSF 

LOA 1 OH 

.SYSTR 

.EXEC 

JRP .AERT 

OH i 1H 

OPSFI .4-1X2 

.TXT /UHOSHEXT.SU/ 

.END CL 
R 


? 1 



TYPE 


56? 

20 


39 

31 

33 

29 


49 

59 


ELURITE 

COMPILER DOUBLE PRECISION 
DIMENSION IDAVO) 

CALL 0PEN(5, ‘ELEMENTS*, 3, IER) 

IFdER.NE.l )TYPE ‘OE’.IER 

TYPE ‘ENTER EPOCH, ND0T9, 19, N0DE9,E9,0NEG9.M9»N9* 
URITEC 19,56?) 

FORMAT (• DO YOU NEED FURTHER EXPLANATION? YES OR NO 

READ (11,29) IANS 

FORNAT(Al) 

IF(IANS.EQ.228l6)C0 TO 299 

TYPE ‘ANSWER VES(V) OR NO(N) TO OK?' 

DO 59 1-1.8 
READUDX 
URITE(19.33)X 
FORMAT (F29. 9,' OR?', 2) 

READ (11,29) IANS 
FORMAT(Al) 

IF (I AMS. 1C. 29999)60 TO 49 


TYPE ‘ENTER IT AGAIN' 
GO TO 31 
WRITECSIX 
CONTINUE 

CALL DATE( IDAY.IER) 
URITE(5)IDAV(3) 
URITE(5)IDAY(1) 
URITE(5)I0AV(2) 

CALL CLOSE (S.IER) 
IFCIER.NE.l )TYPE *CE' 
CALL EXIT 


SNPHI-DSIKPHI ) 
CSPHI-DCOSCPHI 3 
C0RD(4>SNPHI 
C0RD(5 3-CSPHI 
C0RD(1 3-369. -UL jEAST 
SNS-SNPHIXX2 
F- 1/298. 3 

*,Z) . ES-2XF-FXF 

G1-1/DSQRT( 1-ESXSNS) 
G2- ( ( 1-F )XX2 l/DSORT ( t- 
CORD (2)— GltCSPHl 
CORDd)— G2XCSPHI 
CORD (6)— G28SNPHI 
DO 69 J-1,6 

69 URITEtf )CORD(J) 

19 CONTINUE 

CALL CL0SE(8,IER) 

END 


,IER 


TYPE ‘EPOCH - T9, JULIAN DATE OF EPOCH (DAYS) F99.8<15>‘ 

TYPE ‘ND0T9 - FIRST DERIU OF MEAN MOTION, 4 OR - Pl.t <1S>‘ 

TYPE *19 - INCLINATION AT TINE 9, DEQRCES F7.4 <1S>* 

TYPE ‘N0DE9 - RIGHT ASCENSION OF ASCENDING NODE AT T9.DEGS.F7.4 <1S>* 
TYPE ‘E9 - ECCENTRICITY AT T9, FT. 7 <1S>‘ 

TYPE ‘0MEG9 - ARGUMENT OF PERIGEE AT T9, DE0REES.F7.4<15>‘ 

TYPE *M9 - MEAN ANOMALY, DEGREES, F7. 4 <1S>* 

TYPE ‘NO - MEAN MOTION, REUS/DAY. F10.8 <i*>‘ 

GO TO 39 
CALL EXIT 
END 


R 

TYPE KOUR 
C JHB 6/39/75 

COMPILER DOUBLE PRECISION 
SUBROUTINE KOUR 
DIMENSION C0RD(6) 

CALL APPEND(8,*C00RD*,3,IER) 

IFdER.NE.l )TYPE ‘OE‘,IER 
ACCEPT ‘HOW MANY STATIONS? ‘,IJK 
DO 19 I-l.IJK 

ACCEPT ‘LATITUDE! (D,M,S)» ‘.PD.PM.PS 
PHI-PD4 ( PM+ ( PS/69. ))/69. 

ACCEPT ‘WEST LONGITUDE (D.N.S ) I ‘.UD.UM.US 
UL-UD 4 (UM4(US/S9. ))/69. 

C0NU-6 . 2831 8S3972/369 « 

PHI-PNI8C0NU 


LONGITUDE 
SNSXES 3 




TVPE PREDICTLOAD 

PIDR PREDICT OT1 TCALC SGP SRU INCT SEPJI EXANM PEAK * 

GARB/L PORT . LB i DEUETE/’C GARB 
R 

TVPE PREDICT 
C TDB 4 DECTS 

COMPILER DOUBLE PRECISION 
DIMENSION TT(6).ID(3) 

CALL FGTINd, J.K ) 

TT(4)-I 
TT(S)-J 
TT(6 )-B 

CALL DATEdD.IER) 

TT< 1 )-ID(3) 

TT<2)-ID(1> 

TT<3)-ID<2> 

CALL OTl(TT) 

END 

R 

TYPE OT1 C 

C TDB 4 DECTS 

COMPILER DOUBLE PRECISION 
SUBROUTINE OTKT) 

DIMENSION T(S),TS(S> «9 

REAL IB. IM. NB, NOOTB, Mt.NODEB. MODEM. UI.HM. NDOTM, LLONGN* ) 
•LB.NODDT . J2, J3.MU, IS. NODES, LAMBDA E 

EXT E R N AL SEMI 25 

INTEGER VR M 

COMMON EPOCH. VR.MB.NODCB.OMEOB.NDOTB, 
tAM.EM. IM, NODEM. OMEQM, LM , {BJ. NDOTM. EM, NB.IB, LB, AS, C 

BELONG, LLONG, EXLNC, OMEGL . TMUEU . AMAG, ROOT , NOOOT , OMGDT . 
SUX,UV,UZ,RX,RV,RZ,RD0TX,RD0TY,RD0T2 

COMMON/IBU/DAYSINMO( 12 ) I) 

DATA DAVSINM0/31.,2S.,31..30..31.,30.,31..31.,30.,31., 
#30.. 31./ 

CuNU- .B1T453SSRS1 i DEGREES TO RADS 
ICOUNT-B 
J2-.BB1BB84B 
0 — 1 . 

JO-1 

MEM. 

IMANV-B 

X2-S. 

Z-S. 

X-S 

IU-1 

TU0PI-B.2B31BS3BT8 

CALL APPENDt IS, * TLOQ*, 3, IER ) BBS 

IP ( IER.NE. 1 )TYPE 'TLER'.IER 
CALL OPENtS, 'ELEMENTS*, 1, IER) 

CALL DFILUCPTAE'.IER) BBT 

CALL CP ILU( *PTAE* ,2, IER) 

CALL 0PEN(4, 'PTAE',3, IER) CHECK 

IF< IER.NE. DTYPE 'OE'.IER B9B 

C INPUT PROM NORAD 8BB 

READ(S)EPOCH, NOOTB, IB.NODEB.EB.OMEOB.MB.NB 
OMEGS-OMEGSBCONU 
IS-ISXCONU 


NODE0-NODEO«CONU 

M8-MO*CONU 

nb-nb* t uop: 

NDOTO-NDOTBtTUOPI 
CALL CLOSE (5. IER) 

IF (IER.NE. DTYPE 'CEMER 
LB-NB+NODEB+OMEGO 
AB-SEMI (EB.NO.IS) 

TEMP- 1 . S*NB*J2* ( A£/ ( AS* ( 1 . -EBIS2 ) ) )**2 
NODDT— TEMPXDCOSdB) 

OMGDT- TEMP* ( 2 . -2 . 5* ( DS IN (I B ) ) **2 ) 

BDT-3BB. 

SDT-1B. 

DT-BDT 

Xl-DT 

CALL 0PEN(8, 'COORD'. 1, IER) 

IF( IER.NE. DTYPE 'COE*. IER 

READ(S)LANBDA E.G1CSPhLG2CSPHI .SNPHI .CSPHI .G2SNPMI 
LAMBDA E- LAMBDA E-.B9S7* 

E.L. CORR TO COINCIDE WITH NORM) PREDICTION 

IPASS-1 

WRITE (IB, €9) 

URITCdS.SB) 

FORMAT ( 1 SATELLITE RISE TIME PASS B PEAK DURATIO 

DO BS 1-1,6 
TSd)-T(I) 

CALL TCALC(T,TSINCE, THETA. LAMBDA E) 
TSINCE-TSIMCE-CPOCH 

NO. OF DAYS SINCE NOST RECENT NORAD EPOCH 
CALL SGP(TSINCE) 

CALL SRU(TKETA,H,A,GlCSPHI,Q8CBPHI,fNPHI,CSPHI,G2SI#M 

IP(X2.NE.X1 >00 TO SS8 

IQO-B 

Z-B 

X-B 

IMANY-IHAHV+1 

JO-1 

A-B.B 

STEL-BB.BB jSTOU ANGLE 
URITE(4,339)A,STEL 
IFdMANY.NE. IPASS )GG TO 89B 
CALL CL0SEC4, IER) 

CALL CLOSE (i, IER) 

CALL CLOSE (IS. IER) 

CALL 0UERFL0U CS97 .S9B5. 'S' ) 

CALL CHAINCINTERP1 .SU'.IER) 

IFdCR.NE. DTYPE 'LSCE* 

CALL EXIT 

CALL CHAINCTRACK.su*. IER) 

IPdER.NE.DOO TO BBT 
TO SEE IF SATELLITE IS ABOUC THE HORIZON 

xr-b 

IFIH.LT. B. )00 TO 2SB 
IF(Q.QE.O. .OR.IFLAQ.EQ.l )00 TO IBB 
I PL AC- 1 
DO 1SB I-l.S 


is* 


19* 


135 


333 

339 


i 

-C' 

o 


987 


Ttn-rsti) 

DT-SDT 
H-Q 
X2-DT 
GO TO 2* 

Q-H 

IOUER-1 

CALL PEAK(IMANV.H,Z,X,ICOUNT,A.DT.IOUER> 
IFLAG-9 

14- T<4) 

15- T<5) 

16- T<6> 

IFCJO.EQ.nURITEUS. 136)14. IS. 16 
IFC JO. E0.1)URITEU9, 135)14.15, 16 
FowwT(ix.a<ia»f),!2,z) 
jo-a 

TSS-TSIHCE4CP0CH 

IF ( IU.EG . 1 )URITE< 4, 333 )TSS«DT 

FORMAT ( 1X.F14.9/1X.P6. 1 ) 

URITE<4,339>A.H 
F0RMATC2H A,F6.2,1HE,F6.2, ' 
iu-a 

DO 399 1-1.6 
TS( I >-T(I > 

CALL INCT(T.DT) 

GO TO 29 

XFCQ.GT .9. JDT-1DT 
XF(Xa.£0.9W0 TO 887 
IF<Q.GT.9)X2-DT 
O-H 
IU-1 

CO TO 299 
END 


59 


> 



TVJD-6. /THIS VERA'S JULIAN DATE 


199 


N*T<i)-l. | N- LAST VEAR 

DO 199 I-75.N i CHECK FOR LEAF TEARS THRU LAST VR 

IF ( I ''4X4. HE. I )06 TO 199 

XJD-XJD+l. 

CONTINUE 


N»T(2)-1. 

DO 59 i-l.fi iDAVS IN NOHTH THRU LAST RO 
TYJD-TVJD*DAy4iHMO<I ) 




299 

C 


CONTINUE 

N-TC1) ;N0U LOOK AT THIS YR 

IFCN/4X4.NE.N5G0 TO 299 jIS THIS NOT A LEAP V'R? 
IF C T f 2 ) .LE.2. )G0 TO 290 /ARE UE BEYOND 2/297 
ADD A LEAP DAY 


TYJD-TYJDfl. ; YES. 

TYJD-TYJD*TC3) 

JULIAN DATE AT INSTANT 

TSINCE-TYJD+T(4)/24.-MT(S)+T(6)/S0. )/1440. 

DT-T(4)*G0. 4 T(S)4T(6)/69. / (HR, WIN, SECS ) AS MINS 

TU- ( XJD+TVJD-2415929.0 )/3S525. 

THETA G9-DN0D( (99. 69998334 ( 36999 . 7689XTU )4 . 99938708S 
8TUXX2)»369. ) 

THETA G-DNOD< (THETA C94DTXDTHDT>,369. ) 

C SIDEREAL TINE IN RADIANS 

THETA- (DRODCC THETA C+LANIDA E ),369. ) )XTU0PI/369. 
RETURN 
END 
R 

TYRE SCR 

C SGR IV TDB — REV 2 14 75 AT 1439 
CONRILER DOUBLE PRECISION 
SUBROUTINE SOR(TSINCE) 

C THIS ROUTINE COARUTES SATELLITE POSITION USING A SINPLIF 

I ED 


QENERAL PERTURBATIONS METHOD, CLASSICAL MEAN ELENENTS ARC 


C 

E 

C 


INPUT/ AND POSITION, VELOCITY. 6 OSCULATING ELEMENTS AR 


RETAINED 


REAL 16, IH,N9,ND0T9,N6,N0DE6» MODEM, LH.NM.NDOTN.LLONC. 


TVRE TCALC 
C TDB 2^28/75 

COMPILER DOUBLE PRECISION 

SUBROUTINE TC ALCCT.TS I NCE, THETA, LAMBDA E) 

REAL LAMBDAS 
DIMEHSION T(6) 

COMMON /IBlK/DAYSir«0U2) 

EP-2442413.5 )1 JAN 75 

EPVR-75. 

TUOPI -6. 2831853972 

DTHDT- . 25968447 i DTHETA/DT 

XJD-EP4(T(1)-CPVR>*385. /ADD 366 DA/YR THRU LAST VR 


9ND0T6.L9. NODOT, J2,J3,NU, IS, NODES 
EXTERNAL EXAHM 


INTEGER VR 


INPUT PARAMETERS _ 

COMMON EPOCH, VR.M9. NODES, OMEG9.NDOT9, 
•AH,EH,IH t N0DEH / 0HE0H,LH,NH,ND0TH,E6»N9,I#,L9,A9, 
9€L5NG,LL««.ExtNG.OMEGL,TRUEU,RMAG,feot,NODDT,OMGDT, 
9UX. UV, UZ, RX. RV. RZ.RDOTX, ROOTY, RDOTZ 
J2-. 99198241 
J3— .999992562 
NDOT6-6. 

AE-1. 

MU-1. 

TUOPI -6. 2631663972 


COMPUTE TIME UARIANT MEAN ELEMENTS AT TSINCC 
TT-TSINCE iTIME SINCE EPOCH (DAYS) 

C DELETED '♦NDOT6*TTXX3' FROM THE NEXT LINE 

DM-N9XTT4ND0T9XTTXX2 jCHC IN MEAN ANOMALY 
DOMEG • OMGDTXTT ;D ARC PER 
DNODC-HODDTXTT iD ASC NODE 

LM-DM0D((L94DM4D0MEC4DN0DC)»TU0PI) /MEAN ORBITAL LH6T 
UDE 

OMEQM-DMOD ( ( OMEQ9+90MEC ) , TIIOPI ) »A.P. 


1 ,. 








O 


NOD€N-DHOD( (NODEO+DNODE^.TUOPI ) {RA OF AN 
IN-19 {INCLINATION UNCHNGD 
SINI-D$IN(IH) 

COSI-DCOS(IR) 

C DELETED '43.SND0T6STTSX3' FROR THE NEXT LINE 
NR-N048.XND0T0STT 
AN- A** ( ( ( N9/NN ) XX . 3333333333 )XX8 ) 

EN- 1 . -A9/ANI ( 1 . -E9 ) 

IF<EN>19. 19.89 

19 EN-9. 99991 


C 

c 

c 

c 

c 

c 


CONPUTE AND APPLY LONG PERIODIC TERNS CSUISRCPTD ’L’) C 

89 TENPL • (J3/J8)S(AE/Afl)XSINI/(i.-EnXX8) 

AXNL -ERXDCOS ( ONEGN ) 

AYNL-EHXDSIN < ONECN )- . 5XTEHPL R 

ELONG-DSQRT ( AXNLXX84AVNLXX8 ) TYPE 

ONEGL-DNOD< (DATAN8( AYNL. AXNL ) J.TUOPI ) j PRESERVE QUADC 
C LONG PERIODIC ON L IS* C 

ILONG-DRODC ( IN- .8SXTEHPLXAXNIX (3.+S. XCOSI )/ < 1 . 4C0SI ) t, 
iTWOPI) C 

C SOL* KEPLER'S EQUATION AND OTHER TWO-BODY FORNULAE 

C LONG PERIODIC ECC ANORi 

EXLNC-EXANH ( L LONG -ONEGL-NO DEN, ELONC ) 

C TRUE ARG OF LATITUDE* 

TRUEU-8 . SDATAM( DSORT ( ( 1 . +CLONG )/( 1 , -EL0N6 ) )S ( DSIN ( . 5S 
XEXLNC)/DCOS( .SXEXLNG) ) )+ONEQL 

NtlAG- AN* (l.-ELONG tDCOS ( EXLNC ) ) ,R SUB L 
C TRANSVERSE COHPOTCNT OF UCL VECTOR 

RVDT- DSORT (HU*ArtX(l."EL0NGX*8))*a./RNAG) 

C RADIAL CONP OF VEL VECTOR 

RHGDT- DSORT (RUSAN)SELONQ/RRACSDSIN(EXLNG) 

CONPUTE AND APPLY SHORT PERIODIC TERRS 

TENPS>.85XJ8X(AE/(Anxa.-£L0NGtt8>))SS8 j PERTURBATION 
SIN8U-DSIN(8.STRUEU) {CONSTANT 

C0S8U>DC0S< 8 . XTRUEU ) 

RHAG-RHA04TEHPSSSINIXS8XC0S8UX C ARX( 1 . -CL0HGSX8 ) ) 

TRUEU- DR OP < ( T RUEU- .SXTEHPSXC 6.-7. XS INI SX2 )XSIH8U).TU0C ) 
IS- IH+3. fTEHPSXSINIXCOSItCOSaU 
N0DES-N0DEH43. STENPSXCOSIXSIN8U 


UX--SINIXCNODE-COSUISNODEXCOSI 

U*— SIMjX$NODE4COSUSCNODEXCO$I 

VZ-COSUSSINI 

UX-SINIXSNODE 

UY— SINIXCNODE 

UZ-COSI 

RX-RI1AGXUX 

RV-RNAGXUV 

RZ-RNAGSUZ 

RDOTX-RNGDTXUX-RUDTIUX 
RDOTY-RHGDTXUY4RUDTXUY 
RDOTZ-RNCDTXUZ4RUDTXUZ 
ROOT ■ DSORT ( RDOTXSS84RDOTYSX84RDOTZSX8 ) 

RETURN 

END 

SRV 

SRU (SLANT RANGE VECTOR) TDB (S/87/75) 

CONSENTS AFTER LINES IN THIS SUBROUTINE 
ARE EQUATIONS IN APPENDIX OF ESCOIAL. 

’RETHODS OF ORBIT DCTERRINATION* 

CORPILER DOUBLE PRECISION 

SUBROUTINE *RV<THETA,H,A,G1CSPHI,G8CSPHI.SNPHI,CSPHI. 
9G8SNPHI) 

REAL IO,IH,NO,NDQTt,RO,NODEO,NODEH,LR # NR,NDOTH,LLONG, 

♦10, NODDT, LX, LY, LZ, LXH, LVH, L2H 
INTEGER YR 

CORRON EPOCH, YR, HO, NODE#, ONEGO. NDOTO, 

♦AN. EH, IR.NODCR,OREOR,LR,NR,NDOTR.Ef £ NB, 10. LO, AO, 

♦CLONQ , LLONG, EXLNG , OREGl , TRUEU , RRAG, RDOT , NO DDT , ORGDT , 
♦UX, UY, UZ .RX.RV , RZ , RDOTX , RDOTY , RD0T2 
I -S.SS3 1153078 

X-(G1CSPHI )SDCOS(THETA) ilA.68 
V- (G1CSPHI >*DSIN( THETA )| 1A.63 
RH0X-RX4X {lA.Ca 
RMOY-RY+V ilA.69 
RHOZ-RZ+CBSHPHI ilA.79 
RHOH-DSORT ( RH0XSX84RH0YXS84RH0ZXX8 ) 

LX-RHOX/RHOH {UNIT VECTOR FROM SITE 


r lA.71 

A. 73 TO SATELU 


CONPUTE QUANTITIES FOR OUTPUT 
SNODE-DS IN (MODES ! 

CNODE-DCOS(NODCS) 

SINI-DSIN(IS) 

COSI-DCOS(IS) 

SINU-DSIN< TRUEU) 

COSU-DCOS(TRUEU) 

C UNIT VECTOR POINTING TOWARD SATELLITE* 

C SEE P.194 IN ESC01AL, ’RETHODS OF ORBIT DETERRINATION*. 

C TO CHECK VALUES 

UX-COSUSCNODE -SINUS SNODESCOSI 
U Y -COSUSSNODE+S I NUXCNODEXCOS I 
UZ-SINUXSINI 


LY-RHOY/RHOH {DITTO 1A.73 
LZ-RMOZ/RHQHi DITTO 1A.74 
COSTH-DCOS( THETA) 

SINTH-DSIN( THETA) 

LXH-LXXSNPHIXCOSTHHYXSNPHIXSINTH-LZXCSPHI {1A.7S 
LVH— LXSSINTRUVXCOSTH i DITTO 

LZH-LXXCOSTHXCSPHI*LYXSInTHXCSPHI*LZXSNPHI {DITTO 
C0SH-DSQRT(1 .-IZHXS3 ) 

H-DATANOZH/COSH) ilA.TS IH ROOD 
CORRECTION FOR REFRACTION FOLLOWS {COURTESY OF RALPH PASS (GS 
FC) 

H-H+.9997SDC0S(H)/(DSIN(H)4D$QRT(.04-t(DSIN(H))SS8)) 
H-3S9.SH/TU0PI { ♦ OR — DEGREES FROR HORIZON 
A • DATANB ( L YH , -LXH ) ilA.77 

A-300.SA/TW0PI (DEGREES CU FROR NORTH 
IF(A.LT.B.)A-A43SB. {ADJUST COORDINATE SYSTERS 


1 i 

f J 


RETURN 

END 

R 

TYPE INCT 
C TDB 

COMPILER DOUBLE PRECISION 
SUBROUTINE INCT(T.DT) 
DIMENSION T(6) 
C0MM0N/IBLK/DAYSINM0C12) 
IF(DT.GE.6B. )GO TO 6BB 
T(6)-T(6)+DT j INCR SECONDS 
IF(T(6 KLT.6B. )G0 TO 3SB 
T(6)-T(6)-6B. {RESET SECONDS 
CM IF(DT.GE.6B. )T(5) a T!S)*DT/SS. 


IF(DT .GE.6B. )GO TO TBB 
T(5)«T(S)+1 {INCR MINUTES 
TBB IF!T(S).LT.6B. )GO TO 3SB 

T(5> a T(S)-«0. i RESET MINUTES 
T!4)-T!4H1. {INCR HOUR 
IF(T(4 ) .LT .24. ) GO TO 3SB 
T(4)»T(4)-24 .RESET HRS 
I-T(2) jPTR TO NO 
TOl-TOHl. {INCR DRY 
IYR-T(l) 

ILEAP-B 

IF(I.E0.2.AND.IYR/4«4.EQ.IYR)ILEAP a l 
DAYSIf*W(2) a 2B*ILEAP 
IF(T<3).LE.DRVSINB0(I ) XJO TO 35B 
T(3)»l. {RESET DRVS 
DAYSINN0!2) a 2S. {RESET FEB 
T(2)-T(2H1. {INCR NO 
IF!TC2).LE.I8.)Q0 TO 35* 

T(2)-l. {RESET NO 
TUl-TUHl. {INCR YR 
35# RETURN 
END 
R 

TYRE SENI 
C TD» 2/2B/7S 

COMPILER DOUBLE PRECISION 
DOUBLE PRECISION FUNCTION SEMI!EE,XNN,XII ) 
COMPUTES THE MERN (KOZRI) S EM I -NR J OR AXIS OF A SATELLITE 
YY-. 3333333333 
XJ2-.SB1SB241 
XNU-U467.2S29B 
AA a (XNU/XNNW2 )«YY 

DD— 1 .S*XJ2*( ( 1 ./RR)*t2)/< <DSORT( 1 .-EEW2) )«3) 
DD a DDX( l.-l .SKDSINCXII ) )BX2 ) 

SENI a AAI U . ♦YVBDD-YYBDDBB2 ) 

RETURN 

END 

R 

TYPE EXANN 
C TDB 2'28/ 7S 

COMPILER DOUBLE PRECISION 

DOUBLE PRECISION FUNCTION EXANM(XMM,ECC) 

COMPUTES ECCENTRIC RNONRLY USING KEPLER'S EQUATION 
TU0PI a 6.2B31BS3472 



EXANN-DMOD ( XMP1, TUOPI ) 

DO IB I-l. SB 
AA-ECCIDSIN (EXANN) 

DELN a XHM-EXANM+AA 
ZZ« 1 . -ECCXDCOS ( EXANN ) 

D£LE a DELM/(22*( ( .S«DELN )/Z2)*RR ) 

IF ( DABS ( DELE )-l . )3B,3B.2S 
2B DELE-DELE/DABS ( DELE ) 

3# EXANN a EXANM+DELE 

IF ( DABS ( DELE )- .BBBBAi )4B, IB, 1# 

IB CONTINUE 

4B CONTINUE 

return 

END 

R 

TYPE PEAK 
C JHB 6/3B/7S 

COMPILER DOUBLE PRECISION 

SUBROUTINE PEAK!IRAMY,H.Z,X,ICOUNT.A,DT.IO!*R> 
CALL APPEND! IS, *TL0G*,3, IER ) 

ICOUNT-ICOUNTH 
IF(H.LT.Z )GO TO 333 
Z-H 
RETURN 

333 IFIX.EO.l )GO TO 1 

IM-INANY+l 
ZNUN*IC0UNT-1 
ICOUNT-B 

XMUN a ( ZNUNB2BDT )/SB. 

IFIIOUER.EO.BICO TO IB 

IFCZ.EQ.B .0R.X.EQ.B)URITE(IS,222>IN.Z,XNUN 
IB IFIZ.EQ.B .0R.X.E0.B)URITE(1B,222)IM,Z,XNUN 

an F0RMAT(11X,I3,SX,F7.2,2X,F5.B//'/) 

IF(Z.LT.1S.AND.A.LT.B)CALL FCHAN< •UAIT2B.SU* ) 

1 CONTINUE 

CALL CLOSE! IS, IER) 

RETURN 

END 

R 



TYPE LS8DIS 


1 46 ST 

JOHN R.NINEMILE8R 




129*. 


e. 


2 46 

6*3 

43 

84 

147 

243 

375 

546 

•>S6 

1*6* 

3 46 

1426 

1816 

2246 

2781 

3296 

3954 

4766 

5548 

642* 

4 46 

7386 

8416 

9S26 

1474* 

11966 

13444 

14268 

15488 

1669* 

5 46 

17866 

19266 

26666 

22144 

23666 

25244 

26766 

28386 

2990* 

€ 46 

31666 

33266 

35666 

35644 

36266 

36844 

37466 

38666 

3860* 

X 47 ST 

JOHN R. 

DICKEY 




27*4. 


3.5 


2 47 

6.5 

166 

566 

1144 

2626 

312* 

4456 

6656 

78** 

3 47 

9866 

11866 

13866 

16144 

18666 

2144* 

23666 

26466 

2970* 

4 47 

32966 

36266 

39466 

42744 

46666 

54144 

54266 

58766 

632*0 

5 47 

6 47 

1 48 ST 

67766 

72266 

76766 

81244 

85766 

9484* 

94766 

99266 


JOHN R # 

FORT KENT 




5694. 


6.5 


2 48 

1.6 

456 

1286 

228* 

3666 

5314 

7496 

19166 

13*44 

3 48 

16666 

19766 

23766 

28444 

33666 

38444 

44666 

56666 

56544 

4 48 

5 48 

63566 

71666 

79666 

87444 


144444 

113666 

123666 

133444 

6 48 

1 49 PEH08SC0T. 

U. ENFIELD 




657*. 


1. 


2 49 

6.5 

1866 

2666 

3594 

4534 

5754 

6966 

8276 

9744 

3 49 

11266 

12966 

14666 

146*4 

1854* 

847*4 

23666 

25566 

29*44 

4 49 

36766 

33666 

36566 

395*4 

4SS44 

45544 

49666 

52566 

5684* 

5 49 

54444 

63966 

67666 

7184* 

7594* 

84444 

64666 

88266 

98544 

6 49 

97444 

161566 

166266 

111*44 

116444 

181444 

126666 

131666 

136544 

1 56 CARMASSCT 

R, N. ANSON 



354. 


2. 


8 54 

4.5 

4 

54 

188 

874 

494 

795 

1186 

1554 

3 54 

8*34 

8594 

3834 

3884 

4534 

5874 

6676 

6876 

7674 


181*4 


4 5# 

5 5# 17*44 

6 54 272*4 

1 SI SACO R, CORNISH 
8 51 *.5 815 

3 51 384* 4*4* 

4 51 129** 14144 

5 51 845** 85*44 

6 51 367** 3*144 

1 38 PLYHOUTH 


14344 

198*4 


118*4 

843*4 

31**« 


1884* 

81444 

383*4 

751 


1794* 

8994* 

4244* 


*8 38 

3 38 

4 38 

5 38 

6 38 

1 9 

8 9 

3 9 

4 9 

5 9 


.5 
385* 
94*4 


13144 

88644 

189*. 

1464 

8134 

19844 

318*4 

43*44 

688.4 


141*4 

83744 


15444 

84*44 


If 
8644* 


351 531 

4954 5954 

15444 166*4 

87844 89644 

39644 44944 

18 8 1 

5* 19* 439 

375* 489* 486* 548* 6164 

1441* 11364 18354.4 13354.* 14354.* 

1857*. • 1967*. • 8477*.* 8194*. • 83454.* 8484*. • 
891*4.4 34354.4 31654.4 3896*. • 3431*.* 3571*. • 
GOFFS FALLS 4837 3498.* 

1.* 115.* 454.* 111*.* 838*. • 4784.* 

88444.4 86844.4 31644.* 36440.* 4184*. • 4644*. • 
6584*.* 7*44*. 4 74844.4 794**.* 8444*. • 8884*.* 


1444 

9364 


6 9 

153666.6 

6.6 

6.6 

1 52 

IPSUICH R. 
6.2 

IPSUICH 


2 52 

6 

6.76 

3 52 

174 

246 

325 

4 52 

1156 

1246 

1346 

5 52 

2685 

2175 

2296 


1. 7 

193* 8S54 

14644 11744 

8184* 83844 

3444* 35344 

45344 46*4* 48344 

4.* 11.* -1.0 18.6 

1VN 834* 879* 

69*4 7714 8584 

1S37*.* 16489.4 1747*. • 
85444.* 866**.* 87*5*. • 
3711* 3851* 3991* 

*.* 18. 1.* 14. 

88**.* 185**.* 1784*.* 
5484*. • 5564*. 4 6*444.* 
936*4.* 9644*.*1*38**.* 
*143*4*. 4148*4*. • 


58 

53 N. NASHUA R, FITCHBURG 


53 

53 


.8 

416 


17 

53* 


6.6 

6.6 

6.6 

6.6 

6.6 

6.6 



124. 


2.6 


3.39 

17.5 

32 

52 

82 

122 

436 

556 

885 

825 

956 

1655 

1446 

1546 

1646 

1746 

1845 

1955 

2416 

2536 

2876 

2816 





63.6 


2.1 


52 

96 

126 

175 

242 

323 

776 

896 

1616 

1136 

1266 

1466 


31 

654 


> 

1 

Ln 

UJ 


4 53 

154® 

1699 

1189 

1989 

2149 





5 53 

6 53 










1 55 

BRANCH R* FORESTDALE 




91.2 


1.61 


2 55 

1.4 

6 

26 

65 

122 

211 

299 

421 

561 

3 55 

721 

811 

1141 

1215 

1395 

1511 

1781 

1985 

2215 

4 55 

5 55 

6 55 

2425 

2551 

2891 







1 56 PAUTUXET R. 

CRANSTON 




211 • 


3.2 


2 56 

1.2 

19 

47 

89 

141 

217 

287 

368 

451 

3 56 

531 

615 

675 

741 

m 

861 

925 

995 

1165 

4 56 

1135 

1215 

127S 

1345 

1415 

1485 

1555 

1625 

1695 

5 56 

1765 

1835 

1911 

1991 

2171 

2151 

2231 

2311 

2391 

6 56 

2471 

2551 

2631 

2715 

2815 

2895 

2985 

3175 

3165 

l 19 HARTFORD 


4 19 

1 


11489 

9.9 16.9 9.9 

22.1 

2 19 

1.9 

6.9 2HM 

4m. • 

7m. 1 

iNN.e i« 

29*9.9 169*9.9 2* 

m.i 


19 

19 

19 

19 


45 CONN* 
45 
45 
45 


23M1.1 28HM 32111.1 37MM 42M.I 471H.1 S4MM Mttt.f 66M1.1 
?!•§•*• NNM HUM 95M0*tl96#M*tll4tt6'#124#M*tl34Mi*tl43M#*i 
159tM*#152#M*#173#M*ilSMM.0192Mi*08i39M*iC14Mt*MtMM*i237#M*t 

lift MCMAA AantAAA a a a a a a a a a a a X A 

'••♦•Ca riWtVOOVff V»V v»v V*V V»f VtV 

R, RIDDIETOUN II 


5 45 

6 45 

1 57 

2 57 

3 57 

4 57 

5 57 

6 57 
1 51 
a si 

3 51 

4 51 

5 58 

6 51 

1 59 

2 59 

3 59 

4 59 

5 59 

6 59 
1 61 
2 61 

3 61 

4 61 

5 61 

6 61 


WINCHESTER* CONN. 


ST. FRANCIS R. » N.l* 


SHIELDS BROOIC* NE* 


1. 


ALIACASH RIUER.NE. 
1.5 16 

5171 6151 

11211 21311 


263 585 

7261 8541 

25311 


1171 

INN 

27m 


1251. 

1691 

11711 

31311 


1.5 

2471 3321 

tarn um 


4m 

um 


"T"*" 




t ^utnuv 1 . ' x \ 




U1 

-O 


TYPE 

El 

7et2 

7207 

7110 

7273 

7071 

7272 
7356 
712? 
7261 
7021 

7103 

7104 

7105 
7125 
7246 
6063 

7271 
7345 
7254 

7355 

7206 
0 

TYPE 

7012 

7207 
7110 

7273 
7071 

7272 

7356 
7127 
7201 
7021 

7103 

7104 

7105 
7125 
7246 
6063 

7271 
7345 
7254 

7355 
7206 
R 

TYPE 

18 

7012 

7777 

7273 
7071 

7272 

7356 
7201 
7021 
6063 


IMDX 

2 NIMENXLE iR., ME. 

2 ALLAGASH RIUER, HE. 

1 ALLAGASH PALLS, HAINE 

2 FORT KENT, HE. 

2 UEST ENFIELD, HE * 

2 NORTH ANSON, RE. 

2 CORNISH, HE. 

3 SOUTH NT., HUH* 

2 PLYHOUTH, * N.H. 

2 GOFF S FALLS, N.H. 

1 CRREL, NORTH DAKOTA 
1 CRREL, NORTH DAKOTA 
1 CRREL, NORTH DAKOTA 

1 CRREL-SUGARLOAF, HE. 

3 UACHUSETT HT., HA. 

2 IPSWICH, HA. 

2 FITCHiURC, HA. 

2 CRANSTON, R. I. 

2 HARTFORD, CT. 

2 HIDDLETCUH, CT. 

2 HANCHESTER, CT. 

DIRC 


QDIRC 


^271 

7777 

7345 

7254 

7355 

7206 
6504 
7777 

7207 
R 


TYPE STARTERS 


0*0 


0*0 


A 0.00E 

60.00* 

A330.00C 

40.00: 

A300.00E 

20.00* 

A270.00E 

10.00* 

A240.00E 

8.00* 

A210.00E 

6.00* 

A180*00E 

4.00* 

A150.00E 

3.00* 

A120.0O6 

0.00* 

A 90.006 

0.00* 

A 66*006 

0.00* 

A 30.006 

0.00* 

A 0.006 

0.001 

A330.006 

0.00* 

A300.006 

0.00* 

A270.00E 

0.001 

A240.00E 

0.001 

A210.00E 

0.00* 

A100.00E 

0.00* 

A150.006 

0.00* 

A 120. 006 

0.00* 

A 90.006 

0.00* 

A120.00E 

0.00* 

A150.006 

0.00* 

A1M.00E 

0.00* 

A210.00E 

0.00* 

A2 40.006 

0.00* 

A270.00E 

0.00* 

A 300. 006 

0.00* 

A330.006 

0.00* 

A 0.00E 

0.00* 


A999.0iE999.t0t 

R 


xarntojWMt, 




TYPE DTABST 

C TDB 3t JULY ?S PROGRAM TO CONCERT DEC I HAL AZIMUTH AND 
C ELEVATION ANCLES TO BCD SUITABLE FOR INPUT TO INTERFACE 
C TO 1848 DIGITAL COMPARATOR 

COMPILER DOUBLE PRECISION 
DIMENSION X(2 ), IBCDHC8 ), IBCDL (2 ) 

COMMON/IBU/'BCDHIGH ( 18 ), BCDLOW(B ) 

DATA BCDHIGH/ 2 tt.,ltt..it., 4 i.,at.,lt.. 8 ., 4 ., 2 .,l./ 
DATA BCDL 0 W>\B,. 4 ,. 2 .. 1 ,.t 8 ,.t 4 ..ta..tl/ ;BCD TABLE 
CALL OPEN ( S.* STARTERS \ 1 . IER) 

CALL OPEN ( 6 » * START ANGLE* ,3 . IER ) 

IF( IER.NE. 1 )TYPE ‘OE’.IER 
Ml — 1 

5 R£AD(S.lt.END a 5 tt )TSINC£.DT.X jREAD 1 ST RECORD ON FILE 

It F 0 RNAT(F 14 . 9 /F 5 . 1 / 2 ( 1 X,F«. 2 ) ) 

GO TO IS 

WRITE(6,11 )TSINCE»DT 
11 F0RMAT(F14.9«FS.l ) 




IS 


> 

i 

un 

Ln 


2# 

3 t 


St 

ltt 



lit 

12t 

Stt 

R 


DO ltt 1-1,2 |AZ .EL 

IHIGH-t 

ILOU-t 

DO 3t J-l.lt 

IFCBCDNIQHC J ) . LE.X( I ) )GO TO 2t , CHECK BCD TABLE 
GO TO 3t ; SMALLER THAU VALUE IN TABLE 

X( I )-X(I )-BCDHIQH( J ) «>-UALUE IN TABU 

II-J+B-(J-I )*2 iNAPPINO FROM DO LOOM INDEX TO BIT POSITION 

CALL ISETC IHIOH. II ) jBIT ORDER ON P.9-11 OF FORT IU 

CONTINUE 

DO St J-1,8 

IF(8CDL0U(J).U.X(I ))CO TO 4t 
GO TO St 

X<I)-X(I)-BCDLOU<J) 

II-J>€-(J-1 )*2 .BIT NAPPING 
CALL ISETCILOU.il) 

CONTINUE 
IBCDH(I>-IHIQN 
IBCDLCI )*ILOU 
CONTINUE 

WRITE BINARY ( 6 ) IBCDL ( 1 ) , IBCDH (1 ). IBCDL ( 2 ) , IBCDH ( 2 ) 

READC5,llt)X 

FORMAT* 2( 1X.F6.2) ) 

IF ( X ( 1 ) . EG . 999 . )GO TO 12t 
GO TO IS 

WRITE BINARV(t)Nl 
GO TO 5 

CALL CLOSE(S.IER) 

CALL CLOSE (6. IER) 

IF(IER.NE.1)TYPE ’CEMER 

STOP 

END 



Sample output from P3 
A-56 


TOTAL '«UH»€R Of HESSAGES 
NUN 8 ER OF NED NESSAGES • 
PlATFORfl LOCATION 

1 . 0 . t 

7171 NINENILE W.. HE. 

8x3# EST 
3X33 EST 

7384 SHIELDS 9ROOK, HE. 

8*31 EST 
8*34 EST 
8*37 EST 

7373 PORT KENT, HE. 

8*29 EST 
8*3S EST 
8*38 EST 

7147 ALIACASH FALLS. HE. 

8*3* EST U€$— 8.4S1N 

8*33 EST U£S»-#.4S1N 


<4 


7325 NED.LALTw.H, «m>. 
7*87 NED. UAL*- AN, "A. 


5/2«/77 
STQ- a.asrr 

STG- 2.2SFT 
S/2S/7? 
STG- 7 . 9 «rr 
STG- 7.98FT 
STG- 7.98FT 
S/26/7? 
STG- 7.78FT 
STS- 7.78FT 
STS- 7.78FT 
5/26/77 


1248. CFS 
1848. cfs 

I.CfS 
8 . CFS 
8 . CFS 

136*8. CFS 
136*8. CFS 
1368*. CFS 


I.* csh 

1.8 CSN 

8.8 csn 

8.8 CSN 
8.8 CSN 

2.4 CSN 
2.4 CSN 
2.4 CSN 


?#7i 

WEST ENFIELD, HE. 

5/26/77 





8 1 27 EST 

STG- S.15FT 

11718. CFS 

1.8 

CSN 


8*31 EST 

STG- S.15FT 

11718. CFS 

1.8 

CSN 


8*34 EST 

STG- S.15FT 

1171 8. CFS 

1.8 

CSN 


8137 EST 

STG- S.15FT 

11718. CFS 

1.8 

CSN 

7272 

NOGTH ANSON, HE. 

5/26/T? 





8*32 EST 

STG- 3.S5FT 

292. CFS 

l.l 

CSN 

7356 

COGNISH, HE. 

S/26/77 





8*3# EST 

STC- 3. 75FT 

1615. CFS 

1.3 CSN 


8i33 EST 

STG- 3.75FT 

IttS.CFS 

1.3 

CSN 


8:36 EST 

STG- 3.75FT 

1685. CFS 

1.3 

CSN 

7127 

SOUTH HT. , N.H. 

S/26/77 





8*29 EST PGC- 

8.33IN 





8132 EST PGC- 

8.33IN 





81 35 EST PGC- 

8.33IN 





8*38 EST HGC- 

8.33IN 





7228 FORESTDAIE, B-I. 
734S CRANSTON, R.I. 
8*3# EST 
8*36 EST 
8*39 EST 


7254 HARTFORD. CT. 
8*32 EST 
1*35 EST 
I 1 38 EST 

7242 HIDDLETOUN, CT. 
7286 WINCHESTER, CT. 
8*32 EST 
1*35 EST 


72*1 PLYNOUTH, N.H. 5/26/77 

8*3# EST STG- 8.58FT 

8*33 EST STG- 8.58FT 

72*7 GOFFS FALLS, N.H. NO REPORT 

7183 CRREL, NORTH DAKOTA NO REPORT 

7184 CRREL, NORTH DAKOTA NO REPORT 

7185 CRREL. NORTH DAKOTA NO REPORT 

7125 CRREL-SUGARLOAF, HE. NO REPORT 


518. CFS 8.8 CSN 
S18.CFS 8.8 CSN 


7246 UACHUSETT NT.. NA. 

8*34 EST PRC- 8 . 
6463 IPSUICH, NA. 

8*31 EST 
8*36 EST 

7271 FITCH 8 URG, HA. 

8:32 EST 
8*35 EST 
8*38 EST 

7142 CHICOPEE. HA. 

7821 WESTFIELD. HA. 


5/26/7? 

82 IN 

5/26/77 
ST8- 3.82FT 
STG- 3.92FT 
5/26/77 
STC- 6.86FT 
STC- i.##rr 
sto- 6.8#rr 

NO REPORT 
NO REPORT 


291. CFS 2.3 CSN 
291. CFS 2.3 CSN 


19*#. CFS 29.9 CSN 
19##. CFS 29.9 CSN 
19M.CFS 29.9 CSN 


NO REPORT 
NO REPORT 


NO REPORT 
5/26/77 
STG- 3.22FT 
STG- 3.22FT 
STG- 3.22FT 


22. CFS 8.1 CSN 
22. CFS 8.1 CSH 
22. CFS 8.1 CSH 


5/26/77 




STQ. 

5.63FT 

1452#. CFS 

1.4 

CSH 

fTQ® 

5.63FT 

14528. CFS 

1.4 

CSH 

STG- 

5.63FT 

1.528. CFS 

1.4 

CSH 

NO 

GEPOGT 



5/26/77 




STG- 

2.UFT 

• .CFS 

M 

CSH 

STG- 

2.11FT 

• .CFS 

l.l CSN 
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CIRCUIT DESCRIPTION OF THERMOCOUPLE INTERFACE 
Thermocouple Interface 

The thermocouple unit was designed to accept inputs from up 
to 112 copper-constantan thermocouples which are arranged in 16 
"banks" of 7 inputs. Each bank of inputs is recorded with an 
individual update of the DCP. There were 28 thermocouples used 
at the Sugarloaf Mountain installation and, as a result, four 
banks were used. 

The reference thermocouple junction is compensated to within 
+0.3°C over a temperature range of -50° to +40°C using a combin- 
ation of five thermistors. The normal measurement temperature 
range is -34° to +32°C with a resolution of +0.25°C or 10 mv. It 
is possible to trade range for resolution, or vice versa, by 
selecting different groups of data bits from the analog-to-digital 
converter. 

The thermocouple unit is designed to operate with a LaBarge 
Electronics Convertible Data Collection Platform (C/DCP) which has 
a memory capability. Power requirements are +12 volts (nominal) 
at 0.5 ampere and 5 volts at 1.0 mllliampera. The 12-volt power 
is applied only while the C/DCP is acquiring new data. 

Circuit Description 

The electronic components are arranged on nine circuit cards 
that are identified as follows (figure B-l): 

Input multiplexer, group I (MX I) 

Input multiplexer, group II (MX II) 

Input multiplexer, group III (MX III) 

Input multiplexer, group IV (MX IV) 

Amplifier and channel multiplexer (AMPL CHANNEL MUX) 
Analog-to-digital converter (A/D CONV) 

Power supply and reference junction compensator (P.S. 

REF JCT COMP) 

Latch (LATCH) 

Interface (INTERFACE) 

Thermocouple signals are routed through the input multiplexer, 
consisting of cards MX I, MX II, MX III and MX IV (figure B-2 through 
B-5). Seven input signals are read during one DCP update sequence. 
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The 16 banks are numbered 0 through 15. When an update occurs, 
the next highest numbered bank is read. Thus, after 16 updates a 
total of 112 inputs have been sampled. 

Each of the seven inputs of the selected bank occupies an 8- 
bit word of DCP data and are designated as channels 1 through 7 . 
The remaining word, corresponding to channel 0, is the number of 
the bank selected. Four of the eight bits of this word are re- 
quired for this process and the remaining four are always zero. 
There is an available analog input for channel 0. It is used only 
during bench testing and adjusting and has no corresponding DCP 
data word. 

The amplif ier/channel multiplexer, consisting of one card 
(AMPL/CHANNEL MUX), scans the seven channels and amplifies their 
signals to a suitable level for digitizing (figure B-6) . This 
module also controls the associated analog-to-digital converter 
(A/D CONV) and provides signals (called latch strobes) which con- 
trol the transfer of the digitized data to the latch card (figure 
B-7) . 


The A/D CONV is a prepackaged unit mounted on a separate card 
(B-7). Although it has a 12-bit-plus-sign capability, only 7 
bits and the sign appear in the DCP data. The A/D CONV has a 
differential input, in addition to the normal analog signal input, 
and receives the output of the reference junction compensator. 

The reference junction compensator and system power supply 
modules (P.S./REF JCT COMP) are mounted on one card (figure B-8) . 
The reference junction compensator generates a voltage which 
matches that produced by the reference junction itself over the 
temperature range of -50° to ■f40°C. The compensating voltage is 
introduced to the A/D CONV at the same level as the amplifier out- 
put. The compensating voltage is the reference junction voltage 
multiplied by the gain of the amplifier. 

The latch card (LATCH) contains eight 8-bit latches to accept 
parallel data (figure B-9) . A strobe signal is required to enter 
data and an enable signal is needed to make data available at the 
output. If the appropriate signal is not present, the device 
appears to be an open circuit. Thus, the input and output ports 
can operate on common busses. The 8-bit parallel data from the 
A/D CONV is bussed to all of the latch inputs. A strobe signal 
from the channel multiplexer operates the input of each latch 
in succession, according to the analog channel being sampled. 

The data input to the C/DCP is by way of a 16-bit parallel 
bus. The 8-bit latches are paired, so that the output lines of 
latch pairs 0-1, 2-3, 4-5, and 6-7. Each form a 16-bit data word. 


Amplifier and channel multiplexer (AMPl/CHAMEL MUX) 
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Figure B-8 Power supply and reference junction compensator (P.S./REF JCT COMP). 
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The enable signals for all latches (except O^and 1) are re- 
ceived from the C/DCP. The enable signals are the triggering 
signals that would normally be used to start the punch motors in 
the analog-to-digital recorders (ADR’s) operated by the U.S. 
Geological Survey. These are labeled ’’Punch Motor No. 1”, ’’Punch 
Motor No. 2”, etc., in all instructions and diagrams. Punch Motor 
No. 1 initiates a scan by the channel multiplexer. Punch Motor 
signals 2, 3 and 4 enable latch pairs 2-3, 4-5, and 6-7, respec- 
tively, to be started. Since Punch Motor No. 1 arrives before 
latches 0 and 1 can be supplied with data, the input strobe sig- 
nal for latch 7 also serves as the output enable signal for 
latches 0 and 1. 

Latches 1 through 7 carry thermocouple data, which consists 
of 7 magnitude bits and one sign bit. Latch 0 carries a 4-bit 
word to indicate which of the 16 banks of inputs is being scanned. 
The remaining 4 bits in Latch 0 are not presently utilized. 

The 16-bit output from the latches goes through a transistor 
switch array (INTERFACE) to the C/DCP (figure B-10) . The transis- 
tors invert the polarity of the data signals and also simulate 
the grounding-type signal that would come from an ADR. 

The primary power supply is 12 volts DC, supplied from the 
battery that powers the C/DCP, and is turned on during an update 
sequence by a semiconductor switch in the C/DCP. The switched 
12V supply powers a +15V and a +5V supply which are electrically 
isolated from each another to avoid current loops. The bank select 
counter and certain parts of the control portion of the amplifier/ 
channel multiplexer require +5V continuously. The voltage comes 
from a regulated supply within the C/DCP. 

Operating Sequence 

The operating sequence is as follows: 

1. The C/DCP begins an update sequence, by applying a 
12V input to the +15V and 4-5V power supplies, and initializing a 
’’Punch Motor No. 1” trigger pulse to the control portion of the 
amplifier /channel multiplexer. This resets the channel multi- 
plexer to channel 0. At the same time, the channel multiplexer 
sends a ' clock pulse to the Bank Select Counter, which is on MX I, 
advancing it by the count of one. 
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Fifura B-10 Interfaca (IKTERFACE). 
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2. The channel multiplexer scans the seven channels of 
the bank selected. It remains on each channel for about 0.5 
seconds, during which time a pulse is sent to the A/D CONV to start 
a conversion. When the conversion is complete, a pulse (latch 
strobe) is sent to the proper latch to enter the data. This 
process is repeated for all channels. On Channel 0, a conversion 
is performed, but channel 0 data is not entered into a latch. 
Instead, the state of the bank select counter is entered. Thus, 
upon completion of a scan, Latch 0 contains the bank number, 
while Latches 1 through 7 contain temperature data. 

3. The data in latches 0 through 9 is loaded into the 
C/DCP memory. As previously described, the C/DCP "Punch Motor" 
signals, together with the strobe signal for latch 7, enable 
the latches, two at a time. 

4. Upon completion of the update sequence, which re- 
quires about 90 seconds, the C/DCP removes the 12-volt power from 
the unit. The C/DCP continues to supply 5-volt power to the bank 
select counter and certain other circuits that must remain ener- 
gized. 


